:arXiv:cond-mat/0404049vl [cond-mat.mes-hall] 2 Apr 2004; 




Quantum spin dynamics in 
single-molecule magnets 



Proefschrift 



ter verkrijging van 
de graad van Doctor aan de Universiteit Leiden, 
op gezag van de Rector Magnificus Dr. D. D. Breimer, 
hoogleraar in de faculteit der Wiskunde en 
Natuurwetenschappen en die der Geneeskunde, 
volgens besluit van het College voor Promoties 
te verdedigen op donderdag 18 maart 2004 
te klokke 15:15 uur. 



door 



Andrea Morello 
geboren te Pinerolo (Italie) in 1972 



Promotiecommissie 



Promotor: 
Referent: 
Overige leden: 



Prof. dr. L. J. de Jongh 

Prof. dr. G. Aeppli (University College London, Groot Britannie) 

Dr. H. B. Brom 

Prof. dr. G. Frossati 
Prof. dr. P. H. Kes 
Prof. dr. M. Orrit 
Prof. dr. J. Zaanen 



ISBN: 90-77595-20-1 

Printed by Optima Grafisclie Communicatie, Rotterdam: www.ogc.nl 

This work is part of the research program of the Stiching voor Fundamenteel 
Onderzoek der Materie (FOM), which is supported by the Nederlandse Organ- 
isatie voor Wetenschappelijk Onderzoek (NWO). 



Cover: a sequence of nuclear spin echo traces obtained at increasing time intervals after 
the inversion of the ^^Mn nuclear polarization in Mni2-ac. The recovery of the equilibrium 
polarization is driven by quantum tunneling of the cluster spin. 



Said the straight man 
to the late man: 
"Where have you been?" 
"I've been here and 
I've been there and 
I've been ... in between." 

King Crimson, 

"I talk to the wind", 1969. 



Contents 



Introductionl 



2 Exp erimental tec hniauesi 5 
2^^^ ^^^^^^^^^^^^i] ■ . ■ ■ 5 



2.2 ^^clear Magnetic Resonanca 
2.2.1 Basics of pulse NMEI . 



2.2.2 Instruments and circuitrvl 11 

2.2.3 Heating effects! 14 



2J ac 



s^cc^tomctcj . 14 



2.4 SQUID magnetometer 16 



2.4.1 Working princi pl^ 16 



2.4.2 



Design and construction of the pick-up circuiti 23 

2.4.3 ^rtica^noveme'nj _ 24 



2.4.4 Grounding and shielding! 27 

2.4.5 Automation! 28 



2.5 Torque magnetometer 28 



2^5J_^Prin£i 2]gs_of_caji^^ 28 

2'5'^^^esi^^m^^on^fauctionF 30 

2'5'^^^erfonnanc^ 33 



3 Theoretical aspects of molecular maenetism! 


35 


3.1 


Macroscopic Quantum phenomenal 


35 




3.1.1 Extensive differenc^ 


36 




3.1.2 Disconnectivitv 


36 




3.1.3 Macroscopicitv of a single-molecule magnet 


37 


3.2 


Effective Spin Haniiltonian 


37 




3.2.1 Superparamagnetic blocking! 


38 




3.2.2 Spin Quantum tunneling! 


... 39 




3.2.3 Coupling with external fields: coherent and incoherent 


tunneling! 




3.2.4 Parameters for ]\^^7-acj 


44 




3.2.5 Parameters for Mn^ 


46 


3.3 


Hvperfine coupling^ 


... 47 



II 



CONTENTS 



Hvperfine fields on the ^^Mn nuclei in Mni9-ad 48 

3.3.2 NMR spectra in perpendicular fieldl .... 49 

51 



3.5 



N3J£leai_spiiL.dxnainicj 54 



^j^^^^^iclca^gi^a^c^rel^catio^^J^L^I 54 



^^^^^^^^nsvc^^^ir ^Din relaxation fTSSR,)! 56 



3.6 



ProkoTev^tem^^eo^ 57 



3.6.1 Central spin 4- spin bathi 58 



3.6.2 Tunneling rate! 63 



3.6.3 Comparison with the Landau-Zener formalisnJ 64 

3.6.4 Apphcation to SMMsl 65 



^^^^easi^mcn^^n^^a^^nalvsisl 



4.2 



67 

. . . 67 

Nuclear relaxation in moderate parallel fieldd 70 

14.2.1 NST.R, and TSSR, in zero fieldl . . . 70 



4.2.2 Field dependencies 76 



4.2.3 Deuterated sam ple| 81 



4.2.4 Comparison with a Mn'^+ sitel 82 



Effects of strong perpendicular fieldsl 84 



4.3.1 NME, spectral 85 



4.4 



4.3.2 NST.R in perpendicular field j 89 

Ni^c^l^^r ^pin l^^mp^r^^ ^rq 92 



^^^^^Tim^^vohitio^^^tii^iucka^gii 93 



4.4.2 Effect of large magnetic fields and superradiancel 94 

Proposg^fo^^mi^dn^theorvl 98 



45^^^Basi^assumptiona 



4.5.2 Tunneling mechanismi 99 



4.5.3 Coupling with lattice phononsl 101 



4.6 



4.5.4 Nuclear spin-lattice relaxation rat3 102 

Conclusionj 105 



5 Long-range dipolar ordering in Mns 



107 

5.1 Long-range ferromagnetic orderind 108 

5.1.1 ac-susceptibilitvl 108 



Specific heati 110 



5.2 



5.1.3 Monte-Carlo Simulations! Ill 

Sph^a^c^relaxationl . 114 



5.2.1 Field-dependent specific heati 115 

5.2.2 ^^y\r^ NMRI 117 



5.3 



Conchrsionsl . . 121 



Bibliography 



123 



CONTENTS III 
SamenvattineJ 129 



List of Publications! 133 



Curriculum Vit 3 135 



Acknowledgements! 136 



I 



1 



Introduction 



Physicists, chemists and biologists have worked since decades on the study of 
systems of increasing complexity, just as much as engineers and technologists 
have been trying to miniaturize the devices they design for practical applica- 
tions. It seems that these two directions have now come very close to a joining 
point that happens to be located at the nanometer-scale. Accordingly, the name 
"Nanoscience" has been given to the resulting broad field of research. The rate 
at which laboratories and research groups are being converted to activities in 
nanoscience, may raise some fears that it will become a fashion-driven hype. On 
the other hand, it is undeniable that nanoscience is in fact nothing else than the 
natural evolution of last century's science, with a solid history of interests behind 
it. The subject of this thesis, that could be summarized as the quantum dynamics 
of magnetic molecules, is a good example of multidisciplinary topic that combines 
fundamental and practical aspects of the research at the nanometer scale. 

The field of molecular magnetism has its roots in the interest of chemists to- 
wards the synthesis of large molecules, and the assembly of macroscopic amounts 
of them in regular structures. Along this route, it became clear that the synthesis 
of molecules having magnetic ions as constituents could give rise to structures 
where each molecule can be seen as a single-domain magnetic particle often 
called "cluster", coordinated by a shell of organic ligands. Importantly, such 
molecules are stoichiometric chemical compounds that can be packed in a crys- 
talline structure, where the identical magnetic units interact only weakly with 
each other. In this way, each molecule can be treated in first approximation as a 
single nanometer-scale magnet, from where the name "Single Molecule Magnet" 
(SMM), and is characterized by a large total spin S that arises from the combina- 
tion of the atomic electron spins. The possibility of combining magnetic ions and 
organic ligands in the most diverse ways allows to tune the physical properties of 
these systems to obtain a wide range of magnetic behaviors. 

One of the most essential physical properties of SMMs is their magnetic 
anisotropy, meaning that it may be energetically favorable for the magnetic mo- 
ment of each molecule to align along a certain axis. At temperatures much lower 
than the anisotropy energy^^^ , the molecular spin is effectively "frozen" in a cer- 
tain direction along the anisotropy axis, giving rise to single-molecule magnetic 

'^'Throughout this thesis we shall always express the energies in temperature units, i.e. divided 
by the Boltzmann constant fcs. 
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hysteresis 0. This discovery in a molecular cluster containing a core of 12 man- 
ganese ions (Mni2-ac) opened the perspective of using SMMs as the ultimate 
magnetic memory units 

At the same time, people interested in the fundamental aspects of quantum 
mechanics at the macroscopic scale (the "Schrodinger cat" problem ,4 ) realized 
that SMMs could be candidates for the observation of quantum phenomena at 
the macromolecular level, in particular quantum tunneling of the magnetic mo- 
ment [Hlini- This has been indeed recently achieved by observing regular steps 
in the magnetic hysteresis loops of Mni2-ac |Z1E|; which occur when spin states 
on opposite sides of the anisotropy barrier have the same energy, so that the 
cluster spin may invert its direction by resonant quantum tunneling. Again, this 
discovery is important for both practical and fundamental reasons. On the one 
hand, it makes clear that a memory unit based on systems as small as a single 
molecule would be useless if we cannot avoid the memory being self-erased by 
quantum tunneling. On the other hand, it suggests the possibility of using SMMs 
as a test ground for theories concerning the transition from quantum to classical 
physics [SIE], in particular as far as the role of the environment is concerned. 

One of the great advantages of SMMs for fundamental research on quantum 
mechanics at the large scale is precisely that the influence of the environmen- 
tal degrees of freedom on the giant molecular spin can be accurately calculated, 
thanks to the knowledge of the crystal structure and the magnetic couplings be- 
tween the localized moments. Furthermore, the experimental investigations can 
profit of all the best known techniques for solid-state physics. Since the pure 
quantum behavior in SMMs is typically achieved only in the subkelvin tempera- 
ture range, this research requires the use of outstanding low-temperature facilities 
for magnetic measurements. In chapter II we describe the working principles and 
the guidelines for the design of several ultra-low temperature setups that are 
particularly suitable for experiments on quantum magnetism. 

As mentioned above, the magnetic properties of SMMs are determined in 
the first place at the chemistry level, but an essential aspect is that it is also 
possible, given a certain SMM, to tune its quantum mechanical behavior "in 
situ" by applying an external magnetic field. In particular, one can increase the 
tunneling rate of S by applying a field perpendicular to the anisotropy axis of 
the molecule This can be pushed to such an extent that we may expect 

the possibility to observe coherent quantum oscillations of the magnetic moment 
trough the anisotropy barrier ^21) a phenomenon that has no classical analog. 
In this case, SMMs would become suitable qubits for quantum computing I13j . 
with the interesting feature that the operating frequency can be tuned locally by 
the simple application of a magnetic field. 

A large part of this thesis is dedicated to the study of the interaction between 
SMMs and their environment, in particular the nuclear spins. This interaction is 
crucial from many points of view. In the zero- or low-field regime, the tunneling 
probability is very small, in fact so small that, until recently, there were serious 
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doubts about the observability of quantum tunneling of the magnetic moment 
in SMMs. This is due the fact that, in order to tunnel, the total electron spin 
states on opposite sides of the anisotropy barrier must be in resonance within an 
energy A that can be as small as ~ 10~^^ K in Mni2-ac. Any static perturba- 
tion larger than that, would destroy the resonance condition and make tunneling 
impossible. Prokof'ev and Stamp |14|ll5j succeeded in explaining why tunneling 
was actually observed in the experiments, by noticing that the coupling of the 
cluster spin with the nuclear spins in its surroundings, although several order of 
magnitude larger than A, is dynamic. The electron spin energy levels are there- 
fore fluctuating in time with respect to each other, which gives rise to an effective 
"tunneling window" that depends on the strength of the coupling with the nuclei, 
as was subsequently demonstrated by experiments with isotopically substituted 
samples [161. opposite regime, where the fast quantum oscillations of S 

can be obtained by applying a strong perpendicular field, the nuclear spins play 
the role of an intrinsic source of decoherence. Understanding this phenomenon is 
obviously essential for any attempt to use a magnetic qubit for quantum comput- 
ing, since nuclear spins are an unavoidable presence in practically any material. 
The decoherence due to nuclear spins is indeed important also for systems like 
the flux qubit Chapter III of this thesis is dedicated to a review of the 

theoretical aspects of the quantum behavior of SMMs, from the basic notions to 
the details of the Prokof'ev-Stamp theory. 

Despite its relevance, the nuclear spin dynamics in SMMs at ultra-low tem- 
perature remained an essentially unexplored field. In chapter IV we report a 
thorough investigation of the dynamics of ^^Mn nuclei in Mni2-ac. Our results 
uncover many basic aspects, some of which confirm the validity of the Prokof'ev- 
Stamp theory, while others are totally new and have never been considered in 
any theoretical treatment of the coupled system of "quantum spin -|- nuclei" . For 
instance, we demonstrate that the quantum tunneling of the electron spin is a 
mechanism capable of producing nuclear relaxation at an unexpectedly high rate. 
We also show that nuclei in different molecules are coupled with each other, an 
essential ingredient for the creation of a dynamic tunneling window. The anal- 
ysis of the nuclear spin dynamics in large perpendicular fields confirms earlier 
results obtained by specific heat experiments concerning the increase of the 
tunneling rate. Finally, we consider for the first time the issue of the nuclear 
spin temperature in the presence of temperature-independent quantum tunneling 
fluctuations. Surprisingly, the nuclear spin temperature is found to remain in 
equilibrium with the lattice temperature down to T ~ 20 mK! We discuss there- 
fore the need for an extension of the Prokof'ev-Stamp theory to take into account 
inelastic tunneling events. 

In chapter V we present our research on the low-temperature magnetic prop- 
erties of a rather peculiar SMM, containing Mng molecular clusters that are char- 
acterized by a very small magnetic anisotropy because of the highly symmetric 
structure. We flnd that the electron spin-lattice relaxation remains fast down to 
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millikelvin temperatures, which ahows the observation of a long-range magnetic 
ordering of the cluster spins, due only to the mutual dipolar interactions. The 
nature of this ordered state is already interesting in itself, but even more so when 
compared to the situation in the highly anisotropic SMMs, where the freezing of 
the cluster spins has made the investigation of the thermodynamic ground state 
impossible so far. Experiments in disordered rare-earth spin systems have shown 
that the application of a transverse field greatly facilitates the system in reaching 
its ground state JH|) a strategy that seems promising for anisotropic SMMs as 
well. Because of the absence of anisotropy, and therefore of quantum tunneling, 
the nuclear spin dynamics in Mng offers also an interesting comparison with the 
case of Mni2-ac. 

As may have become clear already, this thesis is focused uniquely on the 
fundamental aspects of the quantum spin dynamics of SMMs. The attempts to 
find useful applications for molecular magnets have not stopped meanwhile, and 
we think that the recent progress in patterning ^2] and surface deposition [20] of 
Mni2-ac molecules may represent an actual breakthrough for both the magnetic 
storage and the quantum computing purposes. Wherever this progress will lead 
to, we hope that the work presented here will help the scientific community to 
gain a deeper understanding of how the quantum phenomena at the molecular 
scale are influenced by the interaction with the environment. 
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This chapter begins with a survey of the principles and the design of the setups 
used for the experiments reported in chapters IV and V. In addition, we discuss 
the design and construction of two ultra low-T setups for static magnetic mea- 
surements, a SQUID magnetometer and a torquemeter, which are not involved in 
the research presented in this thesis, but constitute a new and significant exten- 
sion to the experimental facilities of our laboratory: ^2.41 and ^2.51 are therefore 
meant as future reference for the newcomers in our research group, or anybody 
interested in building ultra low-T magnetometers. 

2.1 Dilution refrigerator 

The heart of of our ultra low-T setup is a ^He/^He dilution refrigerator. The 
working principle (see e.g. i2lj for details) is based on the fact that, even at 
T = 0, a liquid mixture of ^He and ^He does not separate completely: the Van 
der Waals force between a ^He and a "^He atom is larger than between two ^He, 
thus a certain fraction of ^He atoms will spontaneously dilute into the ^He phase, 
until the Fermi energy of ^He in the dilute phase equals the difference in binding 
energies between ^He-^He and ^He-^He. If one distils some ^He out of the dilute 
phase, more ^He atoms will migrate from the pure ^He phase into the dilute one to 
reestablish the equilibrium. The increase of entropy associated with this process 
corresponds to the absorption of heat by the ^He/^He mixture, i.e. an effective 
cooling power. The cooling can be made continuous if the ^He extracted from 
the dilute phase is recondensed and injected on the side of the pure ^He phase, 
taking care that the incoming ^He stream is properly cooled by exchanging heat 
with the cold ^He atoms that flow through the dilute phase. 

Our system is based on a Leiden Cryogenics MNK126-400ROF dilution re- 
frigerator, fitted with a plastic mixing chamber that allows the sample to be 
thermalized directly by the '^He flow, specially designed for our purposes in col- 
laboration with G. Frossati. A scheme of the low-temperature part of the refrig- 
erator is shown in Fig. 12. H together with the NMR circuitry (see H2.2.2|) . The 
mixing chamber consists of two concentric tubes, obtained by rolling a Kapton 
foil coated with Stycast 1266 epoxy. The tops of each tube are glued into con- 
centric Araldite pots: the inner pot receives the downwards flow of condensed 
^He and, a few millimeters below the inlet, the phase separation between pure 
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(a) Coaxial cable connected to 
the NMR spectrometer and 
pulse generator. 

(b) 'He distillator ("still") with cold plate. 

(c) Rotating shafts for the tunable 
capacitors, accessible from the top 
of the refrigerator. 

(d) NMR matching capacitor. 

(e) NMR tuning capacitor 

(f) Upper heat exchanger. 

(g) 80 mK pot. 

(h) X-cable, cut to one wavelength 
at 280 MHz. 

(i) Lower heat exchanger 
(j) Araldite mixing chamber 
(k) Upper thermometer. 

(I) Pure 'He phase. 

(m) Double-wall Kapton tail. 

(n) Forced downwards flow of 'He 

in the dilute phase, 
(o) Sample with NMR coil, 
(p) Openings in the inner Kapton tube 

to allow the return of the 'He flow, 
(q) Lower thermometer 
(r) Conical vacuum plugs. 



Fig. 2.1 : Scheme of the low-T part of the dilution refrigerator with NMR 
circuitry. For clarity we omit the vacuum can and the radiation shields anchored 
at the still and at the 80 mK pot. Only the narrow tail of the refrigerator is 
inserted in the bore of a 9 T superconducting magnet (not shown). 



2.1. Dilution refrigerator 
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^He phase and dilute '^He/'^He phase takes place. The circulation of '^He is then 
forced downwards along the inner Kapton tube, which has openings on the bot- 
tom side to allow the return of the ^He stream through the thin space in between 
the tubes. Both the bottom of the Kapton tail and the outer pot are closed 
by conical Araldite plugs smeared with Apiezon N grease. The Kapton tail is 
about 35 cm long and is surrounded by two silver-coated brass radiation shields, 
one anchored at the 80 mK pot, the other at the still. The whole low-T part of 
the refrigerator is closed by a vacuum can, which has itself a thin brass tail to 
surround the Kapton part of the mixing chamber, and is inserted into the bore 
of an Oxford Instruments 9 Tesla NbTi superconducting magnet. In this way, 
only Kapton and brass cylinders (plus the lowest Araldite plug) are placed in the 
high-field region, whereas all other high-conductivity metal parts (heat exchang- 
ers, 80 mK pot, ^He distillator, etc.) are outside the magnet bore and subject 
only to a small stray field. This design was intended to minimize eddy currents 
heating while moving the refrigerator through the pick-up coil of the SQUID 
magnetometer f ^2.4() . while thermalizing the sample directly by the contact with 
the ^He flow. As it turned out, the excellent thermalization obtained is this way 
is also essential for the success of the NMR experiments described in chapter IV, 
whose significance depends crucially on the efficient cooling of the nuclear spins 

(EM- 

The temperature inside the mixing chamber is monitored by measuring with 
a Picowatt AVS-47 bridge the resistance of two Speer carbon thermometers, one 
in the outer top Araldite pot, the other at the bottom of the Kapton tail, just 
beside the sample. We have verified that the temperature is very uniform along 
the whole chamber (except in the presence of sudden heat pulses): even at the 
lowest T, the mismatch between the measured values is typically < 0.5 mK. A 
Leiden Cryogenics Triple Current Source is used to apply heating currents to a 
manganin wire, anti-inductively wound around a copper joint just above the ^He 
inlet in the mixing chamber. In this way we can heat the incoming ^He stream 
and uniformly increase the mixing chamber temperature. 

For the ^He circulation we employ an oil-free pumping system, consisting of a 
Roots booster pump (Edwards EH500) with a pumping speed of 500 m'^/h, backed 
by two 10 m^/h dry scroll pumps (Edwards XDSIO). The main pumping line is a 
100 mm solid tube, fixed at one side with a fiexible rubber joint that allows for 
an inclination of a few degrees, and connected to the head of the fridge by a "T" 
piece with two extra rubber bellows, to reduce the vibrations transmitted by the 
pumping system. In this configuration, the system reaches a base temperature 
that can be as low as 9 mK; with the extra wiring for NMR experiments, the base 
temperature is ~ 12 mK, and the practical operating temperature while applying 
r/-pulses is 15 - 20 mK (see H2.2.3() . The typical ^He circulation rate at the base 
temperature is n ~ 350 /xmol/s, and the cooling power at 100 mK is Q ~ 150 
/iW. When the Kapton tail is replaced by a flat plug, Q can be increased up to 
700 iiW @ 100 mK and n ~ 1200 /nmol/s by applying extra heat to the still; with 
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the tail in place it's more difficult to increase the circulation rate, and Q @ 100 
mK hardly exceeds 250 fiW. 

2.2 Nuclear Magnetic Resonance 

The NMR technique is applied since more than 50 years H2] to widely different 
research fields; in this thesis, we are interested in the use of nuclear spins as local 
probes for magnetic fluctuations in solid-state systems, but also in the dynamics 
of the nuclei themselves, which appears to play an essential role in the quantum 
behavior of single-molecule magnets. 

2.2.1 Basics of pulse NMR 

The Hamiltonian of a nuclear spin / placed in a magnetic field applied along 
a certain axis z is: 

n = hjr-B,, (2.1) 

where 7 is the gyromagnetic ratio, such that fi = h^yl is the magnetic moment 
of the nucleus. The eigenstates \m) are the 2/ -|- 1 projections of the spin along 
i*, having energies m = —I . . . I, separated from each other by the Zeeman 
splitting Em+i — Em = fi'jBz, which corresponds to the Larmor frequencies lo = 
'-fBz for the classical precession. Transitions between adjacent Zeeman levels can 
be produced by introducing a time-dependent field Brf{t), provided that the 
matrix elements {m\Brf\m -|- 1) 7^ 0, i.e. B^f is not parallel to z. Considering 
for simplicity a spin / = 1/2 and taking Brf{t) = B^fX cos{uJi\it), the expectation 
value of the z component of the magnetic moment will oscillate in time at a 
frequency WRabi = T-^r/ (Rabi oscillations): 

(/x,(t)) = (/i,(0))cos(wRabit). (2.2) 

In a classical picture, this means that, if /2(0) = h'yiz, after a time = 
7r/{2'yBrf) the magnetic moment has been turned 90° away from the z axis and 
lies in the xy plane. For this reason, the application of such an alternating field 
of frequency ujjy for a time is called "7r/2-pulse" . The rotation angle can in 
fact take any value, by choosing the appropriate duration and strength of Brf- 
After a 7r/2-pulse the state of the system is {ip) = l/\/2(| -1-1/2) -|-| — 1/2)), and 
the time evolution caused by the Hamiltonian (|2.1j) is equivalent to the classical 
Larmor precession of the magnetic moment within the xy plane, which produces 
a rotating magnetic field that can be detected via the electromotive force induced 
in a pick-up coil with axis || x. In practice the same coil is used to produce Brf 
and to detect the Larmor precession. 

Since we always work with macroscopic ensembles of spins, it is convenient to 
describe the system in terms of its density matrix p. If the state of each spin 
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in the system is described as = X^m^mlw-)) then the matrix elements Pm,m' 
are the average over the sample of amO*^i ■ The diagonal elements of p represent 
the populations of the levels, and the non-diagonal elements account for the 
correlation between different spins. Under thermal equilibrium at a temperature 
T, the non-diagonal terms of p are zero, and the diagonal terms obey: 



which corresponds to a longitudinal (|| B^) nuclear magnetization, Mz{T), ac- 
cording to the Curie law: 



where is the number of spins per unit volume. 

By means of r/-pulses as described above, it is possible to perturb the system 
and change its density matrix. For instance, the effect of a vr-pulse on an ensemble 
of spins 1/2 in thermal equilibrium is to invert the diagonal elements of p, i.e. 
the nuclear magnetization. After the inversion, the thermal equilibrium can be 
reestablished in the presence of perturbations that induce transitions between 
the Zeeman levels, and of a reservoir that can absorb the heat released by the 
nuclear spins. The time constant for the recovery of the equilibrium values of the 
diagonal elements of p is Ti, the nuclear spin-lattice relaxation time. 

Far from equilibrium, it is possible to define a nuclear spin temperature T^nd 
different from the lattice T, provided that the ratio of the populations formally 
obeys a relation like 1)2. 3(1 with T^^ci replacing T, and that Pm,m' = for m ^ m'. 
For example, a vr-pulse on a spin 1/2 produces a negative spin temperature. Ti 
can be then interpreted as timescale for the reequilibration of the spin- and lattice 
temperatures. 

After a 7r/2 pulse, the diagonal elements of p are identical and nonzero non- 
diagonal elements appear. This means that the spins are coherently precessing, 
the correlation between the single Larmor precessions decaying with a time con- 
stant T2 < Ti, which is also the decay time for the non-diagonal elements of 
P- 

In a real experiment, the signal produced by the Larmor precession after a 
7r/2-pulse, called "free induction decay", actually decays on a time T2 < T2 
because of static field inhomogeneity, i.e. the fact that spins in different spatial 
positions may have different Bz, thus different ujn- This static dephasing can 
be recovered by means of spin echo. By applying a vr-pulse at a time r > Tg* 
after the 7r/2-pulse, each spin is rotated 180° about the axis of the coil. In 
this way, the spins that are precessing "too slow" find themselves ahead of the 
fast ones. At time 2r all spins are again precessing in phase (Fig. 12. 2() . giving 
rise to an echo of coherent precession whose amplitude depends on the effective 
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(a) (b) (c) (d) 




Fig. 2.2 : Sketch of the spin-echo sequence: (a) tt /2-pulse; (h) Larmor preces- 
sion and dephasing of the slow (light gray) and fast (dark gray) spins; (c) n-pulse; 
(d) refocusing of the precessions and formation of the echo. 
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Fig. 2.3 : Left: T2 measurement by echo decay. Right: Ti measurement by 
inversion recovery. 



correlation oc exp(— 2T/T2) at that instant. We can therefore measure T2 by 
applying a it/2 — tt pulse sequence with increasing values of r, and observing the 
r— dependence of the echo amplitude. 

As regards Ti, we remark that the coil is sensitive only to the projection of 
the Larmor precession along its axis, i.e. _L z. The effect of a 7r/2-pulse, or of an 
echo sequence, is indeed to project the z-component of the nuclear magnetization 
into the xy plane and make it measurable. For instance, an echo sequence after 
a TT-pulse yields a signal which is proportional to the initial magnetization, but 
with opposite phase as compared to the signal without vr-pulse. By introducing 
a waiting time r between the vr-pulse and the echo sequence, one can access the 
time evolution of the z-component of the magnetization, i.e. the recovery of 
the equilibrium state and thereby the time constant Ti. A sketch of the pulse 
sequences used to measure Ti and T2 is shown in Fig. 12.31 
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2.2.2 Instruments and circuitry 

A scheme of the NMR setup used for the present work is shown in Fig. 12.41 
The non-commercial electronic devices have been developed by R. Hulstman, 
and a computer running a Pascal management program written by J. Witteveen 
[2,Sj controls the generation of the NMR pulses and the detection of the signal. 
The pulse sequences are obtained by programming a home-made device, which 
opens the gate of the pulse amplifier (Kalmus 166UP) at the specified times: 
one typically waits 5 before the full power of the amplifier (typ. ~ 100 W) is 
available, and closes the gate as soon as the power is no longer necessary, to avoid 
leakage currents to the NMR coil. The r/-signal to be amplified is generated by 
the frequency synthesizer (Farnell SSGIOOO), but is first handled by a r/-gate, 
driven by the pulse programmer. The purpose of the r/-gate is to produce pulses 
of the desired length and phase. For instance, to measure a spin echo we first 
produce a 7r/2 — vr sequence with all the pulses in phase with the reference signal, 
and record the result; then we repeat the sequence but now with a 180° phase 
shift on the 7r/2-pulse, and we subtract the measured signal from the previous 
one. In this way, the "true" NMR signal from the second sequence, which has 
opposite sign with respect to the first because of the phase shift in the 7r/2- 
pulse, is actually added to the first one, whereas we cancel out the offset of the 
spectrometer and the ringing (i.e. the fictitious signal that appears because of the 
electrical resonance of the circuit) produced by the 7r/2-pulse j24j. The amplified 
pulses are then fed to the resonant circuit. 

The circuitry we used includes two tunable cylindrical teflon capacitors, mounted 
at the still of the dilution refrigerator (see Fig. 12. 1|) and driven by rotatable 
shafts. The capacitor in parallel to the r/-line is used to match the impedance 
of the circuit to 50 0,, whereas the one in series to the NMR coil tunes the fre- 
quency of the resonator. The peculiarity of our circuit is that, because of the 
distance (~ 50 cm) separating the capacitors from the NMR coil in the mixing 
chamber, we had to introduce a A-cable between coil and tuning capacitor. At 
the frequency where the cable is precisely one wavelength, the circuit is indeed 
identical to the standard concentrated circuit. Away from the A- frequency, the 
cable adds extra inductance or capacitance. More importantly, since the A-cable 
is a low-conductivity thin brass coax for low-T applications, the quality factor 
of the resonator (which includes the cable!) is drastically reduced. Although 
this affects the sensitivity of the circuit, it also broadens the accessible frequency 
range without need to retune the capacitors. We found that, cutting the cable 
for one wavelength at ~ 280 MHz, the circuit is usable between (at least) 220 
and 320 MHz. Once the desired frequency range is chosen and the A-cable is cut 
accordingly, the NMR coil must be constructed by first mounting it just near the 
tuning capacitor, i.e. making a "standard" concentrated circuit. If this circuit 
resonates in the same range as the A-frequency, the coil can be safely moved into 
the mixing chamber. 
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Fig. 2.4 •■ Scheme of the NMR setup. The parts enclosed in the dashed box are 
inserted in the dilution refrigerator. 
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Time (min) 

Fig. 2.5 : Time evolution of the mixing chamber temperature while applying 
12 fis - 24 l^s spin-echo pulse sequences, (a) Response of the thermometer at 
the bottom of the Kapton tail, next to the sample, after two spin-echo sequences 
at 10 min intervals, (b) Like (a), after 20 sequences at 1 min intervals, (c) 
Thermometer at the top of the mixing chamber, close to the ^He outlet, after a 
single spin-echo sequence. 
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The NMR signal induced in the coil is first passed though a low-noise DOTY 
preamplifier placed at the top of the cryostat, then detected by a phase-sensitive 
quadrature detector, and finally stored in the memory buffer of a flash A/D con- 
verter. The results of the positive and negative pulse sequences are subtracted 
and dumped in the measurement data file. The diodes in parallel to the pream- 
plifier and the A/4-cable, together with the diodes in series to the power line, 
are used to separate the small NMR signal from the high-power pulses. For the 
pulse amplifier, the diodes in series are a short circuit, whereas the A/4-cable is 
an open circuit since it is short-circuited at the opposite end by the diodes in 
parallel. In this way all the power goes in the resonator and does not damage the 
preamplifier. Conversely, the series diodes are an open circuit for the very small 
signal picked up by the NMR coil, which is therefore entirely conveyed towards 
the preamplifier. 

2.2.3 Heating effects 

Applying NMR pulses of typically 100 W for several microseconds to a copper 
coil at T ~ 15 mK is potentially very harmful for the thermal stability of the 
sample and of the ^He/^He mixture in which it is immersed. As discussed in ^'2.11 
our setup is specially designed to circumvent this problem. In Fig. 12.51 we show 
three examples of the magnitude of the heating effects that arise when applying 
spin-echo NMR pulse sequences consisting of a 12 ^s 7r/2-pulse followed by a 24 
IIS TT-pulse after 45 fis. Panels (a) and (b) show the temperatures as measured by 
the thermometer placed at the bottom of the Kapton tail, just next to the NMR 
coil and the sample, whereas the temperature in panel (c) is taken at the top of 
the mixing chamber in the outer part of the wide Araldite pot, i.e. 35 cm past 
the sample in the ^He path. In both cases, a very sharp peak in T is visible at 
the instant when the pulses are applied, which is simply due to the direct (and 
instantaneous) heating of the thermometers by the electromagnetic field produced 
by the coil, thus it does not reflect a "real" increase in the temperature of the ^He 
bath. The situation is particularly clear in panel (c), where the radiative heating 
spike is quickly recovered, and only after a few minutes the wave of warmer ^He 
atoms comes along. By comparing (a) and (b) with (c), it is also clear that the 
height of the radiative spike is lower in the upper thermometer, which is indeed 
farther away from the coil. Panel (a) demonstrates that, with an interval of 10 
minutes between the sequences, the temperature of the sample can be kept very 
stable at ~ 15 mK. Even with just 1 minute waiting time (b), the increase of the 
baseline of the ^He bath temperature does not exceed 5 mK! 

2.3 ac - susceptometer 

For the frequency-dependent ac-susceptibility measurements reported in chapter 
V we used a simple and compact home-made susceptometer based on the mutual 
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Fig. 2. 6 : Schematic drawing of the ac-susceptometer. 



inductance principle. A sketch of the instrument is shown in Fig. 12.61 A primary 
coil, consisting of 250 turns of 100 /.im NbTi superconducting wire is wound 
on top of two secondary coils, each consisting of 600 turns of 40 ^m copper 
wire, and connected in series with opposite polarity. When an alternating current 
Ip{t) = IpCos{Lot) is passed through the primary coil, it induces a voltage Vg = 
^1 — V2 in the secondary, where Vi and V2 are the voltages at each section of the 
secondary coil, and the minus sign is due to the opposite winding directions: 

Vs = ujIp{Mi-M2)cos{ujt + 7r/2). (2.5) 

The above relation holds if there is no current flowing in the secondary coils, e.g. 
when they are connected to an ideal voltmeter. The mutual inductances are given 
by M-i = {1 + Xi)f-, where (i = 1, 2) and/ = f{Np, Ng, Ap, As,. . .) is a factor that 
depends on the geometry, the area and the number of turns of the primary and 
secondary coils, and Xi is the magnetic susceptibility of the material filling each 
secondary coil. By leaving coil 2 empty {x2 = 0) and filling coil 1 with the sample 
to be investigated, we obtain: 

Vs = ujIpCos{LJt + -K/2)x- (2.6) 

The susceptibility of the sample is therefore easily measured by connecting the 
secondary coil to a lock-in amplifier (Stanford SR830); the phase-sensitive de- 
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tection allows to discern both the real and the imaginary part of the complex 
susceptibility X = x' + ^x" • 

To apply the primary current Ip we either introduced a 500 kfi resistor in 
series to the 5 V excitation circuit of the lock-in, or we used a home-built V/I 
converter that allows to have more current at low frequencies, since the induced 
voltage is oc u;. Using the V/I converter we may easily measure down to frequen- 
cies ~ 1 Hz, whereas above ~ 1 kHz it is advisable to simply use the 500 kSl 
series resistor, which eliminates the slight phase rotations that the V/I converter 
starts to introduce at high frequencies. With a typical value of Ip ~ 50 /xA, the 
alternating field produced at the sample is noHp ~ 1 ^T. 

The susceptometer is small enough to fit vertically in the tail of the dilution 
refrigerator, or even horizontally in the outer part of the upper Araldite pot (cf. 
Fig. 12. 1() . for measurements without static external field. 

2.4 SQUID magnetometer 

The SQUID magnetometer described in this section is designed to allow high- 
sensitivity measurements on powder samples down to T ~ 10 mK. For example, 
metallic nanoclusters show spectacular quantum-size effects in their thermody- 
namic properties [2^] at very low temperatures, but their magnetic susceptibility 
is very small in this regime. Furthermore, most of those materials are highly air- 
sensitive, thus it is very convenient to introduce and keep the sample in a sealed 
glass tube, which is also ideal to perform preliminary measurements at T > 1.5 
K in commercial SQUID magnetometers. One of the requirements for the design 
discussed here is indeed the compatibility with the sealed sample holders used 
in other experiments, without any need to further manipulate the sample. The 
really challenging part of the design consists in moving the sample through the 
pick-up coils while it is at T 10 mK inside the mixing chamber of the dilu- 
tion refrigerator. For samples with such extremely small magnetic signals such 
a measure is needed because the two components of the pick-up coil will never 
perfectly compensate one another, leaving an empty-coil signal that could wash 
out the signal of the sample to be measured. As we shall discuss below, this has 
required a system that moves the whole dilution refrigerator insert. 

2.4.1 Working principle 

A SQUID (Superconduting QUantum Interference Device) is basically an ultra- 
sensitive flux-to-voltage converter, that exploits the peculiar quantum properties 
of closed superconduting circuits. To understand its working principle, we recall 
that the superconducting state can be described as the condensation of paired 
electrons (Cooper pairs) into an ordered phase characterized by a complex order 
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parameter: 

^P = |V;(r)|e^^(*^, (2.7) 

where is the density of Cooper pairs and (p{r) = p • r/h is the phase. 

In the presence of a magnetic field B = V x A, the generahzed momentum is 
expressed as p = 2mv — 2eA, since the Cooper pairs have mass 2m and charge 
2e (e = -1.6 x IQ-^^ C). 

The current density j can be obtained by introducing (|2.7|1 in the standard 
quantum mechanical expression: 

j = -ih— fv*VV' - i^Vtp*) - (2.8) 
m \ / m 

J=(Mv^-f!!^-)|rf. (2.9) 

\ m m J 

The current in a superconductor can flow only within a surface layer of thickness 
comparable to the London penetration depth, since inside the bulk j = 0. By 
considering a superconducting ring, the integral § j • dl along a closed path deep 
inside the bulk is thus obviously zero, i.e.: 

( Mv(^ - —a] ■ dl = 0; (2.10) 
\ m m J 

recalling the Stokes theorem, ^^•d/ = Jj V x A ■ dS = JJ B ■ dS = <I> yields 
the magnetic flux enclosed by the ring. Furthermore, since the order parameter 
must be single- valued, the total phase accumulated along a closed path must be 
an integer multiple of 2tt: 

V(/7 • dl = 2Tm. (2.11) 

Combining (|2.1U|) and 1)2. llj) yields the quantization of the magnetic flux in a 
superconducting ring: 

^ = = n$o, (2.12) 

where $0 = h/2e = 2.07 x 10~^^ Wb is the flux quantum. This means also that, 
since the externally applied flux <I>ext is a classical variable and can be changed 
smoothly, the quantization of the flux may require a screening current Is to flow 
in the ring to fulfill the condition: 

$ = $ext + Lis = ra$o, (2.13) 



where L is the self-inductance of the ring. 
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Fig. 2. 7 : Sketch of superconducting ring interrupted by a weak link. The dashed 
line is the integration path. 

A SQUID is obtained by interrupting a superconducting ring with one (r/- 
SQUID) or two (dc-SQUID) Josephson junctions, i.e. weak links where the su- 
perconductivity is locally suppressed but the Cooper pairs can cross by quantum 
tunneling. The essential feature added by a weak link is that it is no longer 
possible to find a closed path through the ring such that j = everywhere, since 
j ^ through the weak link (Fig l2.7|) . This adds an extra contribution if' to the 
phase of the order parameter: 



^' = ^ni2 ^-^^^ (2-14) 

eh\W Jw.i. 



where the integral is along the weak link, supposed to be so thin that the magnetic 
flux in negligible, and \ip\'^ in the limit of small currents can be taken equal to 
the bulk value. The condition of single-valued order parameter becomes: 

+ ^ [[ B-dS = ^' + 27r— = ^' + 27r ^""^ ^ = 2^n. (2.15) 

h JJ ^>o ^>o 

It follows that the screening current /«(<&) must be a periodic function of ^ext/'^'o- 
Is = Ic sin(27r^>/$o) = /csin[27r(n - $/^>o)] = /csinv?', (2.16) 

which is the famous Josephson phase-current relation. Ic is the Josephson critical 
current, i.e. the maximum supercurrent that can pass through the link without 
dissipation. 

By applying an external flux which varies linearly in time, a constant voltage 
V = — d^>ext/df develops across the weak link. Substituting in Eq. (|2.16j> yields: 

/, = /,sin (2^ ^^^^^^ )- (2.17) 
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Fig. 2.8 : Model of current-biased Josephson junction (J) with dissipation (R) 
and capacitance (C), and the relative V — I characteristic. 

This implies that, when a constant voltage is applied to the Josephson junction, 
an alternating supercurrent circulates in the ring with frequency v = 2eV/h 
proportional to the applied voltage. Combining (|2.17jl and (|2.16|) we find that 
a constant voltage produces a linear increase of (p' with time, i.e. the voltage is 
proportional to the time derivative of the phase: 

y = A^. (2.18) 

2edt ^ ^ 

This relation is useful to obtain the V — I characteristic of a current-biased 
Josephson junction. When I > 1^ the current through the weak link causes 
dissipation, which can be accounted for by a resistance R in parallel to the ideal 



Josephson junction described by (|2.16|) . Moreover, a realistic junction is formed 
by two superconducting pads separated by a very thin barrier, so that we have 
to add a capacitance C to the model. We may then write: 

^ = ^^^^^^ + + ^ d^ = ^^^"^^ + 2^1 + ^'-''^ 

which leads to aV/I curve of the type shown in Fig. 12.81 

A (ic-SQUID, like the one used in our setup, is operated by biasing the junc- 
tions Ji and J2 with a current /bias > Ic, such that a voltage develops across 
them (Fig. 12.9(1 . The essential feature is that the voltage can be modulated by 
an applied flux, since the critical current Ic also depends on Indeed, after 
some calculation one finds that the screening current Is is related to the critical 
current Ic of the junctions (assumed to be identical) by: 



± o 

2 



/ -tin ^cxt + Lis 

sm — sm I ip^^ — Zir- 



(2.20) 



Ic is obtained by maximizing Ig, yielding the periodic form of /c(*I'cxt) shown 
in Fig. 12.91 By properly choosing the bias conditions, the (ic-SQUID operates 
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Fig. 2.9 : Working principle of a dc-SQUID as flux-to-voltage converter. 

therefore as a flux-to- voltage converter, where the external flux can be applied by 
injecting a current in the input coil. Notice that small a fraction of a flux quantum 
can produce voltage changes of the order of millivolts! In the practice a dc- 
SQUID is used in feedback mode by employing a so-called Flux-Lock Loop (FLL), 
i.e. adding a feedback coil that produces a compensating flux such that V = 
const. This increases the accuracy and the dynamic range of the measurement, 
and allows to implement noise-reducing detection schemes. For more details on 
SQUID sensors, see 

To use a dc-SQUID in an actual magnetometer, it is still necessary^^^ to pro- 
duce a current proportional to the magnetic moment of the sample to be mea- 
sured. Such a current can be injected in the input coil to produce a flux that is 
coupled to the SQUID ring by the mutual inductance Ai. Typically, the input 
current is obtained by constructing a closed superconducting circuit which in- 

Except in a rather radical design like the "microSQUID" |28| 
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eludes the SQUID's input coil on one side, and terminates with a pick-up coil on 
the other side. Once the circuit has been cooled down below the superconducting 
critical temperature, the enclosed flux is constant. Any change in the magnetic 
permeability of the circuit, possibly due to the sample, will result in a screening 
current that, while keeping the total flux constant, produces the required flux in 
the input coil. In our case, the change in permeability of the pick-up circuit is 
obtained by vertically moving the whole dilution refrigerator ( ^2.4.3() . whose tail, 
that contains the sample, is inserted in the pick-up coil. 

Because of the high sensitivity, it is essential to make sure that no sources of 
flux other than the sample may couple with the SQUID. This can be done by 
employing a gradiometer, which in its simplest form consists of two coils wound 
in opposite direction, so that the flux produced by any uniform magnetic field 
cancels out, and only the gradient of B can be detected (first-order gradiometer). 
A further improvement is the second-order gradiometer, shown in Fig. I2.1U1 
which is obtained by inserting a coil with 2A^ windings between two coils with 
N turns each, wound opposite to the central one. This design eliminates also 
the effects of linear field gradients. Furthermore, according to the reciprocity 
principle [^Hl j the flux $ produced in a coil of arbitrary geometry by a magnetic 
moment n at position r is related to the field B{f) produced by the same coil in 
r when carrying a current / such that: 



The field produced by a single-loop coil is equivalent to the field of a magnetic 
dipole, whereas a first-order gradiometer is a magnetic quadrupole, and a second- 
order is an octupole, so that the fields they produce vary in space like r~^, 
and r~^, respectively. From Ea. ()2.21|l it is clear that a second-order gradiometer 
is the least sensitive to magnetic fields produced outside of the coil system. 

The magnetic moment of the sample can be detected by moving it through 
the pick-up coil. The enclosed fiux is easily obtained from the flux induced by a 
dipole with magnetic moment /2 || z at a position z along the axis in a loop of 
radius a placed at z = zq: 



If the upper and lower coils of the gradiometer are placed at z = d and z = —d, 
respectively, then the total picked-up flux is: 



B- fl = '^I. 



(2.21) 




^>loop(2) = y/(^ - ^o)M 
r .2 / . X 21 -3/2 



(2.22) 



(2.23) 



V(^) = N^[f{z -d)- 2f{z) + f{z + d)]^,. 



(2.24) 
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Fig. 2.10 : A second-order gradiometer in Helmoltz geometry (a = 2d) and the 
magnetic flux induced by a dipole moving along z. Notice that the side peaks in 
$(z) do not coincide with the positions of the external coils. 



Typically one adopts the Helmoltz geometry, o = 2d, resulting in a ^{z) as shown 
in Fig. TTWi ^l 

The picked- up flux is related to the screening current in the circuit, Is, by: 

^pu = (-^^pu + Aeads + -^^in)^, (2.25) 

where Lpu, -^^icads and Lin are the inductances of the pick-up coil, the leads and 
the SQUID input coil, respectively. The flux at the SQUID sensor $SQUID is thus 
given by: 

$SQUID = Mis = /tr^'pu, 

/ (2 26) 

Lpu + Lleads + -^in 

where fti is the flux-transfer ratio [HO] . 
'^^The factor a in the a;-axis scale of Fig. I2.1UI is obtained by using the second form of Eq. 
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2.4.2 Design and construction of the pick-up circuit 

The circuitry for our SQUID magnetometer is, with a few additions, based on the 
principles discussed above. The niobium dc-SQUID sensor is part of a Conductus 
LTS iMAG system, which includes a FLL circuitry to be placed just outside the 
cryostat, and is connected to the SQUID controller by a hybrid optical-electric 
cable. To couple a magnetic signal to the SQUID we constructed a second- 
order gradiometer by winding a 100 /im NbTi wire on a brass coil-holder, to 
be inserted in the bore of the Oxford Instruments 9 T NbTi superconducting 
magnet. The gradiometer coils have 36 mm (a = 18 mm) and d = 9 mm, 
due to the 34 mm vacuum can of the refrigerator that contains the sample. 
The leads of the pick-up circuit, which are tightly twisted and shielded by a 
Nb capillary, are screwed onto the input pads of the SQUID to obtain a closed 
superconducting circuit. The input inductance of the SQUID sensor is Lin = 600 
nH; since -Lieads ~ 200 nH and, given the dimensions, Lpu ~ Lm + Lieads already 
with A'^ = 1, it follows from Eq. ()2.26|) that the most convenient choice of windings 
for the gradiometer^^) is 1-2-1. The mutual inductance between SQUID ring and 
input coil is = 10 nH, which means that /tr ~ 1/200. 

The FLL included in the Conductus electronics is able to compensate for 500 
$0 at the SQUID ring, which means that a maximum flux ^I'p™^^^ = 500//tr ~ 10^ 
$0 can be picked up by the gradiometer without saturating the system. The 
maximum magnetic moment 

^(max) therefore [cf. Eq. (ITM and Fig. I^TTTj : 



In order to extend the dynamic range, we have built an extra flux transformer on 
the pick-up circuit, which allows to introduce a magnetic flux from the outside 
to compensate for the flux induced by the sample. By using the SQUID as a 
null- meter, there is the extra advantage that no current circulates in the pick- 
up circuit, thus the field at the sample is precisely the field produced by the 
magnet, without the extra field that would be produced by a current in the 
gradiometer. The flux transformer consists of 16 turns of NbTi wire, wound 
on top of 1 loop of the pick-up wire and shielded by a closed lead box. In the 
same box we glued the pick-up leads on a 100 O chip resistor, which is used to 
locally heat the circuit above the superconducting and eliminate the trapped 
flux. The flux transformer is fed by a Keithley 220 current source, whereas the 
heater is operated by one of the three current sources in the Leiden Cryogenics 
Triple Current Source used for the refrigerator; both sources are controlled by a 
computer program (see ^2.4.5|) 

Despite the efforts to isolate the system from external vibrations, the powerful 
pumps of the dilution refrigerator still provide a non-negligible mechanical noise. 
In particular, the Roots pump has a vibration spectrum with a lowest peak at 

'^'Recall that $pu cx A'' but Lpu oc . 



(max) 



2a 



^^'^^^^ ~ lO^Vs ~ 10"^ emu. 



(2.27) 



0.6/io 
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~ 25 Hz. To prevent the possible flux changes induced by such vibrations when 
coupled to a magnetic field, we added an extra low-pass filter on the pick-up 
circuit |HJ in the form of a 17 mm long 0.3 mm copper wire in parallel with 
the pick-up leads. At T = 4.2 K such a wire has a resistance i? ~ 11 fiQ, which 
together with the input inductance Lin = 600 nH of the SQUID yields a cutoff 
frequency 

Ap. = T^^=^3Hz. (2.28) 

2tT Lin 

The filter is contained in a separate Pb-shielded box inserted between the flux 
transformer and the SQUID. The pick-up leads and the copper wire are contacted 
via Nb pads. 

A picture of the SQUID circuitry is shown in Fig. 12.111 The SQUID, the 
filter and the flux transformer are mounted on a radiation shield of the magnet 
hanging. The whole system is inserted in a 90 liter helium cryostat. The magnet 
hanging is stabilized by triangular phosphor-bronze springs (not shown in Fig 
I2.11|) that press against the inner walls of the cryostat. 



2.4.3 Vertical movement 

The movement of the sample through the gradiometer is obtained by lifting the 
whole dilution refrigerator. For this purpose, the refrigerator is fixed on a movable 
flange, while the top of the cryostat is closed by a rubber bellow. On the flange 
we screwed the nuts of three recirculating balls screws (SKF SN3 20 x 5R), which 
allow a very smooth displacement of the nut by turning the screw^^-*. The base 
of each screw is mounted on ball bearings and is fltted with a gearwheel. The 
gearwheels are connected by a toothed belt (Brecoflex 16 T5 / 1400) driven by 
a three-phase AC servomotor (SEW DFY71S B TH 2.5 Nm) that can exert a 
torque up to 2.5 Nm. In this way, the rotation of the servomotor is converted 
into the vertical movement of the dilution refrigerator insert. The movement is so 
smooth that the consequent vibrations are hardly perceptible and do not exceed 
the vibrations due to the ^He pumping system. 

The motor is driven by a servo-regulator (SEW MDS60A0015-503-4-00) that 
can be controlled by a computer. For safety reasons an electromagnetic brake 
is fitted, that blocks the motor in case of power failure. In addition, a set of 
switches is mounted along one of the pillars that support the screws: when the 
flange reaches the highest or the lowest allowed position, the switches force the 
motor to brake independently of the software instructions. The maximum allowed 
vertical displacement is 10 cm. 

A picture of the top of the cryostat with the vertical movement elements is 
shown in Fig. 12.121 

'■^'The friction between nut and screw is so low that, by placing the screw vertically and leaving 
the nut free, the nut would start to turn and fall off the screw just by gravity! 
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Fig. 2.11 : Picture and scheme of the SQUID circuitry, (a) Superconducting 
magnet, (h) Coil holder, (c) Flux transformer and heater in lead shield, (d) 3 
Hz lead shielded low-pass filter, (e) SQUID sensor, (f) Pick-up leads with Nh 
capillary shield, (g) SQUID cryocahle. 
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Fig. 2.12 : Picture of the top of the cryostat, with the elements for the vertical 
movement, (a) Head of the dilution refrigerator, (h) Movable flange, (c) Rubber 
bellow, (d) Screws, (e) Nut with recirculating balls, (f) Gearwheels, (g) Belt, 
(h) Servomotor, (i) SQUID FLL electronics, (j) Feedthrough to the SQUID 
cryocable. 
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2.4.4 Grounding and shielding 

The shielding of the SQUID sensor and electronics is of course a crucial issue 
for the successful design of a SQUID-based magnetometer. As mentioned in 
N2.4.21 the low-temperature parts are shielded by superconducting Pb boxes or 
Nb capillaries. Strange as it may sound, the electronics outside the cryostat posed 
in fact many more problems. The reason is that the Conductus iMAG system 
uses beautiful hybrid optical-electric cables for the communication between the 
SQUID controller and the FLL electronics on top of the cryostat, but in the cable 
that connects the FLL to the SQUID sensor, the ground is used as return line for 
the signals! This means that any ground loop involving the SQUID electronics 
will completely destroy the functionality of the system. Obviously, all the outer 
metallic parts in the system (e.g. the case of the SQUID sensor, the vacuum 
feedthrough for the cryocable, etc.) are connected to the same electrical ground, 
including the GPIB communication terminals. 

The best way to avoid ground loops in the SQUID system would be to ground 
the SQUID sensor and its cable at the cryostat (and take care that it remains a 
very clean ground), and connect the controller to an isolation transformer [H2j . 
This method is not practicable when communication with the computer via the 
GPIB bus is needed, and the same computer is connected to another instrument 
that requires a common ground with the cryostat (this is the case for the Pi- 
cowatt AVS-47 resistance bridge). The only choice is therefore to ground the 
controller at its power cord and float the whole SQUID circuitry, all the way to 
the SQUID sensor inside the cryostat and the shields of the pick-up circuit (which 
must be connected to the SQUID ground). This is already a rather cumbersome 
operation, but it's not yet sufficient. We found out that, in this configuration, 
the FLL electronics and the room-temperature cables around it are not enough 
shielded from the electromagnetic interference generated by the motor during 
the vertical movement of the dilution refrigerator. This sort of interference does 
not annihilate the functionality of the SQUID like a ground loop would do, nor 
does it simply induce an increase of the instrumental noise: the effect of bad 
shielding is that the SQUID system behaves as if it were connected to a non- 
sup erconduting, inductive circuit! It obviously took some time before we realized 
this, since the inductive behavior of the pick-up system is one of the most ex- 
pectable failures, which can be due to any weakening of the superconductivity in 
the circuit, for instance because of a bad contact on the SQUID input pads. The 
full functionality of the magnetometer was reached by enclosing the FLL and the 
room-temperature SQUID cables into an extra copper shield, grounded at the 
cryostat but separated from the SQUID ground. 



We obviously took care that the motor does not touch the cryostat ground. 
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2.4.5 Automation 

The operation of the SQUID magnetometer can be programmed in a completely 
automatic way by a Delphi software developed by W. G. J. Angenent [SH]) that 
integrates the SQUID controller, the magnet power supply, the servomotor con- 
troller, the thermometers resistance bridge and the current sources for the flux 
transformer and the pick-up heater. 

The software includes a proportional controller that drives the Keithley 220 
current source in order to null the SQUID voltage, in case the user chooses the 
measure in feedback mode rather than just recording the SQUID voltage. An- 
other option is between continuous or step-by-step vertical movement, the latter 
being more convenient for the feedback measuring mode. Furthermore, the soft- 
ware already contains routines for fitting the SQUID voltage or the feedback 
current and extract the magnetization. By programming a sequence of, for in- 
stance, several measurements at constant temperature and increasing fields, one 
can obtain a whole magnetization curve without need of any intervention of the 
user. 

A very user-friendly Windows interface is supplied, by which the user can set 
all the parameters of the measurement, control each single instrument, program 
measurement sequences and analyze the data. A screenshot of the user interface 
is shown in Fig. 12.131 

2.5 Torque magnetometer 

A torque magnetometer (torquemeter, in brief) is the ideal instrument to com- 
plement the SQUID on the high-field side, since its sensitivity grows linearly with 
the field, and it suffers none of the limitations due to the critical field of the super- 
conducting parts of a SQUID magnetometer. We constructed a torquemeter to 
be installed in a different cryostat, fitted with a 18 T Oxford Instruments NbsSn 
superconducting magnet and a smaller dilution refrigerator (Oxford Kelvinox 25), 
which was used previously for specific heat measurements Again we took 

care of preserving the compatibility with other setups, namely the calorimeter: 
the samples used for specific heat experiments can be directly recycled for torque 
magnetometery. 

2.5.1 Principles of cantilever magnetometery 

A sample with magnetic moment /2 placed in a magnetic field B experiences a 
torque T given by: 

f = px B. (2.29) 

The torque is therefore nonzero only if fl is not parallel to B, which can happen 
if the sample has an intrinsic magnetic anisotropy and the anisotropy axis is not 
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Fig. 2.13 : A screenshot of the user interface of the SQUID magnetometer 
management program. 

aligned with the field, or in a single-crystalline isotropic material, provided the 
crystal itself has a certain shape anisotropy. For a powder isotropic sample, or a 
non-oriented powder anisotropic sample, T = 0. Still, such a sample can experi- 
ence a force F = V(/2 • i?) if the magnetic field is not uniform. For instance, the 
magnetic field along the axis (z) of a magnet has, to a very good approximation, 
a quadratic dependence on the distance from the field center: 

B{z) = Bo{l-gz'^), (2.30) 

where ^ is a constant that only depends on the geometry of the magnet coils. 
The magnetic force on a powder sample is therefore F = (di?/dz)/2 = —2QzBQfl. 
If the sample is mounted on a torsion cantilever at a distance r from the rotation 
axis, it exerts a torque: 

f = - i2gzBo)rx fl, (2.31) 

which produces a displacement of the extremity of the cantilever^^^ Az = T/£, 
where £ is an elastic constant. If Az <^ z, then 2Qzr/£ = /C is just a constant 

''''in reality it is a rotation, that can be approximated as a displacement of the extremity when 
the rotation angle is very small. 
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that depends on the cantilever and its position along the magnet^^^ . The magnetic 
moment can therefore be extracted from a measurement of the displacement Az: 

Az = -ICBofi. (2.32) 

The deflection of the cantilever can be measured in several ways, including 
piezo-resistive jSS] and optical [SHI methods, but a simplest effective option is to 
shape the extremity of the cantilever as a platelet, and measure the changes in 
the capacitance formed by the platelet facing a suitable counterelectrode [HZIEHI 
Calling Co = eoA/do the capacitance of the system at rest, where A is the 
area of the platelet and do the distance at rest from the counterelectrode, the 
displacement Az produced by the magnetic force translates into a change in the 
capacitance: 

^ Az fAzV ^ 



For a paramagnetic sample and a cantilever placed above the field center, Eq. 
(|2.32|) simply implies that the sample is pushed down towards the field maximum, 
thus if the counterelectrode is below the cantilever, C increases with field. In 
particular, as long as Az <C do so that we can retain only the linear term in 
H2.33|) . we may write: 

AC = IC^Bfi. (2.34) 
do 

The scheme of a typical measurement configuration is shown in Fig. 12.141 

Notice that, at constant /i, the capacitance change (which sets the instrumen- 
tal sensitivity) grows linearly with field; this is why torque magnetometery is 
essentially a high-field technique. The magnetic moment is typically obtained by 
performing field-sweep measurements of the unbalance of a capacitance bridge, 
then dividing the results by the applied field. Obviously, one should discard the 
data around zero field, since even the tiniest noise would produce a diverging 
effect. It is not advisable to perform temperature sweeps at constant field, since 
the T-dependence of Co and do would affect the results. 



2.5.2 Design and construction 

From the previous discussion, it is clear that there are several possibilities to 
increase the sensitivity of a torquemeter: (i) use a soft cantilever (small £); 
(ii) increase the distance r of the sample from the rotation axis; (iii) build a 
capacitor with large area and very small distance between the plates; (iv) place 
the torquemeter far from the field center. The fourth option must be taken with 

'^'Notice that K, has a sign: fC > above the field center, and fC < below. 




Fig. 2.14 ■' Scheme of a typical configuration for capacitive detection of the 
magnetic force on a paramagnetic isotropic sample in a field gradient. The mag- 
nitude of the field inhomogeneity is largely exaggerated. 



care, and in our case z is determined by the position of the three experimental 
slots available in the refrigerator |40j: one slot is precisely at the center, the other 
two are ~ 25 mm above and below. The longitudinal inhomogeneity factor [cf. 
Eq. (|2.:-{()|) ] of the NbsSn magnet used here is ^ ~ 2.8 • 10~^ T/mm^, which means 
that a torquemeter in the upper slot experiences a field gradient AB /dz ~ 0.035-Bo 
T/mm. Contrary to the setup described in H2.11 in this refrigerator the sample 
is in vacuum, thus it must be thermalized by the cantilever itself. 

The system, shown in Fig. 12.151 lies on an Araldite base, screwed laterally 
between two wide copper slabs that constitute the cold finger connected to the 
mixing chamber (not shown). The torquemeter is designed as a symmetric tor- 
sion balance, to avoid the possible contribution of the platelet to the magnetic 
force. The 20 x 4 mm central strip is held in position by two thin arms, which 
constitute the elastic twistable element. This cantilever is obtained by spark ero- 
sion on a 50 ^um thick copper-beryllium (CuBe) foil, subsequently annealed at 
315 °C for two hours in order to eliminate the mechanical stress and obtain a 
perfectly flat foil. The sample, mixed with Apiezon N grease for thermal contact, 
is placed at one extremity of the balance, which constitutes the upper plate of 
the capacitor. The lower plate is obtained from a square copper pad on a printed 
circuit board. The CuBe foil is held by two brass clamps screwed onto the base, 
with the interposition of 25 ;um thick copper spacers. The clamps have transver- 
sal grooves in correspondence with the torsion arms of the cantilever. To avoid 
the premature touch of the cantilever on the counterelectrode, we typically put 
two spacers on top of each other. Pads with the same shape as the spacers are 
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(a) Brass clamp for thermal post 

(b) Brass clamp for radiation shield 

(c) Copper radiation shield 

(d) Brass clamps for CuBe foil 

(e) Transversal grooves 

(f) Copper strip for thermal contact 

(g) Sample 

(h) CuBe foil 

(i) Torsion cantilever 
(j) Cigarette paper impregnated 

with GE varnish 
(k) Copper spacers 
(I) Soldering pads for coax cables 
(m) Counterelectrode 
(n) Printed circuit board 
(o) Thermal anchoring post 
(p) Copper strip to ground the back 

of the printed circuit board 
(q) Araldite base 



Fig. 2.15 : Drawing of the disassembled torque magnetometer. 
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drawn on the printed circuit board; their function is to thermahze the system 
and to contact electrically the CuBe foil to a coax cable soldered on the board. 
Since the cantilever constitutes one of the plates of the capacitor and may not be 
connected to the electrical ground, a L-shaped copper strip is screwed on top of 
a thermal anchoring post, obtained from teh printed circuit board. A cigarette 
paper impregnated with GE varnish is inserted between the copper strip and 
the anchoring post, to provide good thermal contact while maintaining electrical 
insulation. The copper strip is then screwed laterally between the Araldite base 
and the cold finger. To shield the sample from thermal radiation and electromag- 
netic noise we cover the system with a copper box, coated with thin Kapton tape, 
screwed laterally to the cold finger and therefore also connected to the electrical 
ground. To complete the shielding, a small copper strip connects the front to the 
back side of the printed circuit board, which is completely copper-plated. 

Two coaxial cables are soldered to copper pads on the printed circuit board, 
one connected to the counterelectrode, the other to the cantilever. We can there- 
fore measure the capacitance between cantilever and counterelectrode by con- 
necting the coaxials to a General Radio 1615-A capacitance bridge. A Stanford 
SR830 lock-in amplifier provides the 5 V, ~ 1 kHz excitation, and measures 
the unbalance of the bridge. The measurements are controlled by a Lab View 
program developed by W. G. J. Angenent. Before starting a field-sweep, the 
temperature is set to the desired value, then the voltage due to the unbalance of 
the capacitance bridge at zero-field is accurately measured, in order to subtract it 
from the data afterwards. While sweeping the magnetic field, typically at a rate 
dB/dt = 0.05 — 0.1 T/min, the lock-in voltage Vl is measured and divided by 
B after subtraction of the zero- field offset. This yields the sample magnetization 
M{B) oc fJ,{B), since Vl oc AC [cf. Eq. (|2.34j) ]. The exact conversion factor 
between the lock-in voltage and the magnetization is very difficult to obtain in a 
reliable and reproducible way^*^; we shall therefore express M in the "electrical 
units" /iV/T. 

2.5.3 Performance 

An example of the performance of our torquemeter is given by the experiment on 
a ~ 0.25 mg Cerium Magnesium Nitrate (CMN) sample, which is an ideal Curie 
paramagnet down to at least T ~ 10 mK. 

The measured magnetization, shown in Fig. I2.16r al. exhibits the expected 
Brillouin behavior up to a certain threshold, where the non-linear terms in (|2.33|) 
start to be important. The condition of linear response AC oc Az is obeyed below 
a certain maximum displacement, i.e. a horizontal line in the V^ — B plane, which 
translates into a hyperbola in the M — B plane. 



'^'it depends for instance on the precise position of the sample on the platelet (r), the distance 
at rest between the capacitor plates (do), etc. 
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-Figi. 2.16 : (a) Magnetization of a 0.25 mg CMN sample, at the indicated 
temperatures (in Kelvin). Dashed line: (emprical) threshold of the linear response 
region, (h) Static susceptibility as a function of the inverse temperature, fitted by 
the Curie law. 



By fitting the initial linear part of M{B) at different temperatures, we obtain 
the static susceptibility x{T). As shown in Fig. I2.16f b'). the measured susceptibil- 
ity obeys indeed the Curie law, x ^ ^/T, and confirms the proper thermalization 
of the sample. 

Knowing that the magnetization of CMN is ~ 3.3 emu/g at -B = 1 T and T = 1 
K (as deduced from Fig. 14 in we can estimate that the linear response 

region is roughly ■ B <2 ■ 10^'^ emu-T. The instrumental noise, typically 10 
nV, translates into an equivalent noise of ~ 10~^ emu in the magnetic moment 
at = 1 T. By increasing the field to e.g. 10 T, another order of magnitude 
is gained in the equivalent magnetic noise, but the maximum measurable signal 
is also 10 times smaller. Therefore, the appropriate amount of sample must be 
chosen with care, as a function of the range of B and /x of interest. 
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Theoretical aspects of molecular 

magnetism 



This chapter contains a selection of theoretical approaches necessary to under- 
stand the physics of molecular magnets and their interaction with the environ- 
ment. After an introduction to the relevance of this subject for the extrapolation 
of quantum mechanics to the large scale, we proceed to a "top-down" physical 
description, starting from the model for a single giant spin in its crystalline en- 
vironment, then illustrating the effect of a magnetic field and the coupling to 
nuclear spins and the mutual coupling via dipolar fields. We then discuss the 
standard theory of the nuclear spin dynamics, mainly to show that such a treat- 
ment is no longer adequate in the presence of macroscopic quantum tunneling: 
the Prokof 'ev-Stamp theory is the necessary approach to a unified description of 
a giant quantum spin coupled to a spin bath. One of the most crucial outcomes 
of the present work is that even the Prokof'ev-Stamp theory needs an extension 
in order to account for our experimental results presented in chapter IV. We shall 
therefore come back to this issue in an extra theoretical section at the end of that 
chapter. 



3.1 Macroscopic quantum phenomena 

Besides the perspective of possible applications (magnetic recording media, spin 
qubits, etc.), single-molecule magnets are very attractive systems to study the 
observability of quantum phenomena at the macroscopic scale. The motivation 
behind this interest dates back to the formulation of the Schrodinger cat para- 
dox [1] and the so-called "measurement problem" in quantum mechanics. In fact, 
the "weak point" of quantum mechanics in its present formulation is the rather 
artificial distinction between the (microscopic) quantum system and the (macro- 
scopic) measurement apparatus, which is supposed to obey the laws of classical 
physics, and whose action forces the wave function of the quantum system to 
be projected (collapse) onto one of its eigenstates. Nothing is specified about 
where the border between a microscopic (quantum) and a macroscopic (clas- 
sical) system lies, neither how the collapse of the wave function actually takes 
place. The proposals to solve these problems range from the most pragmatic idea 
of decoherence (Hj, all the way to theories that effectively add extra postulates 
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to the standard quantum mechanics |42ll43j . Before one even attempts to use 
the experimental observations on SMMs to elucidate the fundamental problems 
mentioned above, it is essential to quantify the degree of "Schrodinger cattiness" 
of the quantum states being observed. Here we briefly review the approach of 
Leggett |441I1U| . who usefully introduces two parameters to quantify the macro- 
scopicity of a quantum superposition of states: the "extensive difference" and the 
"disconnectivity" . 

3.1.1 Extensive difference 

Let us consider a system in a quantum superposition of states, i.e. whose wave 
function l^/") can be expressed as a linear combination of two wave functions: 



with the assumption that \ipa) and l-i/^f,) represent states of the system which are 
by some reasonable criterion "classically different" . We may therefore character- 
ize the two branches of the superposition by some extensive quantities {i} (e.g. 
charge, magnetic moment, position of the center of mass, etc.) which should 
be considerably different in the states |V'a) and iV'b). Next, we define Aj as the 
difference between the values of the extensive quantity i in the two branches of 
the superposition, divided by a "typical value" of i at the atomic scale (e.g. 1 
electron charge, 1 Bohr magneton, 1 Angstrom, etc.). The "extensive difference" 
A of the quantum superposition is then the maximum value assumed among the 



This sounds at first as a very natural definition of macroscopic distinctness, 
but in fact it's easy to find examples where A is a very large number, although 
the system is not what one would like to call "macroscopic" . Think of a neutron 
passing through an interferometer with arms e.g. 10 cm apart: one finds A ~ 10^, 
although this is clearly not the sort of examples we are looking for as a challenge 
to the interpretations of quantum mechanics. 

3.1.2 Disconnectivity 

Leggett introduces therefore another parameter, P, called "disconnectivity". The 
precise definition can be found in Ref. [H], but for our purpose it is sufficient to 
say that T> is the number of particles that have substantially different behavior in 
the states \il)a) and \'4'b)- More precisely, V is the number of particle correlations 
that must be measured in order to distinguish the coherent superposition 
from a statistical mixture of iV'a) and \ipb)- For instance, a system of N identical 
particles prepared in a state like 



(3.1) 



Ai's. 



1^) 



N 



(3.2) 
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has D = N, since particles are simultaneously superimposed, whereas a state 



has T> = 1, being just the product of one-particle superpositions. The above- 
mentioned neutron in a diffractometer has obviously I? = 1, thereby reestab- 
lishing the fact that it does not constitute a "true" macroscopic quantum state. 
An interesting example is the one of a superconductor: it's easy to realize that 
V = 2, since the standard BCS wave function is a product of two-particle corre- 
lations. The "charge qubit" !45^, which is one of the most exciting developments 
in solid-state quantum devices of the last few years, again has T> = 2 since it is 
based on the superposition of a Cooper pair (2 particles) being inside or outside 
a superconducting island; the fact that the device itself has a size ~ 1 //m does 
not automatically make it relevant for the problem of quantum mechanics at the 
large scale! 

The observation of high-D quantum superpositions has been realized only very 
recently: for instance, by diffraction of Cqo molecules !46_: a value of P = 1048 has 
been achieved. Even higher values have been reached by the "flux-qubit" |17p47j . 
which features the superposition of counter-rotating supercurrents in a SQUID 
loop: in this case all the current-carrying Cooper pairs (i.e. those within the 
London penetration depth from the surface) are behaving differently in the two 
branches, yielding P > 10^ ! 

3.1.3 Macroscopicity of a single-molecule magnet 

We can now apply the concepts defined above to the case of a SMM. In general, we 
will be interested in phenomena that arise from the quantum superposition of two 
different projections of the total spin along the z axis . For example, for Mni2-ac 
the total spin S" = 10 can give rise to a total magnetic moment of 20 fiB pointing 
along +zor —z. A superposition of such two states would be characterized by an 
extensive difference A = 40. As for the disconnectivity, we anticipate (see ij3.2.4() 
that the total spin of the molecule is the result of the ferrimagnetic coupling of 
8 Mn^"*" ions with spin s = 2, i.e. 4 electron spins per ion, and 4 Mn'^+ ions with 
s = 3/2 (3 electrons). In total, the number of electrons having different states 
when S is along +z or— zisP = 4x8 + 3x4 = 44. Although it doesn't reach 
the macroscopicity of a fullerene molecule or a SQUID qubit, a SMM is therefore 
substantially more macroscopic that an atomic-size quantum system. 

3.2 Effective Spin Hamiltonian 

The magnetic properties of the single-molecule magnets investigated in this the- 
sis are determined in the first place by the net electron spin of each single ion, 
arising from the partial filling of 3d shells. Because of the strong crystal field 




(3.3) 
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effects, the electron angular momentum is quenched, so that no orbital contribu- 
tion needs to be included. Furthermore, the electron spins within a cluster are 
magnetically coupled by superexchange interactions, typically due to the orbital 
overlap through oxygen bridges. In this way, one obtains a magnetic ground state 
which can be characterized by a high value of the total cluster spin S, like in the 
case of Mni2-ac and Mng, and is separated from other spin states by an energy 
determined by the strength of the intracluster superexchange couplings. As long 
as the temperature is kept much lower than the energy separation between the 
ground and the first excited total spin state, it is very convenient to adopt the 
effective spin Hamiltonian description, which means that the cluster is treated as 
an object characterized by just the number of energy levels compatible with the 
total spin ground state, i.e. 25 + 1. The resulting energy levels can be further 
split by the effect of crystal field anisotropy or couplings with magnetic fields. A 
typical effective spin Hamiltonian to describe the crystal field anisotropy is the 
following: 

HcFA = -DS"^, - BSt + E{Sl - Sl) + C{SX + St). (3.4) 

In the next subsections we shall discuss which predictions about the physical 
behavior of a single-molecule magnet can be made on basis of the above Hamil- 
tonian. 

3.2.1 Superparamagnetic blocking 

The terms —DS^ and —BS'^ (with \B\ <^ \D\, typically) represent uniaxial 
anisotropies: if -D > 0, then z is the easy axis of magnetization. In manganese- 
based clusters the uniaxial anisotropy typically arises from the Jahn- Teller dis- 
tortion of the coordination octahedra at Mn^+ sites; the total cluster anisotropy 
is then obtained as the vector sum of the single- ion anisotropies 0H| ■ Notice that 
by considering only these two terms, TLcYA commutes with the Sz operator, thus 
the eigenstates \m) of Sz are also exact eigenstates of TLcFA- The energy levels 
scheme, as shown in the left panel of Fig. 13.11 can be regarded as a set of doublets 
of degenerate states, | + m) and | — m), with energies -E+m = each state 

being localized on the left or right side of the anisotropy barrier of height U . At 
temperatures comparable with the barrier height, thermally activated transitions 
from one side to the other are very fast, and the molecule behaves as a high-spin 
paramagnet. Since the high spin results from several ions behaving collectively 
as a single-domain particle, the system is called superparamagnet [IHl- The re- 
laxation rate by thermal activation, r~^, depends on temperature according to 
the Arrhenius law: 

T-^ = T^^eM-U/kBT). (3.5) 

At sufficiently low temperatures, may become exceedingly long and give rise 
to hysteresis in the magnetization loops ^ and the appearance of a frequency- 
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dependent dissipation peak in the ac-susceptibility ||50J. The temperature, Tb, 
at which the cluster spin can no longer follow an external driving field is called 
blocking temperature, and obviously depends on the timescale relevant for the 
experiment in question. 



3.2.2 Spin quantum tunneling 

The quantum aspects of SMMs are encountered by considering non-diagonal 
anisotropy terms like E{S'^ — Sy), which describes a hard-axis anisotropy caused 
by a rhombic distortion, and C(S^ -|- St.) which is the lowest-order non-diagonal 
term allowed by tetragonal symmetry. These terms do not commute with Sz, 
so the eigenstates j-^) of the complete TCcFA are no longer a set of doublets of 
localized states (Fig. 13. Ij) . Expressing \ip) in the basis of the eigenstates |m) of 

\^p) = '^Cm\m), m = -S... + S, (3.6) 

m 

one finds that each doublet now consists of a state \tps) with symmetric coeffi- 
cients, c+m = C-m, and an antisymmetric state IV'a) with c+m = —C-m, separated 
by an energy gap A = Ea — Eg- In particular for the ground doublet, one finds 
to a very good approximation: 

m = ^{\ + s) + \-s)) 

m = -^i\ + S)-\-S)). (3.7) 

This means that a state localized on one side of the barrier, e.g. | +»S'), must now 
be expressed as a superposition of the actual eigenstates: 

\+S) = ^{\i;s) + \i^A)). (3.8) 

It is clear from basic quantum mechanics jSJ that if one would prepare the system 
at t = in such a state, {iplt = 0)) = 1+5"), the time evolution should obey the 
Schrodinger equation 

\m) = -^e-*^^*(|V^s)e+*^* + |VA)e-^*)- (3-9) 

This is equivalent to having coherent oscillations of the cluster spin between the 
states I + S) and | — 5) with frequency ujt = A/?i, which implies that the spin 
is tunneling back and forth through the anisotropy barrier. For this reason the 
energy gap A is called tunneling splitting. 
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Fig. 3. 1 : Sketch of the energy levels scheme for the Hamiltonian when 
considering only the diagonal terms (left) or including the non- diagonal ones 
(right). The magnitude of the tunneling splittings is largely exaggerated and even 
not to relative scale. 

3.2.3 Coupling with external fields: coherent and incoherent 
tunneling 

It must be stressed that, in real SMMs, the tunnehng sphtting (particularly the 
one of the ground doublet) is a very small quantity when produced by non- 
diagonal anisotropy terms only. As we shall discuss below, the coupling with 
environmental degrees of freedom is many orders of magnitude larger than A, 
which basically eliminates the possibility to observe coherent tunneling oscilla- 
tions as described in Eq. One of the great advantages of single-molecules 
magnets as systems to study macroscopic quantum effects, is that the parameters 
of the spin Hamiltonian, and thereby the whole physical properties, can be easily 
manipulated by applying an external magnetic field B. This introduces an extra 
term in the effective spin Hamiltonian: 

n = ncFA-gf^BS-B. (3.10) 

By choosing the direction of B with respect to the anisotropy axis z, we have 
the freedom to either introduce extra non-diagonal terms (B _L z) that enhance 
the quantum behavior, or to destroy the symmetry of the energy level scheme 
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Fig. 3.2 : If the bias ^{t) crosses the tunneling resonance, the cluster spin may 
tunnel through the anisotropy barrier with a probability given by Eq. \3.11]) . The 
drawing on the right represents the time evolution of a doublet of electron spin 
levels while crossing the resonance. 



(B II z) and thereby reduce the coupling between states on opposite sides of the 
anisotropy barrier. 

Very important for the present discussion is the case where a fongitudinal 
magnetic field Bz{t) fluctuates in time and produces a bias, ^(t) = 2gfiBBz{t)S, 
on the cluster spin doublets. If ^{t) spans a range much larger than A but 
crosses several times through the tunneling resonance, one finds that at each 
crossing there is a nonzero probability for the cluster spin to tunnel through the 
barrier, but there is no correlation between subsequent tunneling events; therefore 
this is called "incoherent tunneling". The probability Plz of a single tunneling 
event upon crossing through the tunneling resonance can be calculated with the 
Landau-Zener formula ^52) : 

A sketch of a Landau-Zener transition is given in Fig. 13.21 This phenomenon has 
been successfully exploited as a way to extract the tunneling splitting by 
sweeping an external magnetic field through the tunneling resonance. As will be 
discussed below, the Landau-Zener formalism can also be used to describe the 
coupling between a single cluster spin and its environment, since very often the 
effect of environmental fluctuations can be treated as an effective magnetic field 
that couples to the cluster spin. 

In the presence of incoherent tunneling fluctuations, it is essential to distin- 
guish between the average time interval tt between subsequent uncorrelated tun- 
neling events, and the so-called tunneling traversal time rtr, which is in our case 
the timescale for the reversal of the cluster spin. A detailed analysis of the tun- 
neling traversal time constitutes a topic on itself |55 ll56p57j . but for our purposes 
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Fig. 3.3 : Illustration of the relationship between tunneling traversal time rtr, 
bounce frequency and tunneling interval tt for incoherent tunneling of the 
cluster spin. 

it is sufficient to estimate the order of magnitude of rtr- In the general framework 
of the semiclassical instanton technique [6|i58j , rtr is related to the inverse of the 
"bounce frequency" or "attempt frequency" VIq, i.e. the frequency of the small 
oscillations at the bottom of each potential well. In SMMs H^q is of the order of 
the energy separation between the two lowest electron spin doublets, as shown in 
Fig. 13.31 For instance in Mni2-ac we have /iOq ~ 10 K, thus r^r ~ ~ 10~^^ 
s. This means that the tunneling events take place in a virtually instantaneous 
way, as compared to the timescale of all other relevant phenomena. 

In order to push the system into a regime where coherent tunneling may 
take place, it is necessary to increase A above the strength of the coupling with 
environmental degrees of freedom. This can be done by applying a magnetic field 
B± perpendicular to the anisotropy axis. Although sophisticated calculations of 
A(Sx) have been proposed [^511001) it is often more practical to obtain A(i?x) 
directly from the numerical diagonalization of the spin Hamiltonian. An example 
of the resulting A(S_|_) in the case of Mni2-ac is shown in Fig. 13.41 A simplified 
expression for A(i?x), based on the Hamiltonian TL = —DS'^—gfiBB±-S, has been 
obtained in the limit g^BB±_ <^ DS and 5 ^ 1 by Korenblit and Shender |61j : 



where e is the base of the natural logarithm. The more general expression for the 
splitting Am of any doublet of levels | + m), | — m) is 



It clearly appears that the simple application of a perpendicular field has an 
enormous influence on A, which allows to study the quantum dynamics of SMMs 
in different regimes. 




(3.12) 




(3.13) 
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Fig. 3.4 ■' Perpendicular field dependence of the tunneling splitting of the ground 
doublet in Mni2-ac, as obtained by numerical diagonalization of the Hamiltonian 
i3.14[ )- The two lines refer to B± parallel to the x or y crystallographic axes, or 
oriented at 45° between them. 



From Eq. (|3.13|) it is also clear that, since g^BB±_ <^ DS, the tunneling split- 
ting increases drastically (exponentially, in fact) when higher excited doublets 
(i.e. with lower |m|) are considered '^^^ Therefore, in the presence of thermal 
excitations, tunneling can more easily proceed through excited doublets than 
through the ground state. This remains true in zero external field, although the 
tunneling splitting is produced by non-diagonal anisotropy terms. The details 
of the mechanism of thermally-assisted tunneling have generated a vast litera- 
ture [63h64. 65,66,67]. since a large majority of experiments on SMMs have been 
carried out in the temperature regime T > 1.5 K, where tunneling through ex- 
cited states is essential. The transition to pure ground state tunneling happens 
indeed at temperatures that, at a first sight, appear surprisingly low. The point 
is that, although the Boltzmann factors for the population of excited doublets 
are exponentially small, grows exponentially upon excitations to higher dou- 
blets ()3.13|) . so the competition between the two exponents cannot be treated in 
a trivial way [HHI- Since the present work is focused on the ultra-low tempera- 
ture regime, the topic of thermally assisted tunneling will be touched upon only 
marginally. 

'^'For the case of a SMM in zero external field, the role of B± in is taken by the non- 

diagonal crystal-field anisotropy terms and by the transverse component of the dipolar field from 
neighboring molecules. By condensing such effects into an "equivalent B±" of constant value, 
it is easy to see that oc /(m) • (const.)^"^ 
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Fig. 3.5 : (a) Structure of the Mni2-ac cluster, with the labelling of the three 
inequivalent Mn sites as described in the text, (b) Energy level scheme for the 
electron spin as obtained from the Hamiltonian i3.14\ ), retaining only the terms 
diagonal is Sz- The non-diagonal terms allow transitions between states on oppo- 
site sides of the anisotropy barrier by means of quantum tunneling (QT). In the 
presence of intrawell transitions induced by thermal excitation, thermally assisted 
quantum tunneling (Th-A T) between excited doublets can also take place . 



3.2.4 Parameters for Mni2-ac 

The structure of the [Mni20i2(02CMe)i6(H20)4] (Mni2-ac) cluster jHH (Fig. 
13. 5|) contains a core of 4 Mn'^^ ions with electron spin s = 3/2, which we shall 
denote as Mn(^\ and 8 Mn'^"*" ions (s = 2) on two inequivalent crystallographic 
sites, Mn*^'^) and Mn*-^-*. The electron spins are coupled by mutual superexchange 
interactions, the strongest being the antiferromagnetic interaction between Mn^^^ 
and Mn^^) j7Uj . The molecules crystallize in a tetragonal structure with lattice 
parameters a = b = 17.319 Aand c = 12.388 A. Numerical studies have shown 
the influence of the magnitude and sign of the exchange constants on the precise 
energy level structure of the cluster and an attempt has been made to verify 
the assignment of the exchange interactions by a spin-wave analysis of the nuclear 
spin-lattice relaxation rate in the thermally activated regime jJS]. The ground 
state of the molecule has a total electron spin S* = 10 and is separated from 
the first excited manifold 5 = 9 by ~ 58 K, according to neutron scattering 
experiments [731 . For the temperature range of interest in the present work 
(r < 2 K), we may therefore describe the electron spin of the cluster by means 
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of the spin Hamiltonian in the S" = 10 manifold: 

n = -DSl - BSt + C{SX + St) + ^bB • g • 5. (3.14) 

We adopt the parameter values D = 0.548 K, B = 1.17 x 10"'^ K and C = 
2.2 X 10~^ K as obtained by neutron scattering data [71], and for the g tensor the 
values (7|| = 1.93 and gf^ = 1-96 from high-frequency EPR^^^ [75j. The uniaxial 
anisotropy terms —DS^ and —BS^ can be attributed to the single-ion anisotropy 
of the Mn3+ ions [T^, which is due to the crystal field effects resulting in the 
Jahn- Teller distortions of the coordination octahedra, where the elongation axes 
are approximately parallel to the c-axis of the crystal. The energy levels scheme 
is then a series of doublets of degenerate states which are separated by a barrier 
with a total height DS^ + BS'^ ~ 66.6 K. The non-diagonal anisotropy term 
C{SX + 5*1 ) arises from the fourfold S4 point symmetry of the molecule, and is 
the lowest order term that would allow the quantum tunneling of the electron 
spin between states Sz = =bm. An external magnetic field B^ parallel to the 
anisotropy axis would introduce an extra term g^^fisBzSz in the Hamiltonian, 
thereby destroying the symmetry of the energy level scheme; at the particular 
values of field S™™" = [n/{g\\iJ.B)][D + B{m'^ + m''^)], with n an integer —S < n < 
S, the states S"^ = m on one side of the barrier come in resonance with the states 
Sz = m' = n — m on the opposite side. In the presence of nondiagonal terms 
in the spin Hamiltonian, one then expects the occurrence of quantum tunneling 
of the magnetization, as indeed observed in the experiments 0. The selection 
rules imposed by the term (7(5"^ -1- S"!) imply that resonant tunneling should 
be observed only for n a multiple of 4, in contrast with the observation that 
steps of comparable height appear in the magnetic hysteresis loops for all values 
of n. There is now solid experimental evidence |77|l78j for the prediction jJO] 
that a disorder in the acetic acid of crystallization is present and gives rise to 
six different isomers of Mni2 cluster, four of which have symmetry lower than 
tetragonal and therefore possess non-diagonal anisotropy terms that would allow 
tunneling transitions for any even value of n. Furthermore, to obtain tunneling 
resonances also at odd values of ri, it is necessary to introduce in (|3.14|) the 
effect of dipolar fields originating from neighboring molecules and the hyperfine 
couplings with the nuclear spins, as will be described below. 

A peculiarity of the Mni2-ac system is the presence in every real sample 
of fast-relaxing molecules (FRMs) 80^, i.e. clusters characterized by a lower 
anisotropy barrier and a much faster relaxation rate, as observed for instance by 
ac-susceptibility and magnetization measurements ISSj- It has been recog- 
nized that such FRMs originate from Jahn- Teller isomerism |53], i.e. the presence 
in the molecule of one or two Mn'^"^ sites where the elongated Jahn- Teller axis 



' 'The anisotropy parameters obtained by EPR seem to depend on the magnetic field range 
used in the experiment [751176) . whereas neutron scattering is a zero-field experiment, and does 
not require assumptions on the g tensor in order to fit the data. 
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Fig. 3.6 : (a) Crystal structure of Mn^OiBr4^(Et2dhm)Q ISI^ . (h) Detail of the 
symmetric octahedral core, containing six ferromagnetically coupled Mr?^ ions, 
yielding a total spin 5 = 12 for this molecular superparamagnetic particle. 

points in a direction roughly perpendicular instead of parallel to the crystalline 
c-axis. This results in the reduction of the anisotropy barrier to 35 or 15 K in 
case of one or two flipped Jahn- Teller axes, respectively , and presumably in 
an increase of the non-diagonal terms in the spin Hamiltonian as well. Further- 
more, the anisotropy axis z of the whole molecules no longer coincides with the 
crystallographic c-axis, but deviates e.g. by ~ 10° in the molecules with 35 K bar- 
rier |S2]. The Jahn- Teller isomerism is very different from the above-mentioned 
effect of disorder in solvent molecules, and produces much more important effects. 
As shall become clear in chapter IV, the presence of the FRMs is essential for the 
interpretation of our experimental data. 

3.2.5 Parameters for Mng 

The molecular core of Mn604Br4(Et2dbm)6, hereafter abbreviated as Mng j85j . 
is a highly symmetric octahedron of Mn^+ ions (with spin s = 2) that are fer- 
romagnetically coupled via strong intra-cluster superexchange interactions (Fig. 
13. 6|) . Accordingly, the ground state is a S = 12 multiplet and the energy of the 
nearest excited state is approximately 150 K higher. The unit cell is monoclinic, 
with space group Pc and contains 4 molecules^^^ bound together only by Van 

'''•'We have used the crystallographic data for a related complex Mne04Cl4(Et2dbm)6 |85|. in 
which the Br~ ions are replaced by Cl~ . 
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der Waals forces. Inter-cluster superexchange is therefore negligible and the only 
relevant interaction interaction between cluster spins is the dipolar one. 

To determine the spin Hamiltonian of the S = 12 manifold, we notice that 
the anisotropy of the sing le Mn3+ ions may be very strong due to Jahn- Teller 
effects. On the other hand, due to the highly symmetric octahedral structure of 
the cluster, the addition of the single-ion anisotropies may lead to a vanishingly 
small net cluster anisotropy, no matter how large is the zero-field splitting of the 
constituting atoms |1H]. Indeed, the analysis of magnetization data on basis of 
the simple Hamiltonian: 

n = -DSl- q^ibB ■ S. (3.15) 

indicated an upper limit \D\ < 0.01-fC for the total cluster anisotropy [HS]) with 
g ~ 1.94 — 1.98. An independent estimate of D was obtained by high-frequency 
ESR data taken in the range 95 - 380 GHz by J. Krzystek (NHMFL Tallahassee, 
USA) (SB]. Owing to the combination of very small D and large S", as well as the 
presence of a signal at g = 2.00 arising from a minute amount of Mn^+ impurity 
(often seen in ESR of Mn'^^ compounds), the interpretation of the spectra was 
not fully conclusive. Nevertheless, signals with a clearly visible structure on the 
low-field end of the spectra could be obtained. It could be identified as fine 
structure originating from zero-field splitting (ZFS), since it was independent of 
field and frequency. Simulations of the spectra performed using Eq. 1)3. 15() agree 
well with the experiment taking \D\/kB ~ 0.03 or ~ 0.05 K, depending on the 
sign of D (which could not be unequivocally determined). Although a smaller 
rhombic component could be present, the data do not justify a more elaborate 
fitting. 

The Mng cluster constitutes therefore a practical example of how the anisotropy 
parameters of high-spin molecules can be tuned by simply modifying the crystal 
structure, while maintaining the same constituents as other clusters having very 
different properties (cf. Mni2-ac in p.2.4|l . 

3.3 Hyperfine couplings 

The coupling between electron and nuclear spins in a single ion can be generally 
expressed in the form: 

"^hyp = A • S = -7Ar^/ • Shyp (3.16) 

where I is the nuclear spin, s is the electron spin of the single ion, A is the 
hyperfine coupling tensor, 77V is the nuclear gyromagnetic ratio and -Bhyp is the 
total hyperfine field at the nucleus. The hyperfine tensor can contain several 
contributions depending on the particular nuclear site under consideration. For 
instance, A for the ^^Mn nuclei in a Mn^"*" ion is a diagonal tensor: the only 
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Fig. 3.7 : ^^Mn NMR spectra in zero field, with assignment of each line to the 
corresponding site in the molecule (adapted from \81^ ). 

contribution comes from the the isotropic Fermi contact term, arising from the 
core polarization of inner s electrons caused by Sd electron spins. For the 
nuclei, located at the ligands of the molecules, the main contribution is instead 
arising from the dipolar field produced by the electron spins. The resulting -Bhyp 
is then typically much smaller than for ^^Mn. 

3.3.1 Hyperfine fields on the ^^Mn nuclei in Mni2-ac 

Thanks to the high anisotropy barrier of Mni2-ac, it is possible to detect the 
Larmor precession of ^^Mn nuclear spins even in the absence of an external field, 
by exploiting the local hyperfine field from the electron spins which, below the 
blocking temperature Tb ~ 3 K of the superparamagnetic clusters, is static on 
the timescale of an NMR experiment. Goto et al. j88j were the first to find 
that the ^^Mn NMR spectrum in zero field consists of three quadrupolar-split 
lines centered around 230, 280 and 365 MHz. The hyperfine fields in each of 
the three inequivalent Mn sites are directly obtained as -Bhyp = '^t^^Mn-, where 
7Ar/27r = 10.57 MHz/T is the gyromagnetic ratio of the ^^Mn nuclei, yielding 
= 21.8 T, = 26.5 T and sgj, = 34.5 T. Subsequently, a detailed 

analysis of the hyperfine couplings has been reported by Kubo et al. [HZI; who 
were able to assign these three NMR lines to specific Mn sites in the molecule. 
The line centered at v^^'^ ~ 230 MHz corresponds to the nuclei in Mn^"*" ions, and 
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is characterized by a relatively small quadrupolar splitting — 0.72 MHz. 

The direction of -Bhyp is antiparallel to the Mn^+ electron spin, thus it coincides 
with the direction of the total spin S" = 10 of the molecule. In the Mn^+ ions, 
the hyperfine interaction contains also a dipolar term, which makes the coupling 
tensor anisotropic. Furthermore, -Bhyp is not exactly antiparallel to the electron 
spin, but there's a small canting angle j87j . Accounting for the effects of a dipolar 
contribution to the hyperfine field in Mn^^, it was then possible to assign the line 
at ~ 280 MHz to the Mn^^) site and the line i/^^) ~ 365 MHz to the Mn^^) 
site jHZI (cf. Fig. 13. 7|) . The quadrupolar splittings of these lines are Az^g ~ 4.3 

MHz and Ai/J^ ~ 2.9 MHz 

The structure mentioned above is refiected in the effect on the NMR fre- 
quencies of an external field applied along the anisotropy axis. The total 
field at the nuclear site becomes -Btot = -^hyp + B^, in particular for the Mn^^^ 
site (but approximately also for the other sites), since -Bhyp is parallel to B^, 
iy{Bz) = 7Ar(-Bhyp + Bz) = J^(0) +jnBz exhibits a slope |89|l87] which is given by 
the gyromagnetic ratio of ^^Mn. 



3.3.2 NMR spectra in perpendicular field 

If the magnetic field is applied perpendicular to the anisotropy axis, the situation 
is complicated by the fact that the electron spins, and as a consequence also 
the hyperfine fields, tend to cant towards the direction of the external field 
Recently, Furukawa et al. |90j have analyzed the ^^Mn NMR spectrum of Mni2-ac 
in perpendicular fields by assuming that the canting angle 9 between the direction 
of the electron spins and the c-axis of the molecules is given by sinO = M_\_/Ms, 
where Mg is the saturation magnetization of the = 10 ground state, and M± 
is the component of the magnetization along the field direction. At the very low 
temperatures involved in the present work the contribution of excited electron 
spin states to the total magnetization is negligible, and the above approach is 
equivalent to expressing the canting angle as^^^: 

sin^ = ^, (3.17) 

where {Sx) = (V'gI'S'xIV'g) is the expectation value of the x-component of the 
spin in the ground state IV'g)) to be obtained by numerical diagonalization of 
the Hamiltonian 1)3.14(1 . There's no difficulty in allowing the field B to lie at an 
arbitrary angle 9b with respect to z, which helps accounting for the possibility 
of misalignments between external field an crystal axes. Assuming that the hy- 
perfine field -Bhyp has always the same strength as deduced from the zero-field 
spectra and a direction given by the canting angle calculated above, the total 



Here and in the following we assume that the field is always applied in the xz plane. 
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Fig. 3.8 : Sketch of the orientation of hyperfine fields in Mn^'^ (^hyp) ^^'^ Mr?^ 
(-^hyp) ions, relative to the crystal axes x and z of the molecule, in the presence 
of an external magnetic field B at an angle 6b from the z axis. The total field at 
the nuclear site, Btot, which determines the NMR frequency, is obtained as vector 
sum o/i?hyp and B. 

field at the nuclear site is obtained as: 

Btot = ^(Bhyp smO + B smOB? + (^hyp cosO + B cos ObY (3.18) 
for Mn^+ ions and 




(3.19) 



for the Mn^ ions, for which the direction of -Bhyp is opposite to the total magnetic 
moment of the cluster (Fig. 13.8(1 . The ^^Mn Larmor frequency can be directly 
calculated as: 

v{B) = ^Btot, (3.20) 

yielding the field dependence of the Larmor frequencies shown in Fig. 13.91 

The condition B^yp(B) = i?hyp(0) for any value of B, which is contained in 
Eq. H3.17() and used to derive Eqs. (|3.18|) and (|3.19() . is actually equivalent to 
assuming that the total spin of the molecule 5 = (Sx)"^ + {Sz)'^ remains equal 
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Fig. 3.9 : Perpendicular field dependence of the Mn resonance frequencies, 
calculated according to Eqs. and 11^. 1 y\) assuming 6b = 90° exactly. For 

comparison, we also plot the Larmor frequency of^H. 

to 10 at all fields. By diagonalizing the Hamiltonian 1)3.14(1 one finds that this is 
the case only if the external field is not exactly perpendicular to z plane but has 
(e.g. due to a small misalignment with respect to the crystal axes) a longitudinal 
component B^. This yields a bias ^ = 2gnBSzBz which forces the eigenstates of 
(|3.14|) to localize at one side or the other of the anisotropy barrier, whereas with 
a perfect perpendicular field {6b = 90°) the eigenstates of (|3.14|) are delocalized 
over both sides of the barrier and yield (Sz) = 0, thus S = (Sx)- As a rule of 
thumb, one obtains 5" = 10 when ^ > A, and (Sz) = 0, S ^ (Sx) otherwise. In 
fact, in zero external field = 0) all the eigenstates of (|3.14|) have 5 = 0, which 
would lead to i3hyp(0) = 0, contrary to the experimental findings. The solution to 
this paradox is in the fact that the Hamiltonian (|3.14|) refers to a single, isolated 
molecule. For a realistic description, we must take into account the fact that the 
clusters are embedded in a crystalline structure and interact with each other by 
intercluster dipolar coupling, producing bias fields that act on the electron spin 
levels. 

3.4 Dipolar fields 

In a real sample consisting of a lattice of identical clusters, each molecule is also 
subject to the dipolar field B^ip created by its neighbors, which adds a term 
—gjjiBS ■ -Bdip in the Hamiltonian. A numerical calculation of the statistical dis- 
tribution of such dipolar fields has been performed by I. S. Tupitsyn (Kurchatov 
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Fig. 3.10 : Calculated distribution of the longitudinal component -B^jp of the 
dipolar fields in a ZFC ( circles ) and FC ( squares ) sample in zero external field. 
The solid line is a gaussian fit with width 2(7dip — 0.042 T, whereas the dotted 
line is just a guide for the eye (courtesy of I. S. Tupitsyn). 



Institute Moscow, Russia) by taking a "real" '^^^ 40 x 40 x 40 lattice of Mni2-ac 
molecules, and assuming a certain value of the total magnetization. Fig. 13.101 
shows the distribution of the longitudinal component S^-p of the dipolar fields 
in the case of a fully polarized, field-cooled (FC) sample and a demagnetized, 
zero-field cooled (ZFC) sample. Whereas the FC distribution is asymmetric and 
obviously depends on the size and shape of the crystal, the ZFC distribution can 
be well approximated by a gaussian with width 2(Tdip — 0.042 T. In practice, this 
means that the vast majority of clusters is subject to a longitudinal bias ^ ~ 0.1 
K, several orders of magnitude greater than the tunneling splitting; diagonalizing 
the Hamiltonian 1)3. 14(1 including Ba\x> yields indeed S = 10, as argued on basis 
of the observed hyperfine fields. Nevertheless, there remains the possibility to 
reduce the total spin to less than 10 by applying a perpendicular field strong 
enough to produce a tunneling splitting larger than the dipolar bias. The exact 
way in which A becomes greater than ^ is a delicate problem that must be solved 
by self-consistent numerical iteration, because as soon as some molecules have a 
small enough bias to cause the spin to reduce from 10 to (S^), also the dipolar 
field they produce on their neighbors will suddenly be reduced, causing a drop in 



'^•'That is, not treating the cluster spin as a point dipole but considering the single atomic 
moments. 
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Fig. 3.11 : (a) Perpendicular field dependence of the full width 2c7dip of the 
longitudinal dipolar field distribution (circles), compared to that of the tunneling 
splitting A (solid line), (h) Calculated NMR spectra for the Mn^^'^ site as a 
function of perpendicular field. Solid line: calculated from ii^.l/A taking 5 = 10 
and independent of B± . Circles: calculated from i3.21\) . with S{B±) as obtained 
by averaging over the dipolar fields distribution (courtesy of I. S. Tupitsyn). 



the width of the distribution of dipolar biases. Fig. l^.llT a) shows the full width 
2crdip of the distribution of longitudinal dipolar fields, calculated self-consistently 
on the real Mni2-ac lattice and expressed in millitesla^^^ . 

Since we are interested in the consequences of such a phenomenon for the 
^^Mn NMR spectrum, we first calculated numerically the values of (Sz) and (Sx) 
as function of by averaging over the self-consistently calculated distribution 
of dipolar biases. We then recalculated the total hyperfine field in Mn^^^ sites by 



'^Notice that we may convert from field to energy units by multiplying by the factor gfisS, 
but S{B±) is not a constant! The vertical scales in Fig. 13. lir a) are such that the field and 
energy units match in the low-field region {B± < 5.5 T) where S = 10, but at higher B± only 
the field units scale (right-hand side) should be used for 2crdip. 
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substituting 




Bhyp{Bx) = -Bhyp(O) 



(3.21) 



10 



into Eq. (|3.18() . The resulting Larmor frequency is plotted in Fig. I3.1ir b) and 
compared with the frequency that would be obtained assuming a constant value 
of -Shyp, i.e. a constant total spin S = 10. It clearly appears that the deviation 
from the /S = 10 spectrum occurs when A approaches the width of the dipolar 
field distribution. An actual observation of such a peculiar spectrum would imply 
that the total electron spin of the cluster is occupying a ground state consisting 
of a superposition of macroscopically distinct states for a time longer than the 
intrinsic NMR timescale, and that such a superposition of states is robust enough 
to survive the interaction with the environment and the measurement performed 
via the nuclear spins. To be completely correct, one should include the hyperfine 
interaction li-Ai-Si in the Hamiltonian ()3.14|) and solve the complete problem, 
but the size of the matrices would become intractable. We notice anyway that the 
hyperfine interactions are of the same order of magnitude as the dipolar fields, 
thus we can empirically assume that their effect on the localization of the electron 
spins should not be much larger than what we calculated so far by considering 
only Bdip. 

3.5 Nuclear spin dynamics 

3.5.1 Nuclear spin-lattice relaxation (NSLR) 

The general theory of nuclear relaxation, which expresses the spin-lattice relax- 
ation rate as a function of the time correlation of a perturbing field [HJI, has been 
successfully applied to interpret the temperature dependence of the NSLR rate 
in the thermally activated regime of molecular magnets, first with proton NMR 
in Fes p2ll93llM] and more recently with ^^Mn in Mni2-ac (89119511^6] . We shall 
briefly review it here, mainly to point out why it breaks down upon approach- 
ing the pure quantum regime for the electron spin fluctuations. The implicit 
assumption that has always been made jn21in31in4 , 89 , 95 , 96 is that the nuclear 
spins have a certain Zeeman splitting, thus a certain Larmor frequency ujn, given 
by a static magnetic field, and that the transitions between the Zeeman levels of 
the static Hamiltonian can be computed within perturbation theory, assuming 
the presence of fluctuations in the hyperfine field which are small compared to 
the static field. In most cases, as in the experiments on Mni2-ac |921l93j (with 
the exception of 94 ), the static field is an externally applied field, much larger 
than the hyperfine fields. For ^^Mn in Mni2-ac however, the reverse is true and 
one may carry out an NMR experiment exploiting the local hyperfine field alone. 
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which is made "quasi-static" by the high anisotropy barrier for the reversal of the 
electron spins. Thus, in the absence of an external field, the static field is just the 
value of -Bhyp obtained when the electron spin is in its ground state. Fluctuations 
of the local field may then arise when the cluster spin is thermally excited to a 
higher level. Calling b{t) the fluctuating component of -Bhypj the NSLR rate W 
can still be computed within the framework of perturbation theory if 6 <C ^hyp- 
W is proportional to the spectral density at the Larmor frequency of the perpen- 
dicular component b±{t) of the fluctuating field, i.e. the Fourier transform of the 
time correlation function, {b_L{t)b±{0)): 

W = ^ J {b±{t)b±{0)) e^p{iujNt)dt. (3.22) 

We may further assume that the correlation function of b±(t) is related to the 
phonon-induced fluctuation rate t~^^^ of the cluster electron spin as: 

{b±{t)b^m = {AbD exp(-f/Ts_ph). (3.23) 

Here (A6^) is the square of the average change in the perpendicular component 
of the hyperflne fleld when the electron spin state is excited from m = ±10 to 
m = ±9. The spin-phonon transition rate t~^^^^ depends on the lifetime, rg, 

and the activation energy, AE = Eg — Eiq, of the m = ±9 states as Tglpjj = 
Tg"^ exp(— Ai^/Zc^T), where: 



1 _ Cs_phAS3 



rg 1 - eM-^E/kBT) ' 



(3.24) 



Tg incorporates therefore the spin-phonon coupling constant Cs_ph [1^11111 ■ For 
the purpose of the present discussion we may neglect the higher lying levels, 
\m\ > 8, which is a reasonable approximation below T ~ 3 K since AE > 10 K. 
Inserting into Eq. (|3.22|) we find: 

which shows clearly the expected exponential temperature dependence of the 
NSLR rate, since rg ~ const, for ksT <C AE. 

Eq. (|3.25|1 has the same form as the expression used by Goto et al. > except 
that the latter authors insert the T-dependence into an effective fluctuating fleld, 
expressed as: 

6cfr = A6^ — , 
no 

1 a_„>,AS3 
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where rio is the Ufetime of the m = ±10 states, which is exponentially T- 
dependent. Goto et al. also pointed out that, in order to obtain a transverse 
component of the fluctuating hyperfine field upon excitation of the electron spin 
from m = ±10 to m = ±9, the hyperfine coupling tensor A (|3.16j) must con- 
tain non-vanishing off-diagonal terms. Then we may write b± = —SszAxz/ilN^), 
where 6sz is the change in the z component of the electron spin in each single ion. 
As mentioned in ^3.3.11 only ^^Mn nuclei in Mn^+ have A^z / [HZlinS, while 
A in Mn^"*" sites is a diagonal tensor. We shall come back to the implications of 
this remark in ^4.2.11 and ^4.2.41 

As will become clear from the experiments presented in ^4.21 although Eq. 
(|3.25|) properly fits the data in the thermally activated regime, it fails completely 
in the quantum regime (T < 0.8 K), where the NSLR is found to become nearly 
independent of temperature. This phenomenon is observed in both the proton- 
NMR experiments in Feg [S^ and in our ^^Mn-NMR results, and clearly calls for 
a different theoretical approach. As will be further argued in i i4.21 the NSLR in 
the quantum regime can be ascribed to the hyperfine field fluctuations arising 
from incoherent quantum tunneling of the electron spin between the | ± 10) and 
I — 10) states. In this context, the crucial point is that the complete hyperfine 
field acting on the ^^Mn nuclei in a certain molecule will be suddenly inverted 
when the electron spin of that molecule tunnels, and will stay like that until the 
next tunneling event occurs. This means that a formalism such as that leading to 
Eq. (|3.25|) is unsuitable for describing the ^^Mn relaxation produced by tunneling 
events, because the effect of tunneling cannot be treated as a perturbation on a 
static Hamiltonian. As discussed in ^3.2.31 the inversion of the hyperfine field 
takes place in a time rtr ~ 10~^^ s, i.e. instantaneous compared to ujj^^ ~ 10~^ 
s, whereas the interval tt between tunneling events is, as we shall show below, 
very long compared to the NMR window. The description of the nuclear spin 
dynamics in the presence of incoherent tunneling requires therefore a completely 
different approach, for which the basis has been laid by the Prokof'ev-Stamp 
theory Jl], which we shall review in ^3.61 

3.5.2 Transverse spin-spin relaxation (TSSR) 

Another essential aspect of the nuclear spin dynamics is the transverse spin-spin 
relaxation (TSSR). Under the same hypotheses as for the NSLR, Goto et al. [HS] 
have derived an expression for the TSSR rate for the ^^Mn, that is applicable 
in the thermally activated regime of Mni2-ac: 



--1 



1 (7jvA6^r9)2 



(3.27) 



no 1 + {-iN^hzTg) 



2 ' 



where A^^ is the change in longitudinal hyperfine field whenever the electron 
spin makes a transition between m = ±10 and m = ±9. In this approach the 
temperature dependence is incorporated in the lifetime tiq of the m = ±10 states. 
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Based on the analysis of the field-dependence of T2 ^ , they also conclude that the 
regime (7ArA6^r9) » 1 holds, so that (|3.27j) becomes: 



{AE/kBT)-l 



(3.28) 



In this model the decay of the transverse magnetization is basically due to 
a phase shift that occurs locally and randomly as a consequence of the sudden 
changes in the longitudinal component of the hyperfine hyperfine field (and thus 
in the Larmor frequency), for the nuclei in a thermally-excited cluster. In fact, 
this mechanism describes the dephasing of nuclear precession without considering 
any interaction between the nuclear spins. 

The nuclear dipole-dipole interaction would add an extra contribution, whose 
magnitude can be calculated by taking into account the "flip-flop" term "Hf-f in 
the dipolar Hamiltonian for two like spins /^■'^ and /'•'^^ placed at a distance rjk 
and an angle 9jk from each other |91j : 

Uf^f = -^(1 - 3cos2 e,k){lfl^^^ + ). (3.29) 



Using the Van Vleck formula one then obtains the contribution of the flip-flop 
term to the secc 
its square root: 



term to the second moment M2 of the resonance line, and the TSSR rate T2 ^ as 



k Jfc 



(3.30) 



Applying this formula to, for instance, the nuclei in Mn^^^ sites, gives Tg"^ = 1790 
s~^ if the Mn*-^-* nuclei in the same cluster are included, whereas restricting the 
summation to nuclei in neighboring clusters leads to = 111 s~^. The latter 
value gives an order of magnitude of the rate of intercluster nuclear spin diffusion 
which, as we shall discuss below, is an essential addition to the Prokof'ev-Stamp 
theory needed to understand our experiments. 



3.6 Prokof'ev-Stamp theory 

The Prokof 'ev- Stamp (PS) theory of the spin bath [97| ll4 ( n3] presents a major 
step forward in understanding the quantum dynamics of a giant spin (in our 
case the cluster spin) coupled to a bath of environmental spins (here the nuclear 
spins). The starting point for its development was to recognize that, in several 
physical systems, the "oscillator bath" theory [HH] cannot be applied, because the 
environmental degrees of freedom cannot be described as a set of non-interacting 
oscillators, and the couplings to the central spin are not weak. Initially [97| [T H . 
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the PS theory was developed for a single central spin. The short-time \/t law 
for the quantum relaxation of an ensemble of spins (namely SMMs) has been 
obtained by including the effect of intercluster dipolar coupling in the form of 
an effective (bias) field [221 • This section contains a review the essential aspects 
of the PS theory, with emphasis on the issues that will be necessary to interpret 
the experimental results presented chapter IV. Fig. 13.121 can be used as a visual 
guide through the ingredients necessary to arrive at the complete description of 
the problem. 



3.6.1 Central spin + spin bath 

The first step for describing the coupled system of "central spin + spin bath" is 
to consider a single giant spin subject to crystal fields and, eventually, a magnetic 
field, and to truncate its total Hamiltonian to the ground doublet, neglecting all 
the physics at energies > (see p.2.Hj) : in this way the central spin can be 
treated as an effective spin 1/2 described by Pauli matrices f , where | ft) and | ij-) 
are the eigenstates of f^. In particular, for an isolated spin having a tunneling 
splitting A in the ground doublet, the truncated Hamiltonian is: 

A C 

Wtrunc = -^T^ COs{ttS) - -f^, (3.31) 

where cos(7rS') expresses the fact that in half- integer spins there is a destructive 
interference of tunneling paths |lUL)lllL)lj . and ^ is a longitudinal bias^''^ If such a 
bias is produced by an external static field, it is labelled = —'^gfJ-sS • B. As is 
well known from textbook quantum mechanics in the absence of any other 
couplings the probability i-*-ft-fi-(t) for S to remain in the state | ft) at time t when 
prepared in | ft) at t = is: 

Pnit) = 1 -^-^sm\Et), (3.32) 



where E = ^{i/2Y + (A/2)2. 

The central spin is further assumed to be coupled to N environmental (e.g. 
nuclear) spins ct^. A; = 1 . . . A^, assumed for simplicity to have spin 1/2, described 
by vector Pauli matrices at- In the notation of Prokof'ev and Stamp, the single 
hyperfine couplings have strength hjJk = 7fc • where 7^ is the hyperfine field 
expressed in energy units^^^ and they are centered around an average value hujQ. 
The energy levels of the central spin are now represented by two hyperfine-split 
manifolds, each containing 2-^ states distributed over a Gaussian with half-width 

'^'in the whole Prokof 'ev-Stamp hterature, A and ^ denote the matrix elements in l]:-!..'-!!^ . 
whereas throughout this thesis we always use those symbols for "real" energy differences. There- 
fore it always holds A(this thesis) = 2A(PS) and ^(this thesis) = 2^(PS), and all equations in 
this section are consistently adapted to the present notation. 

'^'With the notation used in the rest of this thesis, 7^ = — /i7Ar-Bhyp- 
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Fig. 3.12 : "Top-down" sketch of the model used by Prokof'ev and Stamp to 
describe the coupled system "central spin + spin bath". The numbers in brackets 
indicate the equations that describe each further step, (a) Giant spin with crystal 
field anisotropy and external bias, (b) Truncation to the ground doublet, with 
static bias S^b- (c) Coupling with nuclear spins: this adds an extra bias which 
depends on the polarization state AN = — of the nuclei (cf. Fig. \!^.14^ . 
(d) Nuclear spin diffusion: each polarization group is broadened by hT^, so that 
groups on either side of the barrier may overlap and create a tunneling window. 
Calling AN^^^ and AN^'^^ the nuclear polarizations before and after the flip, the 
tunneling probability is governed by the effective tunneling splitting Am, with 
2M = AiV(2) _ AiV(i). 
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Fig. 3.13 : Scheme of the relative orientations of the hyperfine fields before 
('y'k^) '^^^ after ('y^^ ) the electron spin flip, and the consequent sum (hijj\) and 
difference (hujj: ) terms in the hyperfine energies, as defined in Eq. I13.3!:J\) . The 
angle (5^ is typically very small since iJl » uj-^ . 



^ y/NcoQ. Degenerate states within the manifold can be labelled by their "polar- 
ization group", i.e. according to the total nuclear polarization AA^ = — 
(see Fig. 1X11) . 

By calling 7^^^ and 7^.^^ the hyperfine fields before and after the tunneling of 
the central spin, respectively, we can define the sum and difference terms: 

^ II r ^(1) ^(2) 

^M = 7i'^+7f (3.33) 



where 1^ and rhk are unit vectors (cf. Fig. 13.13(1 . fiwj: represents the component 
of the hyperfine coupling energy that does not change during a flip of the central 
spin, whereas fiuJ^ is the variable part. For instance, in the case of ^^Mn nuclei 
in Mni2-ac, one would find uj^ ~ since the hyperfine fields acting on before 
and after the flip are just antiparallel and equal in magnitude^^^ . These terms 
add a static (Zeeman) contribution to the Hamiltonian of the coupled system: 



^static = ^ ( X] ^i^'^ ■ "^fc + XI ^fc -^k] ■ (3.34) 



2 , 

\ A: k 



The first term is a longitudinal coupling that yields an internal field bias acting 



'^In fact, as long as we strictly consider only one electron spin, 7^ — —7^ always. In 
most realistic cases, the nuclei would be subject to the field of different electron spins (or to an 
externally applied field), so that ui^ 7^ is often true, although usually w|! S> uj^t- 
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on Tz that depends on the nuclear polarization, yielding a hyperfine bias'^^*'^ = 
hujQlS.N (see Fig. K-{.14j) . Formally it plays the same role as E,b discussed above, 
with the crucial difference that rnay vary in time, as we shall see below. 
The second term causes the so-called "orthogonality blocking", since when 7^"^^ 7^ 
—7^ , the basis for the eigenstates of Ok (defined taking the quantization axis 
parallel to the instantaneous local field direction) after the flip of 5 is not exactly 
orthogonal to the basis before the flip. This means that after the flip cr^ will start 
to precess around a different axis, which quantum mechanically is equivalent to a 
transition between different nuclear Zeeman levels. The number of nuclear spins 
that would coflip by this mechanism is given by the parameter k, defined as (cf. 
Fig.Einifor /3fc): 



W cos (3k 



■ 2 S fc /^fc 



-(2) 

Clearly, k depends only on the direction of the hyperfine fields, and not on the 
timescale of the central spin flip. 

Furthermore, the presence of environmental spins adds extra Berry phase 
terms to the cos(7rS') discussed above, leading to the so-called "topological deco- 
herence". The detailed derivation of the total phase is not essential here, but we 
mention that it contains terms like = {11/2)00^/^0, which depend on the ratio 
between the nuclear Larmor frequencies and the bounce frequency of the central 
spin (cf. ^3.2.3(1 . aj. is the probability for to coflip with S (see also |lU2j '). but 
now for a different reason than the orthogonality blocking. What matters here is 
only the timescale of the flip compared to the period of Larmor precession, i.e. 
whether or not the nuclei can adiabatically follow the rotation of S. A parameter 

A = l/2j]a2 (3_3g) 

k 

is introduced, that represents that number of nuclear spins that coflip with S 
because of the topological decoherence mechanism alone. 

'^"'Let us recall that, throughout this thesis, we always express the hyperfine field -Bhyp as 
the field created by the electron spin and acting on the nuclear spin, which produces a nuclear 
Zeeman splitting tuu = ft7_Bhyp- Obviously, this energy is exactly the same as the bias 
produced by the nuclear spin on the electron spin levels, which some readers may like to express 
as — — 2(?/ifl_B^yp • S, where B'^^^^ is the field produced by the nuclear spin and acting on 
the electron spin. Therefore, there is no inconsistency in writing, for instance, the Hamiltonian 
13.3411 for the central spin in terms of the hyperfine field acting on the nuclei! Moreover, it 
should be clear from Fig. 13. 141 that the total bias on the electron spin levels can be obtained 
by summing up the Zeeman splittings of the nuclei, provided we account (by means of AA*') 
for their polarization state. The only subtlety is that, in general, the nuclear spins cannot be 
treated as a classical source of external static field acting on the electron spin, since there is the 
possibility that some nuclei may coflip with 5, thereby changing AA'^. 
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Fig. 3.14 : Sketch of the hyper fine- split manifold obtained by coupling the 
central spin to N nuclear spins, with N = 1,2, . . ., assuming that all spins have 
I = 1/2 and Zeeman splitting Hujo. Notice how the bias S^n = HloqAN on the 
electron spin levels depends on the nuclear polarization AN. This picture is valid 
only when all the nuclei have the same loq; in a realistic case, the spread in 
hyperfine coupling makes the energy level scheme much more intricate. 



An additional aspect that must be included for a realistic description is the 
"degeneracy blocking", which consists in allowing the hyperfine couplings to have 
a spread h6uJk around the central value hujQ. This spread can have different ori- 
gins (quadrupolar splitting, nuclear dipole-dipole interaction, transfer hyperfine 
interactions, Nakamura-Suhl interactions |in3j . etc.) and could potentially block 
the tunneling dynamics of the electron spin, since it breaks the degeneracy of 
the hyperfine-split electron spin states within the same polarization group. Be- 
cause of the spread of nuclear frequencies, each polarization group is no longer 
represented by a sharp line, but rather by a Gaussian peak of width 

= ^5L0k. (3.37) 

In most cases, the gap between unperturbed polarization groups is much smaller 
than r^, so that different polarization groups are in fact largely overlapping. 
Because of the mutual dipolar coupling between nuclear spins: 

Wdip({a,}) = Yl ^kk'^k^kn (3.38) 

k^k' 

there is an intrinsic broadening Suj^^^^ ~ Tg"^ of the nuclear Zeeman levels. The 
crucial point is to recognize that this broadening is dynamic, due to nuclear spin 
diffusion. In practice, the hyperfine bias can fluctuate and cover the whole range of 
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each polarization group, i.e. an amplitude hT^, on a timescale ~ T2. The effect of 
spin diffusion is therefore to eliminate the degeneracy blocking. Notice that this 
sort of fluctuations arises from what one would call "intracluster" spin diffusion, 
since the model is written for one single central spin coupled to his nuclear spins. 
The transition between different polarization groups would therefore require the 
exchange of energy with an external reservoir, on a timescale Ti; in the PS theory, 
this mechanism is assumed to take very long times. On the other hand, by taking 
many identical "central spin + spin bath" units, one may allow for "intercluster" 
spin diffusion as well, i.e. between equivalent nuclei belonging to different clusters. 
This possibility is not explicitly discussed in the early PS papers, but will play 
an essential role in the interpretation of our experimental results. 



3.6.2 Tunneling rate 

The strategy adopted by Prokof 'ev and Stamp to derive the incoherent tunneling 
rate of the central spin coupled to a spin bath, is to estimate it from the evolution 
of the probability P-f^-^it) for S to remain in the same state | f^) until time t. Its 
trivial expression for an isolated spin is Eq. (|3.32j) : including the spin bath one 
has to average over the different values Puit) that P^-^{t) may take, depending 
on the change 2M of the polarization state AA^ = A^^ — A^^ of the environmental 
spins caused by the flip of S. Ppiit) is itself obtained by averaging, over the 
"topological" and "orthogonality blocked" processes, the quantity P^^\t), defined 
as: 



Pj^\t) = l-^-^^sin'{EMt) (3.39) 
Ell = (e/2)' + (AAf /2)2 (3.40) 

(X _ \')M/2 



where A' = A in the case of pure topological decoherence and A' = for pure 
orthogonality blocking, and for we take the total bias £, = S,N + which may 
include the effect of an externally applied field. Roughly speaking, this means that 
the effect of the spin bath is incorporated in a renormalization of the tunneling 
splitting, assigning a vanishingly small value of Am to those tunneling transitions 
that would require too many nuclear coflips, as compared to the "natural values" 
of K and A. The really essential point is that ^ in Eq. 1)3. 4U() is no longer a static 
bias, as soon as nuclear spin diffusion is allowed. Prom (|3.38|) one finds that 
contains a time-dependent part which fluctuates on a timescale Tg"^ over a 
range ~ KT^, which is the change in longitudinal hyperfine bias when pairwise 
I Ti) ~^ I it) nuclear flip-flops occur. With the assumption that T2 is much 
shorter than the typical tunneling interval tt, the average of P^\t) over the 
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fluctuating bias can be written as: 



(<'(*.0>«-5 = 5<!-"'»' (3.42) 

The parameter represents therefore the rate for S to tunnel accompanied 
by the coflip of M nuclear spins. In the presence of an external static bias 
r^l starts to be suppressed as soon as £,b is larger than the typical hyperfine 
fleld spread. Combining this with the fact that the largest tunneling splitting is 
obtained for M = 0, one finds that the leading term in the tunneling rate takes 
the form: 

^o_,-|?s|/eo. (3.44) 



2^/i2r 

where ^0 is the width of the "tunneling window" . When some nuclear spins may 
coflip with S by topological decoherence (A) or orthogonality blocking (k), then 
^0 is given by the energy that the central spin can exchange with the nuclear spin 
bath, ^0 ~ (-^ + k)^o- Otherwise, in the absence of nuclear coflips the tunneling 
window takes the width ^0 ~ due to nuclear spin diffusion only. From this 
discussion it is clear that the major role in creating a useful tunneling window 
for the central spin is played by those nuclei that more easily coflip with S. We 
anticipate that this is an experimentally testable prediction ( H4.2.3() . 

The global relaxation, i.e. the precise form of P-ff^{t), is obtained by combining 
all the results discussed above. This leads to very interesting phenomena to be 
observed in the quantum relaxation, but goes beyond the scope of our research. 
We shall mainly make use of the tunneling rates (|3.43)) . (|3.44j) and the concepts 
of topological decoherence and orthogonality blocking to evaluate the numbers 
of coflipping nuclei A and k. 

3.6.3 Comparison with the Landau-Zener formalism 

It is rather instructive to compare the Prokof 'ev-Stamp theory with the Landau- 
Zener formula for the calculation of an incoherent tunneling rate. In the 
presence of nuclear spin diffusion, assuming for simplicity the external bias = 
0, the bias ^ fluctuates over a range hT^ on a timescale T2, thereby forcing the 
energy levels of the central spin S to cross back and forth through the tunneling 
resonance: a sketch of the situation is given in Fig. 13.151 We could therefore 
think of applying Eq. 1)3. simplified by noting that is a very small number 
and expanding the exponential: 

7rA2 /dA"^ , , 

^--^ dJ • ^'-''^ 
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t 



Fig. 3.15 : Time evolution of the bias ^ on the central spin levels as a con- 
sequence of nuclear spin diffusion: the crossing through the tunneling resonance 
happens with slope d^/dt ~ hT^/T2. 

Clearly, in our case d^/dt ~ ^r^/Ta, yielding Plz ~ (A^ /h)[{T2/{hT^)]. To 
obtain the tunneling rate, we multiply the single tunneling probability by the 
number of crossings per unit time, i.e. Tg"^, and obtain 



'LZ 



(3.46) 



which has indeed the same form as Eq. (|3.43|) . Deriving it via the LZ formalism 
clarifies why the tunneling rate does not contain the timescale T2 of the bias 
fluctuations. On the other hand, by this we have completely lost any information 
about the effect of the central spin dynamics on its nuclei, which in Eq. ()3.43|) is 



incorporated in the renormalized tunneling splitting Am, cf. (|3.41|) . An essential 
feature of the PS theory is indeed that the fluctuating bias is not added artificially 
as an "environmental noise" , uncorrelated with the dynamics of S, but is built in 
the model in a self-consistent way, including the back-action of S on the nuclear 
spins [TUl] . 



3.6.4 Application to SMMs 

When applying the PS theory to the combined electron-nuclear spin dynamics in 
SMMs, it is useful to distinguish between nuclei coupled directly to one single ion, 
like ^^Mn in Mni2-ac, and nuclei located at the ligand molecules, like (and a 
small fraction of ^'^C). In the latter case, the strength of the coupling is of course 
much lower than for ^^Mn, since it is due mainly to dipolar fields, but the total 
field is the sum of the dipolar fields from different molecules. If one molecule flips, 
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it is only one component of the total field at the H site that flips, which may 
yield a situation where /3k, and therefore the chance of nuclear coflip, is rather 
large (cf. Fig. 13. 8() . This means that the protons can be the main contributors 
to the width of the tunneling window , despite the relatively weak coupling. 

On the contrary, /3fc ~ in ^^Mn, since the whole hyperfine field is inverted 
by 180° upon tunneling. This yields k ~ and, because of the short tunneling 
traversal time (cf. ^3.2.3() . A ~ as well. Clearly, we can expect the change of 
polarization group for ^^Mn nuclei to be very difficult, unless nuclei in different 
clusters are linked with each other. This possibility is not explicitly discussed in 
the PS theory, but we shall see in the next chapter that it is both experimentally 
observed (i i4.2.Hl . and necessary to explain the observed ^^Mn NSLR rate in 
Mni2-ac. 



IV 
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The nuclear spin dynamics in a single-molecule magnet is one of the most crucial 
aspects influencing its quantum behavior, nonetheless the experimental knowl- 
edge was so far almost totally lacking, in particular in the low-T quantum regime 
for the cluster spin fluctuations. We dedicate therefore this chapter to a de- 
tailed survey of our results on the ^^Mn NMR in Mni2-ac, starting with the 
temperature dependence of the nuclear spin-lattice relaxation (NSLR) and trans- 
verse spin-spin relaxation rates (TSSR) in zero external field. The data indicate 
that quantum tunneling of the electron spin is the mechanism responsible for the 
observed features in the nuclear spin dynamics. The dependencies on small longi- 
tudinal fields, magnetization state of the electron spins, isotope substitution, and 
crystallographic Mn site, yield precious additional information. We then discuss 
the influence on the NSLR of an external perpendicular field, which can be used 
to tune the tunneling splitting. Finally, an essential and original information 
is provided by the study of the nuclear spin temperature. The results are very 
surprising and call for an extension of the existing theory of quantum tunneling 
in magnetic molecules. The chapter is concluded by a suggestion for an approach 
that should unify all the experimental evidences into a new and more complete 
theoretical description. 



4.1 Measurements and data analysis 

About 100 mg of Mni2-ac crystallites were mixed with Stycast 1266 epoxy and 
introduced into a 6 mm capsule, then inserted in the room-temperature bore 
of a 9.4 T magnet and allowed to orient for 24 hours. A two-turns NMR coil was 
wound around the capsule, the whole assembly was inserted inside the mixing 
chamber of the dilution refrigerator f ^2.1() and connected to the NMR electronics 
as described in H2.2.2I 

The ^^Mn nuclear precession was detected by the spin-echo technique. A 
typical pulse sequence includes a first 7r/2-pulse with duration = 12 ^s, 
a waiting interval of 45 ^s, and a 24 ^s vr-pulse for refocusing. At the lowest 
temperatures provided by the dilution refrigerator, such a pulse sequence was 
found to cause an increase in temperature of about 5 mK in the ^He/^He mixture. 
With a waiting time of 600 s between subsequent pulse trains we could easily keep 
the operating temperature around 15 — 20 mK (cf. Fig. I2.5|l . On the other hand. 
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Fig. 4-i ■ example of "real time" echo signals recorded during an inversion 
recovery, i.e. measuring the echo intensity at increasing delays after an inversion 
pulse. In particular, these are single-shot (no averaging) raw data taken at B = 
and T = 20 mK in the Mn^^^ site. The (normalized) integral of the echoes is 
reported in Fig. \4-'A open squares. 

at such low temperature the signal intensity is so high that we could obtain an 
excellent S/N ratio without need of averaging, so that a typical measurement 
sequence (e.g. an inversion recovery) would last less than 12 hours. Above 100 
mK it proved convenient to take a few averages, but then the heating due to the 
r/-pulses becomes negligible, and the waiting time can be reduced to ~ 100 s. 

The NSLR was studied by measuring the recovery of the nuclear magnetization 
after an inversion pulse (cf. Fig. 12., '-{|) . We preferred this technique to the more 
widely-used saturation recovery [HHllHZlinSI because it avoids the heating effects 
of the saturation pulse train, but we checked at intermediate temperatures that 
the two techniques indeed lead to the same value of NSLR rate. An example of 
echo signals obtained as a function of the waiting time after the inversion pulse 
is shown in Fig. 14.11 By integrating the echo intensity we obtained recovery 
curves such as those shown in Fig. 14.21 For the ease of comparison between 
different curves, we renormalize the vertical scale such that M(0)/Af(oo) = — 1 
and M{t > Ti)/M(oo) = 1, even though usually |M(0)| < |M(oo)|, as could be 
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deduced from Fig. 14. Ih ^^. Since tlie ^^Mn nuclei have spin I = 5/2, the recovery 
of the nuclear relaxation for the central line in the quadrupolar split manifold is 
described by jlL)5j : 



M{t) 
M(oo) 



lOQ -30m I ]^^12Wt , ^ -2Wt 

63 45 35 



(4.1) 



^00 -(30V|/t)° I }^ -(12Wt)°' I _^ -(2m)" 

63 45 35 



(4.2) 



where is the nuclear spin-lattice relaxation rate^^^ . All data in zero field and 
moderate parallel fields could be accurately fitted by Eq. 1)4. ip . which indicates 
that the relaxation mechanism can be characterized by a single rate W . Con- 
versely, a proper analysis of the inversion recoveries in strong perpendicular fields 
requires the introduction of a stretching exponent a: 

M{t) _ 
M(oo) 

This can be attributed to a small distribution of orientations of the crystallites in 
the sample, which causes the NMR signal in high fields to consist of a mixture of 
lines from different crystallites, each of them characterized by a different NSLR 
rate. A comparison between an inversion recovery at zero field, fitted by Eq. 
1)4. 1() . and a similar recovery in strong perpendicular field fitted by Eq. 1)4. 2|1 . is 
shown in Fig. 14.21 

The TSSR rate was studied by measuring the decay of echo intensity upon 
increasing the waiting time r between the 7r/2- and the vr-pulses. The decay of 
transverse magnetization M_|_(r) can be fitted by a single exponential 

^'^-(^-) - exp i-p^ (4.3) 



except at the lowest temperatures (T < 0.2 K), where also a gaussian component 
Ti^Q needs to be included: 

^'^^(-)-e.pf-^V^pf-S^) (4.4) 



M±(0) V T^) \ -il^aJ 

As regards the experiments to determine the nuclear spin temperature, the 
measurements were performed by monitoring the echo intensity at regular in- 
tervals while changing the temperature Tbath of the ^He/^He bath in which the 
sample is immersed. Recalling f N2.2.1() that nuclear magnetization is related to 
the nuclear spin temperature Ta_uc\ by the Curie law: 

M(T„„) = A'^^|^. (4.5) 

•J'^^B-inucl 



'^'This is just an artifact that occurs when NMR line is much broader than the spectrum of 
the inversion pulse, and does not mean that the length of the vr-pulse is incorrect. 

'■^'in the simple case of a spin 1/2, where the NSLR is described by a single exponential, W is 
related to the relaxation time T\ by 2W = T^^ . 
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Time after inversion (s) 

Fig. 4-2 : Examples of nuclear inversion recoveries in zero field (squares) and 
strong perpendicular field (circles). The NSLR rate W is extracted by fitting the 
data with Eq. or Eq. i4-S^ , respectively. 

and assuming that Tbath = T'nuci at a certain temperature Tq (e.g. 0.8 K), 
we can define a calibration factor K sucfi tliat M(To) = i^/Tnuci(^o) and use 
that definition to derive the time evolution of the nuclear spin temperature as 
^nucK^) = K/M{t) while the bath temperature is changed. 

4.2 Nuclear relaxation in moderate parallel fields 

In this section we present and discuss the experimental results on the ^^Mn spin- 
lattice and spin-spin relaxation rates in Mni2-ac, starting with the temperature 
dependencies in zero field, and then discussing the effect of small longitudinal 
fields and isotope substitution. All data refer to Mn^^^ sites, except for ^4.2.41 
where a comparison with the Mn^^^ site is presented. 

4.2.1 NSLR and TSSR in zero field 

As mentioned in ^3.5.H the NSLR and TSSR in Mni2-ac have already been 
studied in quite some detail in the thermally activated regime by other groups 
[HnilMlinE]) whereas the purpose of our work has been to extend the analysis 
deep into the pure quantum regime |Hn|- We begin by discussing the temperature 
dependencies of the NSLR and TSSR rates in a demagnetized, zero-field-cooled 
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(ZFC) sample, shown in Fig. 14.31 

What clearly appears is a sharp crossover at T ~ 0.8 K between a roughly ex- 
ponential T-dependence and an almost T-independent plateau. We attribute the 
T-independent nuclear relaxation to the effect of ground-state tunneling fluctua- 
tions of the cluster spins, and we shall dedicate most of this chapter to discuss our 
further results supporting this statement. It should be noted that the crossover 
from thermally activated to ground-state tunneling has been also observed by an- 
alyzing the T-dependence of the steps in the magnetization loops |in6|lin7] . The 
important advantage of our NMR measurements is that the nuclear dynamics 
is sensitive to fluctuations of the cluster electron spins without even requiring a 
change in the macroscopic magnetization of the sample. Clearly, no macroscopic 
probe (except perhaps an extremely sensitive magnetic noise detector) would be 
able to detect the presence of tunneling fluctuations in a zero-field cooled sample 
in zero external field, since the total magnetization is zero and remains so. In 
fact, below T ~ 1.5 K the steps in the hysteresis loops of Mni2-ac can be ob- 
served only at relatively high values of external field |1U6U1U7] . which is indeed 
predicted to lead to a less sharp transition between thermally activated and quan- 
tum regimes [TUHl- We fitted the high-T part of the NSLR rate by Eq. (HT^ . 
as shown by the dotted line in Fig. 14.31 the curve depicted is obtained by fixing 
AE = Z)(102 - 9^) + BiilO"^ - 9"^) = UA K which yields {AbD/Tg ~ 1.1 x 10^ 

s~^. The fit, in particular its slope, may be improved by leaving also AE as 
a free parameter, which then yields AE ~ 12.2 K and {Ab]) /rg ~ 2.2 x 10^ T^ 
s~^. Despite the reasonable quality of the fit, Eq. (|3.25|) is based on a model 
which is seriously questionable for two reasons: (i) it considers only hyperfine field 
fluctuation due to "intrawell" excitation, i.e. it does not account for thermally- 
assisted tunneling through excited doublets, and (ii) it does not explain what is 
the source of transverse fluctuating field b±, since the hyperfine coupling tensor 
in Mn'^^ does not contain non-diagonal elements (cf. ^3.5.1|) . We may already 
notice that including "interwell" fluctuations due to thermally assisted tunneling 
could solve the paradox of the transverse hyperfine field component, provided we 
can demonstrate that quantum tunneling of the electron spin is a possible source 
of nuclear relaxation. We shall come back to this issue in ^4.2.41 

The roughly T-independent plateau in the NSLR rate below T ~ 0.8 K is 
characterized by a value of ~ 0.03 s~^ that is surprisingly high, which at first 
sight may appear like an argument against the interpretation in terms of tunneling 
fluctuations of the electron spin. Experimentally it is indeed well known jl09j 
that the relaxation of the magnetization in Mni2-ac in zero fleld may take years 
at low T, which means that the tunneling events are in fact extremely rare. Based 
on this, we are forced to assume that the tunneling events take place only in a 
small minority of the clusters, and that some additional mechanism takes care 
of the relaxation of the nuclei in molecules that do not tunnel. This is a very 
realistic assumption indeed, since it is well known that all samples of Mni2-ac 
contain a fraction of fast-relaxing molecules (FRMs) [HS1IH2], as anticipated in 
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^3.2.41 Moreover, that the tunnehng dynamics has to be ascribed to FRMs is 
supported by the experimental fact that, when measuring the NSLR in zero field 
in a saturated field-cooled sample (FC), we observe that the sample magnetization 
does not change during the experiment. This can be done by applying a small 
longitudinal field Bz ~ 0.5 T and measuring the NMR signal intensity at the 
shifted frequencies u-^ = + jnBz and z^j = z^(0) — ^nBz, corresponding to 
nuclei in clusters whose spin is respectively parallel or antiparallel to the external 
field |lin| l87]. In a fully saturated sample, the NMR signal at z^j is indeed zero, 
and we verified that it remains that way even after several weeks of experiments 
as long as the sample is kept constantly at low temperature. In other words, 
the macroscopic magnetization as seen by the nuclei that contribute to the NMR 
signal does not change in time, implying that the observable NMR signal comes 
from nuclei that belong to "frozen" clusters^^\ therefore the relaxation channel 
must involve the tunneling in FRMs. 

The question remains what is the mechanism that links the tunneling in FRMs 
to the nuclear spins in frozen molecules. One possibility is to ascribe it to the 
fiuctuating dipolar field produced by a tunneling FRM at the nuclear sites of 
neighboring frozen molecules. In that case we may give an estimate of W using 
an expression of the form: 

^ - ^&L:r^V^ « J^-T-^ (4.6) 



4^°^Pl + a;^T^ 4B 



tot 



where 6dip is the perpendicular component of the fluctuating dipolar field pro- 
duced by a tunneling molecule on its neighbors and t^^ is the tunneling rate. 
Even for nearest neighbors ftjip < 3 mT, which leads to the condition W ~ 0.03 
s^^ =^ tt ^ 10^ s^^. This is completely unrealistic, even for the FRMs with two 
flipped Jahn- Teller axes. It is therefore unavoidable to consider the effect of a 
tunneling molecule on the nuclei that belong to the molecule itself, and to look 
for an additional mechanism that links nuclei in FRMs with equivalent nuclei in 
frozen clusters. It is natural to seek the origin of such a mechanism in the inter- 
cluster nuclear spin diffusion, and we shall provide strong experimental evidences 
to support such an interpretation. 

By magnifying the NSLR rate in the quantum regime (inset of Fig. 14.31 we 
notice that the plateau is in fact not perfectly horizontal but seems to have 
two slopes, with most pronounced one for T > 0.3 K. We may attribute this 
slope to the thermally-assisted intrawell fluctuations in the FRMs, i.e. the same 
phenomenon as observed for T > 0.8 in the whole sample. Since the FRMs have 
a lower anisotropy and therefore lower AE, we may expect their crossover from 
quantum to thermally activated regime to take place at a lower temperature. 



We cannot tell for sure whether also nuclei in FRMs participate to the signal, but the chance 
is indeed small since their T2 would be very short, due to the quickly changing hyperfine field 
they experience. 
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Fig. 4-3 : Temperature- dependence of the NSLR (circles) and TSSR (squares) 
rates for a ZFC sample in zero field and v = 231 MHz. The dashed lines are fits 
to Eqs. \'-!.25\) and \H.2^) . with parameters discussed in the text. The inset is an 
enlargement ofW{T) for T < 0.8 K, with lines as guides for the eye. 
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Magnetization measurements |82| have indicated that the crossover takes place 
at T ~ 0.6 K in the FRMs with just one flipped Jahn- Teller axis, but this value 
could be of course lower for the FRMs with two flipped axis and 15 K barrier. 

Fig. ESI also shows the T-dependence of the TSSR rate T2~^(T). We observe 
that above 0.8 K T2^{T) has the same T-dependence as the NSLR rate, and can 
be fitted by Eq. (|.S.28|) with AE = 14.4 K. This indicates that the fluctuations 
arising from thermal activation of the electron spin are responsible for the TSSR 
as well, at least in this temperature regime. 

Below 0.8 K the TSSR also saturates to a nearly T-independent plateau. By 
fitting the decay of transverse magnetization with Eq. (|4.3|) at all temperatures 
we obtain — 100 s~^. This is just a rough estimate since, as mentioned in 
^4.H the transverse relaxation is not well described by a single exponential decay 
below T < 0.2 K. Nevertheless, we find this estimate to be comparable with the 
value Tg"^ = 111 s~^ calculated on basis of the dipole-dipole coupling of nuclei 
in neighboring clusters ( H'A.^.'Z^ . It is of great interest to point out that in a ZFC 
sample, to which the data in Fig. 14.31 refer, only half of the nuclei in neighboring 
molecules should be taken into account as possible partners for flip-flop transi- 
tions. Since half of the cluster spins point up and the other half point down, 
the nuclear spins in equivalent sites should be split in two groups having Larmor 
frequencies +iON' and —a; at, respectively. For the purpose of calculating the TSSR 
rate that results from intercluster nuclear spin diffusion, this is equivalent to hav- 
ing diluted the system by a factor 2. Carrying out the sum 1)3. 30() over half of the 
spins only, reduces the calculated TSSR rate by a factor \/2, yielding T2 = 78 
s~^ for a ZFC sample, whereas the above-mentioned T2~^ = 111 s~^ remains the 
correct estimate for a FC sample. Moreover, the dipole-dipole coupling should 
lead to a roughly gaussian decay of the transverse nuclear magnetization [UJ . 

We have indeed been able to verify both these predictions by measuring the 
TSSR at T = 20 mK in a FC and a ZFC sample, as shown in Fig. The 
decay of the transverse magnetization is indeed best fitted by Eq. 1)4.4(1 . whereby 
the Gaussian component, T^q, is separated from the Lorentzian one, T^j^. From 
the Gaussian component of the decay we can extract directly the effect of the 
nuclear dipole-dipole interaction, whereas the other mechanisms of dephasing 
(e.g. random changes in the local field, cf. ^3.5.2(1 contribute mainly to the 
Lorentzian part. The fit yields r2;^^(FC) = 104 ± 3 s"^ and T^gCZ^FC) = 77 ±3 
s~^, respectively; the ratio r2^(FC)/T2^^(ZFC) = 1.35 of the TSSR rates in 
FC and ZFC samples agrees indeed with the prediction of a -v/2 reduction of 
the intercluster flip-flop rate in a ZFC sample^^^ and constitutes solid evidence 
for the presence of intercluster nuclear spin diffusion. This experimental finding 
points to the need of an extension of the basic assumption of the Prokof'ev- 
Stamp theory jl4j to the case of an assembly of interacting "central spin - spin 



'^In reality it is likely that, even in the ZFC sample, some clusters may have locally the same 
spin orientation, thereby reducing the ratio to something less than y^, as observed. 
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Fig. 4-4 ■ Decay of transverse nuclear magnetization in FC (open circles) and 
ZFC (full circles) sample at T = 20 mK, zero applied field and u = 231 MHz. The 
lines are fits to Eq. yielding the ratio T'^iYC) /T:2^{ZFC) = 1.35 ~ \/2. 

The inset sketches represent pictorially the fact that intercluster spin diffusion 
is possible in a FC sample since all the nuclei have the sam Larmor frequency, 
contrary to the case of a ZFC sample. 

bath" units. It also supports our previous interpretation of the mechanism of 
nuclear spin-lattice relaxation in terms of the combined effect of tunneling in 
FRMs plus nuclear spin diffusion. We stress that, in view of our quantitative 
analysis of the TSSR, the fact that the NSLR and the TSSR are both roughly 
T-independent below 0.8 K does not mean that they have the same origin in 
this regime. Rather, we attribute them to two different mechanisms, both T- 
independent: the quantum tunneling of the electron spin (for the NSLR) and 
the nuclear spin diffusion (for the TSSR). It is only for T > 0.8 that the same 
thermal fluctuations are responsible for NSLR and TSSR at the same time. 

We point out that, in spite of the even quantitative agreement with the simple 
picture of intercluster nuclear spin diffusion as described by Eq. (|3.30() , a puzzling 
point remains. Due to the dipolar field produced by neighboring clusters, ^^Mn 
nuclei in different molecules will have different Zeeman energies, and the typical 
magnitude of the energy spread involved is actually much greater than the inter- 
cluster nuclear dipole-dipole coupling, and should lead to a strong suppression of 
the flip-flop probability. We shall suggest a solution to this problem in the next 
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subsection, when discussing the field dependence of the TSSR rate. 
4.2.2 Field dependencies 

Further insight in the interplay between quantum tunneling and nuclear spin dy- 
namics is provided by the study of the dependence of the NSLR and TSSR rates 
on a magnetic field Bz applied along the anisotropy axis. It is clear from the 
Hamiltonian ()3.14|) that, in the absence of other perturbations, such a longitudi- 
nal field destroys the resonance condition for corresponding electron spin states 
on both sides of the barrier and therefore inhibits the quantum tunneling. In the 
presence of static dipolar fields, i?dipi by studying the tunneling rate as a function 
of Bz one may in principle obtain information about the distribution of longitu- 
dinal -Bdip) since at a given value of Bz there will be a fraction of molecules for 
which -Bdip = —Bz and will therefore be allowed to tunnel just by the application 
of the external bias. We measured the longitudinal field dependence of the NSLR 
rate W{Bz) while shifting the measuring frequency i^{Bz) = + "InBz with 
z/(0) = 230 MHz, in order to stay on the center of the NMR line that corresponds 
to the molecules that are aligned exactly parallel with the applied field. 

The results for the ZFC sample at T = 20 mK are shown in Fig. 14.51 Since for 
a ZFC sample the magnetization is zero, the field dependence should be the same 
when Bz is applied in opposite directions, as is observed. The data can be fitted 
by a Lorentzian with a full width at half maximum lS.Bz = 0.12 T, thus different 
both in shape and in width from the calculated dipolar bias distribution (Fig. 
I3.1()|) . A possible reason for the extra broadening observed in W{Bz) could be 
the presence of misaligned molecules in the sample, combined with quadrupolar 
broadening of the NMR line. The manifold of quadrupolar-split Mn*^^) lines can 
be closely described as a gaussian with total width 2cj^ ~ 2.4 MHz, i.e. ~ 0.23 T 
in field units, thus already broader than W{Bz). Although we follow the center 
of the line for the well-aligned molecules, the misaligned ones (which experience 
a smaller longitudinal field) could still be in condition to tunnel and contribute 
to the NMR signal by means of their quadrupolar satellite lines. An obvious 
source of misalignment is the fact that our sample consists of an ensemble of 
crystallites, thus there is the chance for part of them to be misoriented. In fact, 
the misalignment is in any case unavoidable for fast-relaxing molecules, which 
we consider to be the actual source of quantum tunneling events. It is known 
(cf. ^3.2.4|) that the FRMs with one flipped Jahn- Teller axis also have a local 
anisotropy axis for the cluster spin which deviates ~ 10° from the crystallographic 
c-axis, and the misalignment could be even greater for the FRMs with two flipped 
J-T axes. A confirmation for the role of misaligned molecules in the width of 
W{Bz) comes from the fact that, if we repeat the experiment taking i/(0) = 231 
MHz, thus slightly higher than the center of the line, W{Bz) appears indeed much 
narrower [Fig. I4.6f a)] since the tails of the NMR lines in misaligned molecules 
are now farther away from the measuring frequency. The central frequency for 
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Fig. 4^.5 : Longitudinal field dependence of the NSLR rate W in ZFC sample 
at T = 20 mK. The measuring frequencies are v{Bz) = 230 + jnBz MHz. The 
solid line is a Lorentzian fit with full width at half maximum AB^ = 0.12 T. 

the data in Fig. 14.51 is i/(0) = 230 MHz, which is just shghtly below the center of 
the Une, and this is sufficient to cause the maximum of W{Bz) to be at B^ ^ 

Fig. I4.6f a') also shows a comparison between W{Bz) at T = 20 and 720 mK, 
thus just at the edge of the thermally activated regime. It appears that both 
W^(0) and the width of the resonance increase by raising the temperature. In 
addition, at T = 720 mK an extra small peak appears at B^ — 0.5 T, which 
value coincides with the first level crossing for slow-relaxing molecules. Since the 
peak is not visible at 20 mK, we may attribute it to tunneling of slow molecules 
through higher excited doublets, e.g. Sz = —9 — > +8, which becomes favorable 
as soon as the thermally activated regime is approached f ^3.2.3|) . Moreover, the 
fact that W{Bz) at T = 720 mK and fields < < 0.5 T is largely increased 
as compared to the data at 20 mK, would confirm that the thermal activation in 
FRMs becomes an important source of relaxation already much below 0.8 K, as 
already argued from the analysis of the slope of W{T) in the quantum regime. 

In Fig. I4.6r b) we compare the ZFC data of Fig. 14.51 with a similar field 
scan of W{Bz) at T = 20 mK for an FC sample. We notice that in the FC 
case it is possible to define both positive and negative longitudinal fields, i.e. 
corresponding to whether B^ is applied parallel or opposite to the magnetization 
of the sample. One may observe in Fig. I4.6r b) that the shape of the resulting field- 
dependence is asymmetric, in a way that resembles the asymmetry of the dipolar 
field distribution (Fig l3.ld|) . On the other hand, the width in the FC is much 
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Fig. 4-6 : (a) Longitudinal field dependence of the NSLR rate in ZFC sample 
at T = 20 mK (down triangles) and T = 720 mK (up triangles). The measur- 
ing frequency in these datasets is u = 231 + ^nBz MHz. (b) Longitudinal field 
dependence of the NSLR rate in FC (open circles) and ZFC (full circles) sample 
at T = 20 mK, with measuring frequency v = 230 + ^nB^ MHz. The insets are 
sketches of the electron spin transitions responsible for the observed features. 

larger than in the ZFC sample, contrary to the prediction for i?dip- Moreover, 
as mentioned before, the dipolar field distribution is one order of magnitude 
narrower than the peaks in W{Bz). At increasingly large negative fields (where 
negative means opposite to the magnetization the sample) W{Bz) stabilizes to a 
rather high plateau with a small anomaly around B^ — —0.5 T. We shall suggest 
a possible explanation for the high values of W{Bz) in FC sample at negative 
fields when discussing Fig. 14.191 in ^4.4.21 

The peak at B^ — 0.5 T in ZFC sample at 720 mK and the anomaly at 
Bz — —0.5 T in FC sample suggest that tunneling in slow-relaxing molecules 
may also give a contribution under special conditions like high T or inverted 
magnetization. An important and interesting observation is that, contrary to the 
TSSR (cf. Fig. WM > at zero applied field the NSLR is faster in the ZFC than in the 
FC sample |lllj . Since we argued above that intercluster nuclear spin diffusion 
should be faster for the FC than for the ZFC sample, this would indicate that 
there are more tunneling events in a ZFC sample. As mentioned before, such an 
observation is only possible thanks to the fact that NMR is sensitive to magnetic 
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Fig. 4.7 .■ Decay of transverse magnetization at T = 20 niK in ZFC sample, 
with an applied field = 0.2 T. Notice the chaotic behavior after about 5 ms. 
The solid line is a fit to Eq. \4-^ first part of the decay. 

fluctuations rather than to macroscopic changes in the magnetization. 

We also investigated the field dependence of the TSSR rate, but we encoun- 
tered an unexpected and rather intriguing phenomenon: already for I-B2I ^ 0.05 
T, the decay of transverse magnetization is characterized by a "chaotic region" 
that appears after a few milliseconds, as shown in Fig. 14.71 In other words, 
after an initially well-behaving decay, the echo intensity fiuctuates at random, 
although it always remains below the envelope that can be extrapolated from 
the first part of the decay. In the same way, by measuring several times the 
echo intensity with the same delay between pulses, one would randomly find any 
intensity between zero and the height of the extrapolated envelope. Also, the 
extension of the initial well-behaving decay tends to reduce with increasing field. 
Such a chaotic behavior is observed both in FC and ZFC samples, and both in 
Mn'^"'" and Mn'^"'" sites (see i i4.2.4j) . We stress that this has nothing to do with 
instrumental noise: at any point we see a clean echo signal, just the amplitude 
and the phase are totally random. We don't have an explanation for this phe- 
nomenon, but it's intriguing to remark that 0.05 T is just the maximum spread 
of the dipolar bias distribution (Fig. I3.10j) . 

By fitting just the initial part of the decay we can still tentatively extract 
the field dependence of the TSSR rate T2^{Bz). Unfortunately, on such a short 
interval it is not possible to separate the gaussian and exponential contributions 
to the decay, thus we fitted all the data to Eq. 1)4. 3|) which for the zero-field 
dataset gives just an average value. Given the circumstances it's not easy to 




Fig. 4-8 ■ Comparison between the field dependencies of the NSLR (circles) 
and TSSR (squares) rates, normalized to the zero-field value, in ZFC (left panel) 
and FC (right panel) sample, at T = 20 mK. 



evaluate the error bars, but we may assume them to be rather large. Fig. 14.81 
shows a comparison between W{Bz) and T2^{Bz), both normalized to the zero- 
field value. As it appears, the field dependencies of the NSLR and TSSR rates 
are quite similar, although curiously the similarity seems closer for the positive 
field direction, relative to the FC sample. Moreover, the observation that the 
TSSR is faster in a FC than in ZFC sample remains valid when ^ 0. 

We suggest that there is a possibility to explain at the same time the resem- 
blance of the field dependencies of NSLR and TSSR rates and the efficiency of 
the intercluster nuclear spin diffusion. The tunneling in fast-relaxing molecules 
produces a fiuctuating dipolar field, such that the Zeeman energy of ^^Mn nuclei 
in equivalent sites of neighboring molecules also fluctuates in time; in this way, 
there is a chance that at some instant the nuclei in different molecules will have 
precisely the same Zeeman splitting and will therefore easily undergo a flip-flop 
transition. If we accept this mechanism as responsible for the efficiency of the 
intercluster spin diffusion, then it's clear that by applying a longitudinal field 
the drastic reduction of tunneling rate should also decrease the probability that 
two nuclei in neighboring clusters would come at exact resonance, resulting in a 
reduction of the TSSR rate with the same field dependence as the NSLR rate, as 
indeed observed. 
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Fig. 4-9 ■' Comparison between the nuclear inversion recoveries (a) and the 
decays of transverse magnetization (h) in the "natural" Mni2-ac (circles) and in 
the deuterated sample (squares), at T = 20 mK in zero field and ZFC sample, for 
the Mn^^^ site. The solid lines in (a) are fits to Eq. \4-l^ - 

4.2.3 Deuterated sample 

The role of the fluctuating hyperfine bias on the incoherent tunneling dynamics 
of SMMs, predicted by Prokof'ev and Stamp f ^3.6|) . has been very clearly demon- 
strated by measuring the quantum relaxation in Feg crystals where the hyperfine 
couplings were artificially modified by substituting ^^Fe by ^''Fe or by jl6j . 
For instance, the time necessary to relax 1% of the saturation magnetization be- 
low 0.2 K was found to increase from 800 s to 4000 s by substituting protons 
by deuterium, whereas it decreased to 300 s in the ^^Fe enriched sample. We 
already mentioned in ^3.6.21 that the "tunneling window" for the cluster spin is 
most effectively enlarged by those nuclei that easily coflip with the cluster spin. 
In SMMs, this is indeed the case for the nuclei located at ligand molecules (like 
^H) because their local field is the sum of the dipolar fields from several clusters, 
whereas the ^^Mn or the ^^Fe nuclei are less likely to coflip since the hyperfine 
fields they experience before and after the flip are almost exactly antiparallel. 

In this sense, it is understandable that the largest effect on the tunneling rate 
is seen by isotope substitution of by ^H. The deuterium nuclei have spin 1 = 1 
but their gyromagnetic ratio is 6.5 times lower than for protons. This is the reason 
why, in the framework of the PS theory, one expects indeed a reduction of the 
incoherent tunneling rate due to the weaker (fluctuating) hyperfine bias. Since 
in Mni2-ac the only possible isotope substitution is — >^H, we performed a 
short set of measurements on a deuterated sample. The sample consists of much 
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smaller crystallites than the "natural" ones used in all other experiments reported 
in this chapter. Despite having been prepared following the same procedure as 
described in ^4.11 the orientation of the deuterated sample turned out to be 
almost completely random, probably due to the too small shape anisotropy of 
the crystallites*^^\ We report therefore only experiments in zero external field, 
where the orientation is in principle irrelevant. 

The results are shown in Fig. EHl the ^^Mn NSLR rate at T = 20 mK m zero 
field and ZFC sample is indeed reduced to VFdeut — 0.0035 s~^, i.e. precisely 6.5 
times lower than in the "natural" sample! This clearly points to a reduction of 
the electron spin tunneling rate of the same type as observed in the deuterated 
Feg jn]) and confirms once more that the NSLR is determined precisely by such 
tunneling fluctuations. As regards the TSSR, the results is very intriguing: slow 
but rather ample oscillations are superimposed to the decay of transverse mag- 
netization, somehow reminiscent of what is observed in the natural sample upon 
application of a small field (Fig. 14. 7j) . but without the chaotic features. Fitting 
the average decay constant [Eq. (|0|) ] yields T^^ ~ 38 s"^ i.e. a TSSR ~ 2.5 
times slower than in the natural sample. Such reduction is consistent with the 
hypothesis that the intercluster nuclear spin diffusion is helped by the fluctuation 
of the bias on the nuclear Zeeman levels when tunneling events occurs f M.2.2|) . 
As a matter of fact, both the NSLR and the TSSR in the deuterated sample are 
rather similar to the data in the "natural" sample with an applied field Bz — 0.05 
T. The oscillations in the decay of transverse magnetization remain an obscure 
phenomenon, but observing them even in this zero-field experiment seems to rule 
out instrumental artifacts, and suggests that something bizarre happens in the 
nuclear spin-spin interactions as soon as the tunneling fluctuations slow down 
below a certain threshold. 

4.2.4 Comparison with a Mn'^+ site 

Some rather intersting results come from the analysis of extra measurements 
performed on the NMR line of the Mn^^) site, i.e. a Mn^+ ion, which at low 
temperature is centered at v^"^^ ~ 283.7 MHz. We didn't follow the whole tem- 
perature evolution of the spectrum, but comparing our data with the previous 
experiments |89U87j at higher T it seems that the resonance line has slightly 
shifted to higher frequencies upon cooling. 

Fig |4.10l shows a comparison between the decay of transverse magnetization 
and the recovery of the longitudinal magnetization in Mn(i) and Mn^^) sites, at 
T = 20 mK in FC sample and zero external field. The TSSR is very similar 
in both sites, although a closer inspection evidences that the gaussian nature of 
the decay is less pronounced in the Mn^^^ sites, which leads to = 84 s~^ 
instead of the T^q = 104 s~^ found in Mn^^^. More importantly, the NSLR is 

'^'At room temperature, the shape anisotropy is indeed the only mechanism by which the 
crystaUites may orient in a magnetic field. 
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Fig. 4-iO : Comparison between the decay of transverse magnetization (a) 
and the recovery of longitudinal magnetization (h) in Mn^-^^ (circles) and Mn^"^^ 
(diamonds) sites, at T = 20 mK in FC sample and zero external field. The solid 
(Mn^"^^ ) and dashed (Mn^^^ ) lines are fits to Eg. \4--4^ in panel (a) and Eg. 
in panel (h). 



three times slower in the Mn*-^-* site, as seen in Fig. I4.10r b). This is opposite 
to the high-T regime 

jSniES], where the Mn^+ sites were found to have much 
faster relaxation. Furthermore, the field dependence of the NSLR rate appears 
sharper in the Mn^^^^ site, as shown in Fig l4.1ll The asymmetry in W{Bz) for 
a FC sample is still present, but less evident than in the Mn^^^ site due to the 
more pronounced decrease of W already for small applied fields. These results 
suggest that the tunneling in fast-relaxing molecules is less effective in relaxing 
the nuclei in a Mn^"*" site than in a Mn*^^ site, although the intercluster nuclear 
spin diffusion appears to work comparably well. In fact, the measured TSSR 
rate is mainly determined by intercluster spin diffusion between nuclei in slow 
molecules, since the nuclear spins in FRMs are just a small minority. On the 
contrary, the NSLR rate is determined by the combination of tunneling in FRMs 
and nuclear spin diffusion between fast and slow molecules. Since the FRMs differ 
from the majority clusters precisely in the electron spin arrangement of the Mn^"*" 
sites, it is conceivable that this would be an obstacle for the coupling between the 
nuclei in Mn^+ sites of fast and slow molecules, which are no longer equivalent, 
whereas the Mn'*^ sites remain relatively unaffected by the distortion in FRMs. 
This would explain both the much slower NSLR because of the "broken link" 
between nuclei in fast and slow molecules, and the slightly slower TSSR where 
we just notice the "absence" of the small fraction of nuclei located at the FRMs 
from the participation to intercluster spin diffusion. 

When the temperature is increased such that also the slow-relaxing molecules 
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Fig. ■ Longitudinal field dependencies of the NSLR rates in Mn^^^ (circles) 

and Mn^"^^ (diamonds) sites, normalized at the zero- field value. The data are taken 
at T = 20 mK in FC sample with central measuring frequencies i^^"'^^(0) = 230 
MHz and u^'^\0) = 283.7 MHz. 



start to have considerably fast tunneling fluctuations via higher excited levels 
(T > 0.8 K), the problem of bad nuclear spin diffusion at Mn^+ sites between 
fast and slow molecules becomes irrelevant. We may then expect that the NSLR 
in Mn^+ sites is produced by the combined action of the fluctuating transverse hy- 
perflne field due to intrawell excitations [the model on which Eq. (|3.25|) is based] , 
plus the effect of thermally-assisted tunneling. In Mn^+ sites the first contribution 
should be absent, which could explain why the NSLR in the thermally-activated 
regime is lower than in Mn^"*" sites. 

We may therefore conclude that most if not all the different data presented in 
^4.21 can be consistently explained within our model of tunneling in FRMs plus 
intercluster nuclear spin diffusion. 



4.3 Effects of strong perpendicular fields 

In this section we describe the effects on the NSLR rate of a strong magnetic 
field B±, applied perpendicular to the anisotropy axis. The interpretation of the 
NSLR is tightly related to the NMR spectra, which provide useful information 
about the orientation of the crystallites and the crossover between ^^Mn and 
Larmor frequencies. 
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4.3.1 NMR spectra 

In a single crystal having the c-axis exactly perpendicular to the external field, the 
three ^^Mn NMR lines should evolve as in Fig. 13.91 neglecting the possibility of 
delocalization of the spin wave functions f H3.3.2|) . Our sample consists of a large 
number of crystallites, oriented (within a certain spread) along an axis which is 
orthogonal to both the external field and the axis of the NMR coil. In this way, 
whatever is the canting angle of the electron spins towards the external field, the 
total field at the nuclei is always perpendicular to the axis of the NMR coil. 

Fig. 14.121 shows field-sweep spectra taken at T = 20 mK. The sample can be 
considered as saturated, since we applied 9 T before starting the sweeps. By this 
procedure we expect most of the cluster spins to be aligned in the same direction 
within each crystallite, except eventually for those crystallites that are oriented 
perpendicular to the field within a fraction of a degree. The signal belonging 
to the protons (a) is recognizable by the chopped top; in fact, the line is so 
intense that it saturates the receiver on this scale, and to observe it properly 
we had to attenuate the input by at least 30 dB. Focusing on the Mn^^) line, we 
observe that at low frequencies and fields the signal is characterized by a peak (b), 
plus a more spread-out background that we attribute to misoriented crystallites. 
The peak fits indeed with the calculated position for crystallites with easy axis 
oriented at an angle 0b = 88° — 90° with the field (see below), as explained in 
^3.3.21 By increasing the field, line (b) crosses and is retrieved on the right hand 
side of line (a), but now two peaks are clearly visible. The peak on the right, 
labelled (c), tends to increase in intensity at the expenses of (b) when moving to 
higher fields and frequencies. For instance, the ratio between the amplitudes of 
the (b) and (c) peaks is 0.73 at = 303.2 MHz and 0.67 at v = 311 MHz. An 
extra series of small and sharp peaks (d) is also visible on the right hand side of 
the ^H line for u < 290 MHz. 

In Fig. 14. 131 we report the positions of all these peaks and compare them with 
the calculations. It is then clear that, of the two broad peaks at B± > 7 T and 
v > 295 MHz, the left must indeed be the continuation of line (b) corresponding 
to nuclei in crystallites oriented at 88° from the field. Remarkably, if line (b) 
at low fields contained contributions from crystallites at 9b — 90° (it's hard to 
believe that this would not be the case), they seem to have disappeared at higher 
fields. The right peak (c), of increasing height, could have two origins. On the 
one hand, the calculated spectra for Ob = 70° — 80° tend to accumulate precisely 
around the position of (c), thus the signal intensity that is found in the region 5 T 
- 270 MHz should in principle pile up just around (c) for 7.5 < Bj^ < 8.5 T. On the 
other hand, for B± > 8 T the tunneling splitting becomes so large (Fig. 13. 4 j) as to 
overcome the bias ^ originating from the longitudinal component of the dipolar 
fields, but also that arising from the misalignment of the crystallites as long as 
(90° — 9b) is not too large. One may expect f ^3.3.2|) the NMR signal intensity 
to start being displaced towards the calculated line for delocalized electron spins. 
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Fig. 4-^^ ■ Field-sweep NMR spectra taken at T = 20 mK. Several features, 
marked by letters, are visible and discussed in the text. The high and chopped 
peak observed in the range 6 — 7 T is due to the protons in the sample. 
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Fig. 4-13 : Positions of the peaks in the field-sweep NMR spectra: the symbols 
correspond to lines (h), (c) and (d) as specified in the inset. For the protons (a) 
only the calculated central value is given, plotted as a dashed gray line. Solid lines 
are calculated spectra according to Eq. liiS.18\) . with the indicated orientation of the 
z axis relative to the magnetic field. The dotted line is the calculated spectrum 
for Ob = 90°, taking into account the reduction of the total spin by quantum 
superposition [Eq. iS. 21\) ] 
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beginning with the best ahgned crystaUites and involving more and more amounts 
of sample as B± increases. This may explain why, upon increasing B±, line (c) 
seems to grow at the expenses of line (b), and why there is no trace of signal 
at 9b — 90°. Most likely, a detailed description of the NMR spectra involves a 
combination of the phenomena mentioned above. 

One may ask to what extent can we trust such fine-detailed comparisons be- 
tween measured and calculated spectra, in view of the fact that the parameters 
of the spin Hamiltonian are affected by some uncertainty. We found out that the 
calculations are in fact very robust, since the effect of varying the spin Hamilto- 
nian parameters, even over a range much wider that the accepted error margins 
for the commonly used values, is hardly visible on the scale of Fig. 14.131 

The extra small line (d) is rather fascinating: its position is in fact close to 
the calculated spectrum for a perfectly oriented crystallite, where the total spin 
falls into a symmetric superposition of | ff) and | JJ-) as soon as A > ( ^3.4[ Fig 
13.11(1 . This calculation, contrary to those performed assuming S* = 10 (solid lines 
in Fig. I4.13|) . is extremely sensitive to the values of the parameters of the spin 
Hamiltonian. One could easily allow the dotted line in Fig. 14.131 to be drawn 
just on top of the data, by varying the parameters in the calculation within a 
reasonable margin. To further investigate this issue we measured field-sweep 
spectra on another sample, consisting of deuterated Mni2-ac crystallites having 
an almost random orientations f ^4.2.3|) . Line (d) is still visible but its intensity 
is even smaller than shown in Fig. 14.121 Moreover, it continues linearly down to 
at least 230 MHz, i.e. much lower than expected from the calculation (cf. Fig. 
13.11(1 . without showing any turning point. The most plausible origin of line (d) 
is therefore a set of ^H nuclei very close to the magnetic core and experiencing 
an extra dipolar field ~ 0.6 T; in the deuterated sample, the intensity is further 
reduced because of the substitution of most of the ^H by ^H. The ^^Mn signal in 
the deuterated sample does not yield well-resolved peaks because of the random 
orientation of the crystallites. 

The NMR lines corresponding to nuclei in Mn^"*" sites, whose calculated fre- 
quencies for 6b = 90° are shown in Fig. 13.91 are broader and less intense than the 
Mn(i) line [HHIEZI. The distribut ion of 9b broadens their spectrum even further, 
such that their contribution to the measured signal is basically just a uniform 
background. 

Finally, we remark that all the discussion made so far is relative to ^^Mn 
nuclei in Mn^^) sites of standard, slow relaxing molecules. Nuclei in fast-relaxing 
molecules have their resonance frequency shifting much more quickly to higher 
values, because the lower anisotropy implies that the electron spin can cant more 
easily towards the field. Moreover, the direction of the local hyperfine field in 
FRMs is likely to deviate from the crystallographic c-axis even n zero applied 
field ( p. 2. 4(1 . yielding an extra contribution to the spread of the resonance line. 
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4.3.2 NSLR in perpendicular fields 

In view of the results described in ^4.21 it appeared very interesting to check 
whether the possibihty to increase the tunnehng sphtting by means of a perpen- 
dicular field can give visible features in the NSLR rate W{B±). One would expect 
an acceleration of the tunneling rate and therefore an increase of VF; in reality, 
the problem is made more complicated by the misalignment of the crystallites 
with respect to the field, and by the fact that, at low fields, the tunneling dy- 
namics is provided only by FRMs. Whatever the orientation of a crystallite may 
be, it is very unlikely that many FRMs will be exactly perpendicular to the field, 
since their local anisotropy axis does not coincide with c. As a consequence, one 
would expect any external field to effectively destroy the tunneling dynamics of 
most of the FRMs. 

Fig. 14.141 shows W{B_i_) at T = 20 mK; the values of W are obtained by 
fitting the inversion recoveries to Eq. 1)4. 2() . i.e. by a stretched exponential. As 
expected, as soon as B± > 0.2 T we observe a strong reduction of W, due to the 
suppression of tunneling in FRMs. For B± > 4 T, the NSLR starts to accelerate 
again: at this point A ~ 10~^ K, already larger than in most FRMs, but the effect 
of misalignment can still be very important. At higher fields it becomes essential 
to specify at which position we take the data. The main results (square symbols) 
obey the criterion of maximal signal intensity at a given field: in practice, line 
(b) for S_L < 6.5 T and line (c) for B^ > 6.5 T. 

At B± ~ 6.7 T, the Mn^^) line crosses the protons line, which is 30 dB more 
intense that the ^^Mn line because of the much larger number of protons present 
in each molecule. This fact makes it is very easy to recognize when protons are 
involved, but it also means that the Mn*^^^ signal is in fact completely masked by 
the signal when the crossing of the Larmor frequencies occurs. The measured 
NSLR (crosses in Fig. I4.14j) must therefore be interpreted as the NSLR of the 
protons. 

By increasing B± above 7 T, W{B±) starts to decrease again. As shown in Fig. 
14.151 the tunneling splitting is at this point much higher than both the thermal 
and the nuclear Zeeman energies. Under these conditions, due to depopulation of 
the upper level of the tunnel-split ground doublet, the electron spin fiuctuations 
should indeed slow down, consistent with the observation of a Schottky anomaly 
in the field-dependent specific heat Cm{B±) [TT]lll2j which occurs at the crossing 
between thermal energy and tunneling splitting. Although the transition rate 
between the symmetric ground state and the antisymmetric excited state is here 
measured less directly than in a specific heat experiment, the largely similar field 
dependencies of W{B±) and Cm{B±) for 5 T < B± < 9 T points to a common 
physical origin for the two phenomena. The crossing between Mn^^^ and lines 
also gives rise to an apparent increase in the NSLR, but this is simply due to 
the fact that the protons line is itself characterized by a higher NSLR rate, and 
completely masks the Mn^^^ line. On the other hand, as seen in Fig. I4.12[ the 
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Fig. 4-14 ■ NSLR rate as function of perpendicular field at T = 20 mK, follow- 
ing the highest signal intensity (squares), at the crossing between and ^^Mn 
(crosses), at B± > 7.25 T as function of frequency (triangles), and along line 
(b) (circles). The inset shows the measuring frequency for each set of datapoints 
(same symbols as in the main panel), compared with line (c) (solid), line (b) 
(dashed) and line (a) (dotted). 
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Fig. 4-^5 : Perpendicular field dependencies of the Zeeman energies in 
(solid) and (dashed), the tunneling splitting (dotted) for the indicated orien- 
tations in the xy plane, and the thermal energy (grey shadow). 

crossing range should be limited to the 6.5 — 7 T range and is therefore definitely 
not wide enough to justify the whole peak in W{B±). 

We remark that the mentioned possibility to observe the field dependence of 
the NSLR from time-dependent specific heat experiments has interesting implica- 
tions. In order to measure the hyperfine contribution to Cm, the nuclei must be 
able to absorb heat quickly enough as compared to the experimental time window 
of the specific heat experiment (~ 100 s). Indeed, this condition was observed to 
hold only for S_l > 5 T [TT^ . 

Two extra sets of data points have been measured. One set (circles in Fig. 
I4.14|) is taken along line (b) on the high-field side, i.e. where crystallites with 
9b — 88° are expected to show up, and shows an apparently opposite behavior 
compared to line (c), since W{B±) strongly decreases just where it would go 
towards a maximum on line (c). The other dataset (triangles) is taken at constant 
= 7.25 T while shifting the frequency towards lower values, i.e. moving from 
line (b) to line (c). Indeed W grows along that direction. These observations are 
difficult to understand if we interpret line (c) as the bunching of the signal coming 
from crystallites oriented at 70° — 80° from the field, since with such misalignment 
in the field range 7 — 9 T the two lowest energy levels are separated by at least 10 
K, due to the large longitudinal bias. One would then expect the electron spin 
fluctuations to be totally frozen at T = 20 mK. As mentioned when discussing 
Fig. I4.12| it is also possible that line (c) contains (at least for a part of its total 
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Fig. 4- 16 : Perpendicular field dependence of the stretching exponent a for the 
inversion recovery \4-^ j at T = 20 mK. 

intensity) a signal from a growing fraction of molecules that are very well oriented 
perpendicularly and whose electron spin occupies stably the symmetric ground 
state. The excitations to the antisymmetric state would then represent tunneling 
oscillations of the cluster spin, and, given the magnitude of A in this regime, 
we may even expect such oscillations to be coherent for fairly long times ^21- If 
we adopt this picture, (c) should be a highly inhomogeneous line, consisting of a 
mixture of signals from very well oriented {Ob = 88° — 90°) and badly misoriented 
{Ob < 80°) crystallites. In this case one expects the inhomogeneity of the line to 
show up in the stretching exponent a of the inversion recovery (' ^4.1|) . Fig. 14.161 
shows that a{B±) is indeed a good indicator for the inhomogeneity of the signal: 
it first decreases smoothly up to ~ 6 T, then shows an abrupt dip where 
also the protons participate to the signal, finally it keeps decreasing quite steeply 
along line (c), possibly because of the increasing contribution of nuclei from well 
aligned molecules in their symmetric ground state. 



4.4 Nuclear spin temperature 

A crucial issue in the nuclear spin dynamics of single-molecule magnets is to 
which extent the nuclear spin temperature jll3j is able to follow the lattice tem- 
perature (as determined by the phonons in the molecular crystal), and what is 
the mechanism that links the two systems. To our knowledge there is neither 
experimental nor theoretical work published on this topic so far. 




Fig. ; Comparison between bath temperature Ti,ia,h (solid lines) and nuclear 

spin temperature T^uci (circles), while cooling down the system (a) and while 
applying step-like heat loads (b). The inset illustrates how Tnuci is obtained from 
the NMR data. The waiting time between NMR pulses was 60 s in (a) and 180 
s in (b). Both datasets are at zero field in ZFC sample. 



4.4.1 Time-evolution of the nuclear spin temperature 

We have investigated this problem by coohng down the refrigerator from 800 mK 
to 20 mK while monitoring simultaneously the temperature Tbath of the ^He bath 
in the mixing chamber (just next to the sample) and the NMR signal intensity 
of the Mn(^) line, checked by an echo pulse sequence every 60 s, in zero external 
field and on a ZFC sample. The nuclear spin temperature Tnuci is obtained as 
described in M.ll and plotted in Fig. 14.171 together with Tbath- We find that 
the nuclear spin temperature indeed strictly follows the bath temperature, with 
a small deviation below ~ 200 mK that we attribute to the heating effects of 
the NMR pulses. In fact, the pulse rate used in this experiment is ten times 
higher than in all the measurements shown in ^4.21 so it is likely that not all the 
energy provided to the nuclear spin system can dissipate to the ^He stream at 
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such a rate. We therefore repeated the experiment at the lowest temperatures 
applying a pulse train every 180 s and step-like heat loads. As shown in the inset 
of Fig. KT7\ 

TqucI indeed follows quite accurately Tbath except for a remaining 
discrepancy at the lowest T, but now much less evident than in the previous 
experiment with high pulse rate. 

The application of step-like heat loads is a useful method to investigate the 
effective time constant rth for the thermalization of the nuclear spin system with 
the helium bath, in particular the relationship between rth and the NSLR rate W 
as obtained from the inversion recovery technique and discussed in M.2L This can 
be done by monitoring Tj^^^i in different longitudinal fields and different Mn sites, 
by which we can easily change the NSLR rate. Since also the NMR signal intensity 
changes under each different condition, we must redefine every time the conversion 
factor K between signal intensity and Tnuci- Moreover the monitoring pulse rate 
must be adjusted as well, and we found that it is indeed the combination of pulse 
rate and NSLR rate that determines the discrepancy between the observed lowest 
spin temperature and the base temperature of the refrigerator. Upon heating up 
to, say, 200 mK, this effect becomes irrelevant as seen in the inset of Fig. 14.171 
It is therefore a natural choice to define K in such a way that the asymptotic 
value of Tnuci matches the measured Tbath at the end of the heat step. Fig. 14.181 
shows four examples of the time evolution of Tnuci under the application of a heat 
load for ~ 2 hours, in Mn^^^ and Mn^^^ sites, with or without an applied field, 
and with an increased ^He flow rate. We fitted the data to the phenomenological 
function: 



where Tnuci (oo) is set by definition equal to Tbath at the end of the step, 
Tnuci(O) follows automatically from the above constraint, and to is the time when 
the heat pulse is applied. We find that rth is always much longer than the nuclear 
relaxation time Ti ~ 1/2W, and that the dependence on Mn site and applied 
field is not proportional to the change in Ti under different conditions. What 
really seems to make the difference is the ^He flow rate: within the errors, the 
ratio of heat transfer from the ^He stream to the nuclear spins is proportional 
to the ^He flow rate, given the same conditions of nuclear site and external field. 
This confirms that the thermal contact between nuclear spins and lattice phonons 
is quite fast, as we demonstrated by the values of the NSLR rate, and that the 
main bottleneck is at the interface between the crystallites, the embedding epoxy 
and the '^He stream. 

4.4.2 Effect of large magnetic fields and superradiance 

To conclude this section we discuss an issue which may be relevant to most of 
the experiments on quantum tunneling studied by magnetization loops or fast 



Tnuci (i) = Tnuci (0) + [Tnuci (oo) " Tnuci (0)] 




(4.7) 
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Fig. 4-18 : Time evolution of the nuclear spin temperature (open symbols) and 
the hath temperature (dotted lines) upon application of a step-like heat load. The 
solid lines are fits to Eq. \4- yielding the thermal time constants reported 
in each panel. All data are for an FC sample, with the Mn site, external magnetic 
field Bz, NSLR rate W , NMR pulse waiting time, ^He flow rate and applied heat 
load as given in the annexed table. 
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field-sweeps. While sweeping the magnetic field in such a way that many cluster 
spins are forced to flip, also the hyperfine fields at the nuclear sites are being in- 
verted at an abnormal rate. One may then ask whether the nuclear spin system 
remains in thermal equilibrium, or the macroscopic inversion of the electron spin 
polarization leads to a remarkable increase of Tnuci- To answer this question we 
performed a full magnetization loop — > —8.5 — > -1-8.5 ^ T at a sweep rate 
IdS^/dt] = 0.5 T/min on a sample with initial positive magnetization and mea- 
sured Tnuci as soon as the loop was finished, with the field back to zero. Because 
of the continuous shift of the resonance frequency it would have obviously been 
impossible to monitor T^^d during the field sweep. As shown in Fig. 14.191 we find 
that Tnuci at the end of the field sweep is much higher than Tbath- Interestingly, 
the evolution of Tnuci is characterized by an initially very fast decay, whereas 
after a few minutes Tnuci relaxes to its equilibrium value with precisely the same 
time constant rth — 57 min as already found in the heat transient experiment 
in Fig. I4.18r a). By extracting Tnucl(c«) from the fit to Eq. (|4.7j) and imposing 
that it should be the same also for the initial fast decay, we refitted just the very 
first datapoints obtaining Tnuci (0) ~ 350 mK and ^^^^^^^^^ ~ 300 s. The latter 
value is still a bit longer than Ti = {2W)~^ ~ 40 s in FC sample at T = 20 mK 
[see Fig. I4.1()r b). open circles] but not incompatible with the idea that the very 
first part of the relaxation is the real spin-lattice relaxation, whereas afterwards 
the equilibrated lattice and spins relax to the bath temperature at the usual rate 
mainly determined by the ^He circulation. 

In the inset of Fig. 14.191 we also report and Tbath as a function of time 
during the field sweep. Tbath remains rather stable around 35 mK except at two 
special points, —1.9 T sweeping down and -1-1.9 T sweeping up, where it suddenly 
jumps to ~ 100 mK. We notice that this occurs only when is opposite to the 
instantaneous magnetization, and that in fact \Bz\ = 1.9 T is a level crossing 
field for the electron spins in slow-relaxing molecules, which suggest that the 
phenomenon may be due to magnetic avalanches. A similar heat emission was 
already observed by specific heat experiments jll4j . but at T ~ 2 — 3 K. On the 
other hand, it is rather surprising to find that the main heat emission takes place 
already at B^ = ±1.9 T, i.e. the fourth level crossing, whereas it is known jl06j 
that small Mni2-ac crystals at very low T start to appreciably invert their mag- 
netization only at B^ ~ ±4 T. Moreover we know from experience that, because 
of the large volume of our refrigerator, it is virtually impossible to heat the ^He 
bath so abruptly as it appears in the inset of Fig. 14.191 In fact, such a sudden 
temperature increase resembles very much the effect we find when applying a 
NMR pulse (cf. Fig. 12. 5|) : we first observe a sharp peak in the measured T, 
which is due to the direct heating of the thermometer by the r/-pulse, and only 
in the subsequent minutes we observe a smooth temperature increase, when the 
thermometer indeed senses the heating of the ^He bath. In a recent work jll5j . 
Tejada et al. report the observation of radiation emitted by a large amount of 
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Fig. 4-19 : Time evolution of Tnuci ( open symbols ) measured with NMR pulse 
intervals of 180 s on the Mn^^') line just after a magnetic hysteresis loop performed 
with a field sweep rate \dBz/dt\ = 0.5 T/min. The sample is FC with initial 
positive magnetization. The dashed and solid lines are fits to Eq. l\j for the fast 
(initial) and slow part of the relaxation, respectively, having imposed Tnuci(oo) = 
49 mK for both fits. The inset shows and Tbath during the field sweep: notice 
that the time scale in the main panel is the continuation of the one in the inset. 



Mni2-ac crystals when measuring magnetic hysteresis loops. In their experiment 
the radiation bursts are observed at Bz = ±1.4 T, thus the third level crossing, 
but the measurements are performed at T = 1.8 K, much higher than in our case. 
The proposed explanation is based on the phenomenon of superradiance jll6j : 
when the wavelength of the radiation emitted by each electron spin flip is greater 
than the size of the sample, large numbers of spins may tunnel simultaneously 
because they interact with a common radiation field, thereby producing radiation 
avalanches. A prerequisite for this phenomenon to occur is indeed a population 
inversion, in this case a total electron spin magnetization opposite to the direction 
of the applied field, which is precisely what we observe. 

In fact, even though the thermometer detects only the radiation bursts pro- 
duced by large avalanches, we may expect that weak hints of superradiance could 
occur even at small values of B^, provided that the magnetization is opposite to 
the field. An example is given by the data in Fig. I4.6r b). where we investigated 
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the NSLR rate in a FC sample with "positive" magnetization, and observed in- 
deed that W{Bz) shows a surprisingly high plateau for Bz < 0. An argument 
against the attribution of that plateau to weak hints of superradiance could be 
that we didn't notice any substantial changes in the sample magnetization while 
performing those measurements. On the other hand, we are still speaking about 
relatively low values of W, which would not require an extremely high num- 
ber of tunneling events. A pragmatic solution to this puzzle would be given 
only by continuously repeating the same nuclear relaxation experiment (e.g. at 
Bz = —0.5 T) and checking whether W changes in time, in the same way as the 
total magnetization. 

4.5 Proposal for a unifying theory 

We are now at a stage to gather all the experimental results discussed so far 
and use them as basic elements to build a complete description of the quantum 
dynamics of the coupled system "cluster spins -|- nuclei". The starting frame- 
work is the Prokof'ev-Stamp theory ( ^3.6|) . but we shall supplement it with some 
extra assumptions in order to account for the observed intercluster nuclear spin 
diffusion and the thermal equilibrium between nuclear spins and lattice phonons. 

4.5.1 Basic assumptions 

The model we are going to set up is based upon the following (experimentally 
justified) starting points: 

(i) the tunneling dynamics takes place only in fast-relaxing molecules, whereas 
the neighboring slow molecules can be safely considered as frozen during the 
timescale of interest. For our purpose, the slow molecules serve simply as a 
"reservoir of nuclear polarization" ; 

(ii) for the ease of discussion, we assume that every cluster is coupled to N 
nuclear spins / = 1/2, and that the hyperfine field -B^yp always lies along the 
z axis; upon tunneling, the direction of -Bhyp is turned by 180°. The latter 
assumption simulates indeed the real situation for ^^Mn nuclei in Mni2-ac. We 
also define a local nuclear polarization for the i-th cluster AN{i) = nJ — N^ along 
the -t-i* axis, and the associated local hyperfine bias C^f(^) = /i7ArAA^(i)Shyp " 
Notice that the sign of (,Ni'i) depends on the instantaneous direction of B^yp, 
whereas the sign of AN{i) does not. 

(iii) nuclei in neighboring clusters are coupled by dipole-dipole interactions 
that lead to intercluster nuclear spin diffusion. Experimentally, the spin diffusion 
rate is quite fast: we may therefore assume that any sudden change in local 
hyperfine bias is quickly redistributed among all clusters in order to attain the 
equilibrium within the nuclear spin system. 
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For all other aspects we rely on the Prokof'ev-Stamp model for the descrip- 
tion of the coupling between central spin and nuclei. We assume therefore that 
the hyperfine interactions split the electron spin levels into a manifold of states, 
indexed by AA^ and broadened by intracluster nuclear spin diffusion, spanning 
a total width = ^^^^^^ — ■^^^^ = Nh-f]\fBiiyp. Prom assumption (ii) it fol- 
lows that, since B^yp only takes two antiparallel directions (=^ io^ = 0), the 
hyperfine-split manifolds on either sides of the anisotropy barrier are simply the 
mirror of each other. In the absence of external bias, to any state with AA^*^^-* on 
the left side of the barrier corresponds a state with — AA^(^) on the right side^^\ 
Also, LO'^ = implies that the tunneling of the central spin does not produce any 
coflipping of nuclear spins by the "orthogonality blocking" mechanism f ^3.6.1|) . 
i.e. K ~ 0. Furthermore, the "topological decoherence" is negligible as well; 
the typical nuclear Larmor frequencies are uj^ ~ 10^ s~^, whereas the tunneling 
traversal time f ^3.2.3|) is rtr ~ ^ig ^ ~ 10~^^ s. In the PS language, this leads to 
A ~ N{u;n/^o)'^ — 0. All this means that the only allowed tunneling transitions 
are those that do not require any coflip of nuclear spins, i.e. AA^ = const.; the 
relevant tunneling splitting is therefore Aq. On basis of the assumption (iii) about 
the internal equilibrium in the nuclear spin system we may write the nuclear spin 
temperature Tnuci in terms of the populations of nuclear Zeeman levels. For this 
purpose it is convenient to introduce the populations of the nuclear Zeeman levels 
relative to the local direction of B^yp, A^_|_ and A^_, where A''+ is the number of 
nuclei per cluster that are in the Zeeman ground state, and A^_ counts the nuclei 
in the excited state: 

where ( ) indicates the average over all the clusters in the sample. Calling n = 
N+ - Ar_,(7) {N+) = {N + (n))/2, (N^) = {N - (n))/2 and (^tv) = -^jv(n), 
(|4.8j) becomes: 



Sat + (^Af) _ (Jtwn_ 
Sat - {^n) V ksTancl 

4.5.2 Tunneling mechanism 

In the absence of an external bias on the cluster electron spin, tunneling would 
take place only when (^^r = as well, since the polarization groups with AA^ = 

''^'in H:-{.().ll we used the superscripts and for the hyperfine fields and the nuclear po- 
larizations before and after the fiip, respectively. The notation used here is consistent with the 
above as long as we assume that the the electron spin is in the left well ai t — and subsequently 
tunnels to the right. 

^'^'in other words, n = AA'' if -Bhyp is aligned along +z, and n = — AN otherwise. 
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Fig. 4-20 : Sketch of a possible configuration where tunneling may occur in the 
presence of a dipolar bias ^£,. The gray areas represent the nuclear polarization 
groups on opposite side of the harrier having the same AN , such that tunneling 
does not require any nuclear coflips. In order to populate one of those groups, 
an initial hyperfine bias = Co must be achieved via intercluster nuclear spin 
diffusion. 

are the only ones having equal energy and equal nuclear polarization. In the 
more general case where we include the quasi-static dipolar bias arising from 
neighboring clusters, the tunneling transition with AA^ = const, can still take 
place, but requires now an initial hyperfine bias ^at = ^d- After the tunneling, 
since -Bhyp is now reversed, the hyperfine bias becomes — ^d! in other words, the 
effect of tunneling is in fact to convert the hyperfine energy into dipolar energy 
or vice versa. A sketch of the situation is depicted in Fig. 14.201 starting from 
a cluster where the electron spin is favorably aligned with respect to the dipolar 
field produced by the neighbors, a tunneling transition requires "charging" the 
local nuclei up to polarization group where £,N = S,D, such that on the opposite 
side of the barrier there is another polarization group having identical AN. After 
tunneling, the cluster spin is oriented opposite to the local dipolar field, and the 
required energy has been taken from the nuclear reservoir, which is now in a 
favorable polarization state with respect to the new direction -Bhyp- 

It is clear that intercluster nuclear spin diffusion is an essential ingredient 
to attain the condition ^tv = at all. Flip-flop processes between nuclei in 
different clusters conserve the total energy of the nuclear system, but allow the 
local hyperfine bias to wander through the whole hyperfine-split manifold until 
the tunneling resonance is found. A tunneling event like in Fig. I4.2UI would 
have the effect of "locally extracting" an energy S^£) from the nuclear spin system; 
by the same mechanism of intercluster spin diffusion, this energy can then be 
redistributed among the neighbors, thereby attaining an effective reduction of 
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the nuclear spin temperature. On the other hand, the tunnehng process can 
work in the opposite direction, reducing the dipolar energy at the expenses of an 
increase of Tnuci • 

Let us recall that this picture, in particular as a consequence of assumption 
(ii), accounts only for nuclei like ^^Mn in Mni2-ac, i.e. nuclei where the hyperfine 
fields flips by 180° when a tunneling event occurs. We have shown in ^4.2..Sl that 
actually the most important role in triggering the incoherent tunneling dynamics 
is played by nuclei such as protons, precisely because they have a large coflipping 
probability due to the non-antiparallel directions of and ^^yp- In that case, 
the condition AA^ = const, may be relaxed to |AA^(i) - AiV(2)| < 2M, where 
M is the number of coflipping nuclei. Although such an extension is essential to 
properly account for the tunneling rate of the electron spin, we shall not include 
it in our discussion as long as we focus on the tunneling-driven nuclear relaxation 
of ^^Mn (which is the main subject of our experiments) rather than the hyperfine- 
driven tunneling rate of the electron spin. 



4.5.3 Coupling with lattice phonons 

So far we have not introduced any thermostat to set the temperature to which 
the equilibrium Tnuci should relax. The only natural way to obtain a certain 
equilibrium nuclear spin temperature from the tunneling mechanism described 
above, is to assume detailed balance between the rate of the transitions that 
reduce Tnuci at the expenses of dipolar energy (like in Fig. I4.2()|) . and the rate 
for transitions in the opposite direction: 

4 = e~«^/'=«^. (4.10) 

On the other hand, as long as the tunneling rates are derived from the standard 
Prokof'ev-Stamp expression 1)3. 44(1 . there's no way to distinguish between and 

We propose therefore that, in order to account for the observed equilibrium 
between nuclear spin and lattice temperatures, the PS theory must be extended 
to include the possibility of inelastic tunneling events, accompanied by creation 
or annihilation of lattice phonons. In particular, we argue that an essential role 
could be played by the time-dependence of the dipolar bias that arises when 
the distance between neighboring clusters is modulated by a phonon. Due to 
the peculiar structure of single-molecule magnets, containing very rigid clusters 
bound by soft organic ligands, we may expect a clear separation between the 
intracluster, high-energy phonon modes, and the intercluster low-energy modes, 
corresponding to the relative displacement of the clusters considered as rigid 
objects. This approach has been already successfully used to account for the 
Mossbauer recoil- free fraction in metal cluster molecules |117| . which are bound 
by similarly soft ligands, yielding a Debye temperature @d ~ 10 — 20 K for 
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intercluster phonon modes. In some sense, this idea is a "molecular version" of the 
well-known Waller mechanism pTlRj for the spin-lattice relaxation of paramagnetic 
ions. The problem is that, if used to estimate the nuclear spin-lattice relaxation 
rate in the standard way, i.e. by calculating the spectral density of the fluctuating 
dipolar field directly at the nuclear sites |119j , the Waller mechanism would lead to 
extremely long relaxation times. What we suggest is that the Waller mechanism 
can be regarded as a source of fluctuating bias producing a similar effect as the 
nuclear spin diffusion in the PS theory (i jS.G.ll) : similar conclusions on the effect 
of the Waller mechanism could be drawn by using the Landau-Zener formalism 
(see ^3.2.31 and H3.6.3() to calculate the probability of incoherent tunneling while 
the dipolar bias sweeps through the tunneling resonance. We already remarked 
in ^3.6.31 that the tunneling rate t^^ obtained in this way is independent of the 
bias sweep rate, in this case the frequency of the phonon modulating t^^ 
does depend on the range of fluctuation hT^, whose role would here be taken 
by the amplitude A^/j of the change in dipolar bias arising from lattice strain 
associated with the phonon modes. By inspecting Eq. 1)3. 44() one finds that the 
prefactor oc A^/F^ — > A'^/A^d in r^^ is inversely proportional to the amplitude 
of fluctuation. On the other hand, the extra factor e"'^'/^" may lead to at least a 
partial cancellation effect if A^/j contributes to .^o as well, which could justify the 
approximately T-independent NSLR rate found below 0.8 K. Since the fluctuation 
of dipolar bias has to be combined with the fluctuating hyperfine bias, deriving 
the exact expression for the inelastic tunneling rate constitutes a very challenging 
but formidable theoretical task, that largely exceeds the scope of this thesis. 

4.5.4 Nuclear spin-lattice relaxation rate 

Irrespective of the precise expressions for and w^, their deviation from the 
average value t^^ = {w^ + w^)/2 will be very small, unless the temperature 
becomes much lower than the typical range of the dipolar interaction ~ 0.1 K. 
The reason why the relevant energy scale is and not Hujm is that, in our model, 
the nuclear spin system can change its global polarization only by converting 
hyperfine energy into dipolar energy (^jy = (,d)- In the extreme case where 

= in all clusters, all tunneling events would take place at AA^ = 0, with 
no effects on the nuclear spins^^^. In fact, it would be very intriguing to find a 
system where HuJiy » S^d (e.g. by diluting the magnetic clusters) and check if the 
thermalization of the nuclear spins works as efficiently as we find here. 

Furthermore, by focusing on the NSLR rate W(T) only in the quantum regime, 
we have noticed that the "roughly" T-independent plateau has not only a slope, 
but also a slight and broad maximum around T ~ 0.18 K, below which W{T) 
tends to decrease somewhat more steeply. As shown in Fig. 14. 211 by measuring 
W{T) in the quantum regime in two independent runs, we got the impression 

'^'This argument obviously does not hold for nuclear spins like protons, which have a non- 
negligible coflipping probability. 
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Fig. 4-21 : NSLR rate in the quantum regime, for Mn^^'^ site in ZFC sample and 
zero external field. Circles are the same data as in Fig. \4-^ squares are from 
an experiment carried out in a different run several months before: the arrow 
indicates the position of the weak maximum in W{T) around 180 mK. 



that the maximum around 0.18 K is a weak albeit reproducible feature. The 
value of 0.18 K may not be a coincidence in view of the fact that, as we discuss 
in chapter V, a system like Mng orders ferromagnetically at Tc = 0.16 K purely 
by the effect of dipolar interactions. It would be very interesting to see whether 
a fully developed theoretical description can account for this peak in W{T), and 
if this is indeed related to the dipolar interaction energy. 

An estimate of the NSLR rate in the quantum regime can be obtained by 
writing a master equation for the populations of the nuclear Zeeman levels relative 
to the local hyperfine field direction, and N- (cf. ij4.5.1() . For simplicity, 
since the internal equilibrium is reestablished within a time T2 <C tj- after each 
tunneling event, we assume that just before tunneling all clusters have the same 
values of and A^_, neglecting fluctuations around the mean values. If a fast- 
relaxing molecule tunnels at time t, the polarization of its own nuclei is abruptly 
inverted. Each time a tunneling transition lowers the energy of the local nuclei, 
which occurs at a rate wi, then nuclei have been added to the total number 
of nuclei in the Zeeman ground state. After a time T2 this decrease in local 
hyperfine bias has been redistributed over the sample: calling x the fraction of 
FRMs over the total, then the tunneling event has increased A^-|- to + xN-. 
The same reasoning holds for transitions that increase the hyperfine bias. The 
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master equation is therefore: 



dt 



xN^wi - xN+w-^. (4.11) 



From here the NSLR rate can be obtained by standard textbook calculations 
[nnj . Writing N+ = {N + n)/2 and iV_ = (iV - n)/2, (imi) becomes: 

— =xN{wi-w^)-xn{wi + w^), (4.12) 
which can be rewritten as: 

^ = 2W{no-n), (4.13) 



dt 



no = N^^, (4.14) 



(4.15) 



where uq is the equilibrium nuclear polarization and W is the desired NSLR rate, 
since the solution of (|4.13|) is precisely of the form n{t) = n(0) — [n(0) — ?io][l — 
exp(-2m)]. 

Comparing Eq. (|4.15|) with the experimental value of ~ 0.03 s~^ reported 
in ^4.2.11 one would deduce that, for instance, 1% of FRMs tunneling at a rate 

^ ~ 3 s"^ is sufficient to explain the data. Such a tunneling rate is high 
but not unreasonably high, especially considering the FRMs having two flipped 
Jahn- Teller axes and a barrier ~ 15 K ( ^3.2.4|) . We recall that for an easy-axis 
spin 5 S> 1, described by a Hamiltonian Ti = —DSl + ES"^ with D ^ E, the 
tunneling rate is proportional to D{E/D)^ _58^. By simply reducing the height 
of the barrier by a factor 4, we may already expect an increase of five orders of 
magnitude in the tunneling rate with respect to a normal molecule! FRMs are 
also likely to have higher values of E, thereby increasing the tunneling rate even 
further. 

Let us remark that the interplay between electron spin (tunneling) fluctuations 
and nuclear spin diffusion is completely different from the standard picture of 
relaxation by paramagnetic impurities jl21j . In the presence of an impurity spin, 
the closest nuclei are relaxing very fast because of the strong fluctuations, but 
the static component of the field produced by the impurity drives them out of 
resonance with respect to the nuclei that are farther away: one finds that nuclear 
spin diffusion starts to work only outside a sphere whose radius is determined by 
the strength of the nuclear dipolar coupling compared with the field produced by 
the impurity. In this sense, intercluster nuclear spin diffusion in SMMs is a very 
peculiar phenomenon, since all the nuclei are equivalent and there is no minimum 
radius for the spin diffusion. 
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4.6 Conclusions 

We now briefly summarize the results presented in this chapter, in order to place 
them in a complete framework of observations and interpretations, and in par- 
ticular to make clear what the breakthroughs of our research are. 

(i) Incoherent quantum tunneling of the electron spin provides a relaxation 
channel for the nuclear spin system, and its effects are much more pronounced 
than expected TA . This can be deduced from the roughly T-independent plateau 
in the ^^Mn NSLR rate below 0.8 K, and from the coincidence between maxima 
in the NSLR rate and tunneling resonances for the electron spins. The observa- 
tion of the effect of tunneling on the nuclear spins can be done even in zero field 
and demagnetized sample, thus without any changes in the macroscopic magne- 
tization of the cluster spins. We also suggest that the correct interpretation of 
the observed NSLR for T > 0.8 K should include the effect of thermally-assisted 
tunneling as well. 

(ii) The intercluster nuclear spin diffusion is fast on the timescale of the NSLR 
and of the electron spin tunneling rate. Its existence can be inferred from the T- 
independent value of the TSSR rate below 0.8 K, and from the ~ \/2 ratio between 
the TSSR rates in a fully polarized and a demagnetized sample. This is the 
first direct experimental observation of such a phenomenon in SMMs. Moreover, 
the field dependencies of NSLR and TSSR suggest that the fluctuating dipolar 
field produced by tunneling events may play a role in facilitating the intercluster 
nuclear spin diffusion. 

(iii) The tunneling dynamics in zero external field is ascribed to the fast- 
relaxing molecules present in any real sample of Mni2-ac. This follows from both 
a quantitative analysis of the NSLR rate, and from the experimental observation 
that the macroscopic magnetization of the sample does not change during the 
measurements. Furthermore, the comparison between NSLR and TSSR in Mn^"*" 
and Mn^+ sites indicates that the intercluster nuclear spin diffusion between 
FRMs and normal molecules is less efficient between nuclei in Mn'^"*" sites, which 
are precisely those where the FRMs are different. 

(iv) The isotope substitution of by leads, as expected, to a significant 
decrease of the incoherent tunneling rate, which manifests itself in a congruent 
reduction of the NSLR and TSSR rates. 

(v) The application of strong perpendicular fields produces a broad peak in the 
NSLR rate which is consistent with the observed increase of excitations within the 
tunneling-split ground doublet, as previously observed by field-dependent specific 
heat experiments |llj . 

(vi) Even in the regime (T < 0.8 K) where the nuclear spins are relaxed 
only by ground state quantum tunneling of the electron spins, the nuclear spin 
temperature remains in equilibrium with the bath temperature. This is a new 
and unexpect result, and its implication have never been included in any theo- 
retical discussion. We attempt to account for this observation by a suggesting a 
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model that combines the Prokof 'ev-Stamp theory with the effect of modulation of 
the intercluster dipolar fields (Waller mechanism), in order to elucidate how the 
cluster spin tunneling can relax the nuclear spins, and why it is necessary to at- 
tribute different rates to tunneling transitions that increase rather than decrease 
the nuclear spin energy. 

(vii) Experiments where a longitudinal magnetic field is applied opposite to 
the direction of the sample magnetization could be interpreted by assuming the 
presence of magnetization avalanches, which also manifest themselves in the form 
of electromagnetic radiation jll5j . 
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Few examples of long-range magnetic order induced by purely dipolar interactions 
are known as yet |122|I123| . Therefore, the possibility to study such phase tran- 
sitions and the associated long-time relaxation phenomena in detail in high-spin 
molecular cluster compounds, with varying crystalline packing symmetries and 
different types of anisotropy, presents an attractive subject ^124, 125 . However, 
for the most extensively studied molecular clusters sofar such as Mni2, Fes and 
Mn4 jl2Hpi27pi()9pil2j the uniaxial anisotropy experienced by the cluster spins is 
very strong. Consequently, the electronic spin-lattice relaxation time Tf becomes 
very long at low temperatures and the cluster spins become frozen at tempera- 
tures of the order of 1 K, i.e. much higher than the ordering temperatures Tc ~ 0.1 
K expected on basis of the intercluster dipolar couplings |124pi25pi28j . Although 
quantum tunneling of these cluster spins has been observed jl26lll27lll091lll^ . 
and could in principle provide a relaxation path towards the magnetically ordered 
equilibrium state [12411125] . the associated rates in zero field are extremely small 
(< 100 s~^). For these systems, tunneling only becomes effective when strong 
transverse fields B± are applied to increase the tunneling rate. Although the 
tunability of this rate and thus of Tf by B± could recently be demonstrated for 
Mni2, Feg and Mn4 jll2j . no ordering has yet been observed. 

The obvious way to obtain a dipolar molecular magnet is thus to look for a 
high-spin molecule having sufficiently weak magnetic anisotropy and negligible 
mier-cluster super-exchange interactions. Here we report data for 
Mn604Br4(Et2dbm)6, abbreviated as Mng j85|, whose net magnetocrystalline 
anisotropy (cf. ^l-{.2.5j) proves to be sufficiently small to enable measurements 
of its equilibrium magnetic susceptibility and specific heat down to our lowest 
temperatures (15 mK). 

In this chapter we report the observation that the magnetic clusters do undergo 
a transition to a long-range ferromagnetically ordered state at Tc = 0.161(2) K. 
The transition can be observed both by ac-susceptibility and by specific heat 
experiments, thanks to the fact that the magnetic relaxation remains very fast, 
even at very low temperatures, as compared to the intrinsic experimental times 
(0.1 — 5 ms for ac-susceptibility and 1 — 100 s for specific heat). The magnetic 
relaxation has been investigated by field-dependent specific heat and ^^Mn NMR. 
By combining the two techniques we are able to study the interplay between 
electron- and nuclear- spin-lattice relaxation. 
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5.1 Long-range ferromagnetic ordering 
5.1.1 ac-susceptibility 

The ac-susceptibility data in Fig. 15.11 were taken using the susceptometer de- 
scribed in ^2.3[ varying the frequency v between 230 to 7700 Hz and recording 
the data while slowly cooling down the refrigerator. Additional measurements 
at T > 1.8 K were performed in a commercial SQUID susceptometer. For this 
zero-field experiment, the susceptometer was mounted in the upper Araldite pot 
of the dilution refrigerator (cf. Fig. 12. 1|) . after having replaced the Kapton tail 
by a closed plug. This has the advantage of allowing a higher ^He circulation 
rate, and assuring a the perfect coincidence of the sample and thermometer tem- 
peratures even for T > 0.8 K, when some rather sudden temperature jumps may 
occur. 

The real part x' shows a sharp maximum at Tc = 0.161(2) K. We found 
that x' ^-t Tc is close to the estimated limit for a ferromagnetic powder sample, 
l/{pNs + PsamNc) ^ 0.14 lb 0.02 emu/g, where p = 1.45 g/cm^ and psam — 0.45 
g/cm^ are respectively the densities of bulk Mng and of the powder sample, 
Ng = Att/2> is the demagnetizing factor of a crystallite, approximated by a sphere, 
and Nc — 2.51 is the demagnetizing factor of the sample holder. This indicates 
that Mng is ferromagnetically ordered below Tc- 

The temperature Tpeak ^-t which the maximum value of x' is found, as shown 
in detail in the inset of Fig. 15.11 is seen to vary only weakly with z/, which 
we attribute to the anisotropy. The total activation energy of Mng amounts 
to DS^ ~ 1.5 K, i.e. about 45 times smaller than for Mni2. Accordingly, 
one expects the superparamagnetic blocking of the Mng spins to occur when 
T ~ TB(Mni2)/45, that is below ~ 0.12 K. Since this value is close to the ac- 
tual Tc, one may expect that for T ^ Tc the approach to equilibrium begins 
to be hindered by the anisotropy of the individual molecular spins. We stress 
however, that the frequency dependence of x' observed here is quite different 
from that of the well known anisotropic superparamagnetic clusters. A way to 
quantify the frequency dependence of the peak in x' is by means of the parameter 
^^peak/[^peak^(log z^)], which gives the variation of Tpeak per decade of frequency. 
We find here ATpeak/[TpeakA(log i/)] ~ 0.04, to be compared with the typical val- 
ues of ~ 0.20 for superparamagnetic blocking. In fact it is closer to the value 
~ 0.06 found for certain types of spin glasses |129j . but the peak found here is 
much higher and sharper. 

Below Tc, x' decreases rapidly, as expected for an anisotropic ferromagnet 
in which the domain-wall motions become progressively pinned. The associated 
domain-wall losses should then lead to a frequency dependent maximum around 
Tc in the imaginary part, x") as seen experimentally indeed [cf. Fig. I5.ir b)]. 
In fact, although the Mng spins can be considered as nearly isotropic at high 
temperatures, the anisotropy energy is still large compared with the dipolar in- 
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7(K) 

Fig. 5.1 : Real (a) and imaginary (b) parts of the ac- susceptibility at the 
indicated frequencies. Inset: magnification of x'(T) to evidence the frequency 
dependence of the peak. 
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Fig. 5.2 : T -dependence of the inverse of the real part of the ac- susceptibility. 
Circles: low-T susceptometer, v = 7700 Rz. Squares: SQUID susceptometer. 
The inset shows a magnification of the temperature range around Tc- 

teraction energy /i^/r^ ~ 0.1 K between nearest neighbor molecules. Thus the 
ordering should be that of an Ising dipolar ferromagnet. 

Susceptibility data Xi corrected for the demagnetizing field {xi = x' /[^~{p^s + 
Psam-^c)x']) follow the Curie-Wciss law x = C/{T — 9) down to approximately 0.3 
K, with C = 0.034(1) cm^K/g and 9 = 0.20(3) K. The constant C equals, within 
the experimental errors, the theoretical value for randomly oriented crystals with 
Ising-like anisotropy NAg'^fJ,%S{S + l)/3kBPni = 0.0332 cm^K/g, where S = 12, 
g = 2, and the molecular weight Pm = 2347.06. The positive 9 confirms the 
ferromagnetic nature of the ordered phase. From the mean-field equation 9 = 
2zJeffS{S + l)/3kB, we estimate the effective inter-cluster magnetic interaction 
Jeff ^ 1.6 X 10"'^ K, and the associated effective field H^ff = 2zJeffS/gfiB = 
3.5 X 10^ Oe coming from the 2; = 12 nearest neighbors. 

5.1.2 Specific iieat 

Additional evidence for the ferromagnetic transition at Tc is provided by the 
electronic specific heat Cg, shown in Fig. 15.31 [3l| . The specific heat of a few 



5.1. Long-range ferromagnetig ordering 



111 



milligrams of sample, mixed with Apiezon grease, was measured at low-T in a 
home-made calorimeter [34^1 12) that makes use of the thermal relaxation method. 
An important advantage of this method is that the characteristic time Tg of the 
experiment (typically, Tg ^ 1 — 100 seconds at low-T) can be varied by changing 
the dimensions (and therefore the thermal resistance) of the Au wire that acts as 
a thermal link between the calorimeter and the mixing chamber of the dilution 
refrigerator. 

The electronic specific heat was obtained by subtracting from the total spe- 
cific heat the contribution of the lattice, which follows the well-known Debye 
approximation for low temperatures Qattice oc {T/Qd)^, where ~ 48 K, as 
well as the contribution Cnuci arising from the nuclear spins, as discussed in ^5.2.11 
below. Ce(T) reveals a sharp peak at 0.15(2) K, indeed very close to the peak 
in the real part of the ac-susceptibility. Numerical integration of Ce/T between 
0.08 K and 4 K gives a total entropy change of about 2>AkB per molecule, which 
is indeed very close to the value ks ln(25 + 1) = 3.22/cb for a fully-split 5" = 12 
spin multiplet. We may therefore attribute the peak to the long-range order of 
the molecular spins. We note, however, that at Tc the entropy amounts to about 
IkB per spin, showing that only the lowest energy spin states take part in the 
magnetic ordering. 

5.1.3 Monte-Carlo Simulations 

In order to simulate the zero-field specific heat data, Monte-Carlo (MC) calcula- 
tions were performed by J. F. Fernandez (univ. of Zaragoza, Spain) for an S" = 12 
Ising model of magnetic dipoles on an orthorhombic lattice with axes = 15.7 
A, ay = 23.33 A, and = 16.7 A, which approximates the crystal structure of 
Mng. The Hamiltonian includes dipolar interaction term as well as the anisotropy 
term —DS^ given in Eq. (|3.15|) . When dipolar interactions are neglected, the 
zero field splittings of the S = 12 multiplet produced by this crystal field term 
lead to the Schottky anomaly in Cg shown by the dotted curve in Fig. 3. The fit in 
the range above 0.5 K yields D/ks — 0.013 K, in good agreement with our above 
estimates. For T < 0.5 K the intermolecular dipolar interactions come into play 
and remove the remaining degeneracy of the lowest-lying | it m) spin doublets. 
The MC simulations show that the ground state is ferromagnetically ordered, as 
observed, and predict a shape for Ce that is in very good agreement with the 
experiment. In Fig. 15.31 we show Cg calculated assuming all molecular easy (z) 
axes to point along o^, i.e. one of the two nearly equivalent short axes of the 
actual lattice. Similar results were obtained for other orientations chosen for the 
anisotropy (z) axis. We note that the Ising simulations give Tc = 0.22 K, which 
is slightly higher than the experimental = 0.161(2) K. This difference may be 
due to the Ising approximation taken for the intercluster dipolar interaction. 

In order to pursue this point further, additional MC simulations were per- 
formed for the same crystal lattice, but with classical Heisenberg spins replacing 



112 



Long-range dipolar ordering in MNg 




r(K) 

Fig. 5.3 : T -dependence of the electronic specific heat at zero applied field 
(circles). The dotted line is the Schottky anomaly calculated with D = 0.013 K. 
The full line is the Monte Carlo (MC) calculation for an orthorhomhic lattice of 
1024 Ising spins with periodic boundary conditions. For each point, 2 x 10^ MC 
steps per spin were performed. 

the {S = 12) Ising spins. To investigate the sensitivity to the type of anisotropy, 
the sign and magnitude of D were varied. These calculations resulted in the 
phase diagram shown in Fig. 15.41 which we include here since it illustrates the 
complicated way in which the nature of the actual ground state and the value 
of Tc may depend on the combination of the long-range dipolar interaction and 
the anisotropy parameter. Although the ground state for this lattice is always 
ferromagnetic, it can be either "uniaxial" , with strong preference for the ay axis, 
or "planar", in the sense that the — plane becomes an easy plane with 
a weak preference for a given direction in the plane, as indicated in the figure. 
Interestingly, the switching point between these orientations is not at D = but 
at L> = 3 mK. The reason for this is as follows. Because the crystal lattice un- 
der consideration is far from being cubic, the dipolar interaction energy is rather 
anisotropic. The energy is minimized when the magnetization M points in the 
direction (in the — flz plane) shown in the inset, on the left-hand side of Fig. 
15.41 Therefore, M points along this direction for D < 0. On the other hand, the 




energy minimization for > is a competing process. Clearly, dipolar inter- 
action must become dominant for sufficiently small values of D. The numerical 
results show that this occurs if < D < 3 mK. The numerical datapoints in 
Fig. 15.41 also show that Tc varies sharply within the —0.01 < D < 0.01 K range. 
Because of the discreteness of the MC simulations we can not tell whether Tc 
vanishes at -D ~ 3 mK. The lowest numerical value obtained is as small Tc — 0.03 
K at D ~ 3 mK. Outside this narrow range of D, Tc is already almost equal to 
the limiting values of ~ 0.7 K and ~ 0.3 K, reached for infinite negative and 
positive D, respectively. Such a variation of Tc with anisotropy, as well as the 
form of the calculated and observed specific heat ordering anomaly, appear to 
be specific for dipolar interactions. They differ widely from the corresponding 
behavior known for 3-dimensional (Heisenberg, Ising, XY) ferromagnetic lattices 
with nearest-neighbor interaction only. For instance, in those models the varia- 
tion of Tc with anisotropy is restricted to about 20% and the Tc value is highest 
for the Ising case. 
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Fig. 5.5 : T -dependence of the magnetic specific heat at the specified tem- 
peratures. The lines are the calculated electronic (dotted) and nuclear (dashed) 
contributions. The solid lines are the total (electronic + nuclear) Cm, calculated 
accounting for the nuclear Ti (see text). Inset: detail of Cm{T) at B = 2 T, for 
long (^ 100 s) and short (^ 1 s) experimental times. 



5.2 Spin-lattice relaxation 



The electron- and nuclear- spin-lattice relaxation in Mng is an interesting topic 
to investigate and compare with the typical behavior of the highly anisotropic 
SMMs like Mni2-ac. As we shall see, provided that a Zeeman splitting for the 
electron spin is created by the application of an external field, the nuclear spin 
dynamics in Mng closely resembles the behavior of the nuclei in Mni2-ac in the 
thermally activated regime, where the electron spin energy levels are split by 
crystal-field effects instead. 
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5.2.1 Field-dependent specific heat 

We begin by discussing the specific heat data obtained in varying magnetic field 
B pij, plotted in Fig. 15.51 Even for the lowest B value, the ordering anomaly 
is fully suppressed, as expected for a ferromagnet jl3Uj . Accordingly, we may ac- 
count for these data with the Hamiltonian (|3.15|1 neglecting dipolar interactions. 
The Zeeman term splits the otherwise degenerate | it m) doublets, and already 
for B ~ 0.5 T the level splittings become predominantly determined by B, so 
that the anisotropy term can also be neglected. As seen in Fig. 15.51 the calcu- 
lations performed with D = reproduce the data quite satisfactorily at higher 
temperatures (dotted curves). 

However, when the maxima of the Schottky anomalies are shifted to higher 
T by increasing B, an additional contribution at low T is revealed. It is most 
clearly visible in the curves for lT<i3<2.5T, and varies with temperature as 
cT'^/R = 4 X 10~^. We attribute this to the high-temperature tail of the specific 
heat contribution Cnuci arising from the Mn nuclear spins (I = 5/2), whose energy 
levels are split by the hyperfine interaction with the Mn^^ electronic spins s. 
This interaction can be approximated by 7ih / = AI ■ s, where A is the hyperfine 
coupling constant (cf. ^3.3j) . At high temperatures (i.e., when As <^ ksT) 
Cnuci/R - ^A'^s'^I{I + l)r-2 [ni- Taking A = 7.6 mK as used previously to 
simulate ESR spectra measured on a Mn4 cluster [Wj . we obtain the dashed line 
shown in 15.51 This contribution was subtracted from the zero-field data shown 
in Fig. E21 

A remarkable feature of the experimental data that is not reproduced by these 
calculations is that, at the lowest T, the nuclear specific heat drops abruptly to 
about 10^^ J/gK. The temperature T„ where the drop occurs depends on B but 
also on the characteristic time constant Tg of our (time-dependent) specific heat 
experiment: the deviation from the (calculated) equilibrium specific heat is found 
at a lower T when the system is given more time to relax, as shown in the inset 
of Fig. 15.51 We conclude that the drop indicates that the ^^Mn nuclear spins can 
no longer reach thermal equilibrium within time Tg. We may write: 

Cnuclire) = C^nudi^ " exp(-re/ri )] , (5.1) 

showing that the transition should occur when the nuclear spin-lattice relaxation 
time Ti becomes of the order of Tg. These transitions to non-equilibrium provide 
therefore direct information on the temperature and field dependence of Ti . 

As already discussed in ^3.5.11 Ti{T,B) can be related to the fluctuation of 
the transverse hyperfine field, as produced by the phonon-induced transitions (at 
a rate rs_ph) between different levels of the electronic spin. In the case of Mng 
these are simply the Zeeman levels of the 5 = 12 cluster spin in the magnetic 
field B, split by an energy AE = gfisB. 

At low T and high B, only the ground and the first excited states, m = +12 
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Fig. 5.6 : Field- dependence of the magnetic specific heat at T = 0.3 K. The 
dotted line is the calculated electronic specific heat, the dashed line is the equilib- 
rium nuclear specific heat, and the solid line is the calculated total Cm, accounting 
for the nuclear Ti. Inset: Bn{Tn) obtained from T-sweeps, as in Fig. 15.51 (full 
circles), or from B -sweeps, as in the present figure (open circles). From the linear 
fit (solid line) we extract the average tq — 3 x 10^^ s. 



and m = +11, need to be considered. We may therefore rewrite Eq. (|3.25j) as^^-*: 

where rj~^ = (A6^)/(2i?^Q^rii) incorporates the hfetime 

1 - cxp{-gfiBB/kBT) 
Til = 7 ^ (5-3) 

of the first excited state m = +11, and the ratio between the average value of 
the transverse component of the fluctuating hyperfine field Ab± and the square 



(1) 



Recall that l/Ti = 2W. 
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of the total field at the nuclei -Btot- The approximated expression of (|5.2() is valid 
when UNTs-ph 3> 1, which is practically always the case, except at very low field. 

It is easy to see from Eq. 1)5.2(1 that the nuclear spins can be taken out of 
of equilibrium either by decreasing T down to T„ at constant field, as in Fig. 
I5.5|l . or by increasing S up to a given value i?„ at constant T. This is indeed 
observed experimentally, as shown in Fig. 15.61 The effect of the field is just to 
polarize the electron spins, which reduces the fluctuations of the hyperfine field, 
thus effectively disconnecting nuclear spins from the lattice. It also follows that, 
for a given Te, is linear in T„ as long as the field-dependence of rn can be 
neglected as compared to the exponential dependence of Tj~^ on T and B. The 
inset of Fig. 15.61 shows Bn{Tn) obtained either from T-sweeps at constant B (as 
in Fig. 15.5(1 or from i?-sweeps at constant T (as in Fig. 15.6(1 . The two methods 
prove to be indeed perfectly consistent with each other. The slope of the linear fit 
of Bn{Tn) gives an average value of tq — 3 x 10~^ s in the high-field regime. We 
observe that at lower fields i?„(T„) tends to deviate from the linear fit, which we 
attribute to the decrease of th. Using the average value of tq we have calculated 
the time-dependent Cnuci from Eq. ()5.1(1 . which has been added to the calculated 
electronic specific heat to yield the solid lines in Figs. 15.51 and 15.61 and seen to 
be in good agreement with the experimental data. 

5.2.2 s^Mn NMR 

We have now the possibility to compare the NSLR rate as inferred from the 
specific heat experiments, with direct measurements of this quantity by NMR 
experiments. The ^^Mn NSLR and TSSR rates were obtained following the same 
procedures and analysis as described in ^4. II (Fig. 15.7(1 . but using a standard ^He 
system instead of the dilution refrigerator. The advantage of the ^He cryostat is 
that, because of its construction, we do not need to insert a A-cable in the NMR 
resonant circuit (cf. ^2.2.2(1 . This increases the quality factor of the resonator 
and allows to measure large enough signals without need to reach millikelvin 
temperatures, where the very slow NSLR would make the measurements almost 
unfeasible. Unfortunately, the useful frequency range becomes accordingly much 
narrower. 

The inversion recoveries [Fig. I5.7f a')] are best fitted by a stretched exponential 
[Eq. ((4.2(1 ] with a ~ 0.5, although the choice of the stretching exponent does 
not strongly influence the extracted value of W. The field-sweep NMR spectra 
in Fig. 15.81 clearly show that it is impossible to determine whether there are 
inequivalent sites in the molecule, as regards the hyperfine coupling (compare 
with the case of Mni2-ac in Fig. 13.7(1 . This may be due to the large quadrupolar 
splitting expected in Mn'^"^ sites, plus the fact that our sample is an unoriented 
powder. As expected from the internal ferromagnetic structure of the cluster 
electron spins (cf. ij3.2.5(l . the ^^Mn spectrum shifts to higher fields by lowering 
the frequency. The spectra can be fitted by a Gaussian shape with total width 
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Fig. 5.7 : Inversion recovery (a) and decay of transverse magnetization (h) for 
the ^^Mn nuclei at T = 0.9 K, B = 5 T and v = 251.5 MHz. The lines in (a) 
are fits to Eq. \4.1\j (solid), and Eq. i4-^ with a ~ 0.5 (dashed). The solid line 
in (b) is a fit to Eq. \4-^ - 



2aB — 2.2 T; if this width were due to quadrupolar sphtting only, it would imply 
even larger than the highest Al^q ~ 4.3 MHz in the Mn*^^^ 
sites of the less symmetric Mni2-ac cluster. We expect therefore that the random 
orientation of the crystallites and, eventually, the presence of inequivalent Mn 
sites as regards the hyperfine coupling, should contribute as well to the observed 
broadening. Indeed, when decreasing the frequency by 5 MHz the maximum of 
the spectrum shifts only by 0.24 T, instead of the 0.47 T that would be expected 
when all the local hyperfine fields are antiparallel to B. We conclude that the 
observed spectrum, as well as the NSLR and TSSR data, should be considered 
as obtained from a mixture of nuclear signals arising from randomly oriented 
crystallites with largely overlapping and quadrupolar-split NMR lines from all 
the Mn sites in the cluster. Extrapolating to i? = the field-dependence of the 
peak of the spectrum, one would obtain zv(0) ~ 360 MHz ^ Bhyp ~ 34 T, very 
similar to the value found in the Mn^^^ site of Mni2-ac. 

Fig. 15.91 shows the field-dependencies of the NSLR rate W and the TSSR rate 
Tg"^, measured at constant frequency u = 251.5 MHz and temperature T = 0.9 
K. From the discussion above it is clear that these data must be interpreted with 

'^'This estimate is obtained by comparing to the Mn'^' line in Mni2-ac, where Auq = 0.72 
MHz yields 2a^ = 2.4 MHz. 
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Fig. 5.8 : ^^Mn NMR spectra at T = 0.9 K and measuring frequencies v = 256.5 
MHz (black circles) and v = 261.5 MHz (gray squares). The gap in the data 
around B ~ 5.8 T is due to the presence of the ^H line. The lines are Gaussian 
fits with total width 2(7b ~ 2.2 T. 

a certain caution, since shifting B at constant v means that we are samphng each 
time a different portion of the NMR signal, which means different quadrupolar 
satehites, different orientation of the crystahites, etc. Nevertheless, the agreement 
with the estimate of Ti obtained by specific heat data (inset of Fig. 15. 5|) is quite 
satisfactory. We can directly compare the NSLR rate W{B) as obtained by 
NMR with the estimated Ti from specific heat data by plotting the line W{B) = 
l/{2Ti{B)) with Ti{B) calculated from Eq. El fixing tq = 3 x 10"'' s and 
T = 0.9 K (solid line in Fig. 15. 9|) . The agreement is very good on the high-field 
side of the data, whereas at low fields the NSLR rate is seen to level off instead of 
increasing exponentially. This is actually the same phenomenon that we already 
observed in the low- field part of i?„(T„) in the inset of Fig. 15.51 due to the fact 
that To is not constant. We fitted therefore W{B) including the field-dependence 
of Tq: 



As shown by the dashed line in Fig 15.91 this function indeed fits much better the 
low-field part of the data. 

In Fig. 15.91 we also show the TSSR rate T2 (B), rescaled in order to evidence 




(5.4) 
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Fig. 5.9 : Field dependencies of the Mn NSLR (black circles, left scale) 
and TSSR (grey squares, right scale) at T = 0.9 K and v = 251.5 MHz. The 
vertical scales are adjusted to evidence the similar field dependencies. Solid line: 
calculated W{B) according to Eq. \5.'^) with tq = 3 x 10~^ s. Dashed line: fit to 
Eq. JO). 



the practicahy identical 5-dependence as found for the NSLR rate. This is indeed 
to be expected (cf. ^3.5.11 and ^3.5.2|) if the TSSR is induced by the random 
change of local hyperfine field when the electron spin makes thermally-assisted 
transitions between the ground- and the first excited states. Eq. (|3.28j) can be 
immediately readapted for Mng by substituting the energy gap /S.E = gfisB and 
the lifetime of the ground state T12 = C^-ph^E'^ /[{IS.E /kBT) — 1], which yields: 

^2 " = .M9^^BB/kBT) - 1 " ^-p^(^^/^^^) [-T^J' (^-^^ 

i.e. the same dependence as W{B)^^\ 

The similarity between NSLR and TSSR rates produced by thermal fluctua- 
tions of the electron spin was already observed in the T-dependence of W and 
in Mni2-ac for T > 0.8 K (cf. ^4.2.1jl . It is interesting to observe that in 
Mne we find T^'^ /T^^ ^ 3800, i.e. larger than the value T^^/T^'^ 1000 found 
for the Mn^"*" site of Mni2-ac and much larger that T2^/T^^ ~ 200 for its Mn^^ 
sites |95| . This in fact confirms the hypothesis that a complete description of 
the NSLR in the thermally-assisted regime of Mni2-ac should include the effect 

'•^'Notice that in Eq. (15.411 the factor 1 — exp(— g/xsB/fcBT) ~ 1 at the denominator is a 
constant except at very small fields. 
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of thermally-assisted tunneling, yielding a contribution to the NSLR which is 
obviously absent in Mng. 

5.3 Conclusions 

In conclusion, our experiments on Mng demonstrate that dipole-dipolc interac- 
tions between molecular magnetic clusters may indeed induce long-range mag- 
netic order at low temperatures, provided that the anisotropy is sufficiently small. 
The spin-lattice relaxation then becomes fast enough to produce equilibrium con- 
ditions down to the low temperatures needed. We should add that similar condi- 
tions could in principle also be reached in the highly anisotropic cluster systems, 
for which it was shown that by applying magnetic fields perpendicular to the 
anisotropy axis, the spin lattice relaxation can be tuned and made similarly fast 
through the process of magnetic quantum tunneling. However, magnetic ordering 
phenomena have not been seen in those systems so far. This is most probably 
due to the fact that, given the magnitude of the fields needed to have a consid- 
erable increase of the relaxation rate {B± ^ 1 T), any longitudinal component 
of the field would create a Zeeman splitting that is much larger than the energy 
involved in the magnetic dipolar ordering. We found indeed that in Mng the 
ordering transition is removed already for relatively small fields (~ 0.5 T). 

We have also studied the nuclear spin dynamics of Mng , either directly by NMR 
experiments, or through the hyperfine contribution to the field-dependent specific 
heat. The agreement between the two techniques is very good, and provides an 
interesting comparison with the nuclear spin dynamics in the anisotropic SMM 
Mni2-ac. We find that the splitting of the electron spin levels in Mng, as produced 
by the application of an external field, yields similar features in the nuclear spin 
dynamics as the crystal field splitting in Mni2-ac. On the other hand, since 
the quantum tunneling fluctuations are strictly absent in Mng, the quantitative 
comparison of the data confirms that, even in the thermally activated regime, a 
correct description of the ^^Mn NSLR of Mni2-ac should include the effects of 
thermally-assisted tunneling. 
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Samenvatting 



van het proefschrift met de titel: 

"Kwantum spin dynamica in Enkel-Molecuul Magneten" 

Terwijl de fundamentele wetenschap steeds meer bezig is met onderzoek aan 
complexe systemen, blijft de Industrie investeren in de voortdurende miniaturise- 
ring van, bijvoorbeeld, clcktronische circuits en gchcugcns. Deze twee ontwikke- 
lingsrichtingen gaan elkaar ontmoeten op de nanometer-schaal, met gevolg dat 
de nanowetenschap een zeer belangrijke onderzoeksgebied eordt voor zowel de 
fundamentele als de toegepaste fysica. 

Ecn van dc mccst bclangwckkende verschijnsclcn die op de nanometer-schaal 
kunnen worden bestudccrd is de overgang tussen de klassieke- en de kwantum- 
fysica. Natuurkundigen zijn gewend aan de toepassing van de kwantummechanica 
op het atomaire niveau, maar zodra het te beschrijven systeem groter wordt, kun- 
nen sommige aspecten van de kwantummechanica (met name het meet-postulaat) 
tot paradoxalc conchisics Icidcn (dc "Schrodingcr kat" paradox). In dit opzicht 
is de theorie nog niet afgerond, het is dus essentieel om systemen te vinden die 
geschikt zijn om experimenteel te testen welke factoren van belang zijn voor de 
overgang tussen kwantum en klassiek gcdrag. 

Dit proefschrift beschrijft een onderzoek aan het kwantummechanisch gedrag 
van Enkel-Molecuul Magneten (EMMen). Dit zijn chemische verbindingen die 
bestaan uit een centrale cluster van magentische ionen, omringd door organische 
liganden. De magnetische wisselwerking tussen dc clcktronspins van dc ioncn in 
eenzelfde cluster resulteert in een "gigantisch" magnetisch moment (gigantisch 
t.o.v. de typische waarden voor atomaire momenten) voor de cluster als geheel. 
De moleculaire magnetische clusters zijn geordend in een kristalstructuur, waarin 
de moleculaire magnetische momenten (spins) slechts een zwakke onderlinge wis- 
selwerking vcrtonen, zodat het systeem beschouwd mag worden als een ensem- 
ble van identieke nagenoeg onafhakelijke moleculaire spins. Hiermee wordt het 
mogelijk om de bekende experimented technieken voor vaste stof fysica toe te 
passen op de studie van het gedrag van de EMMen. Hoofdstuk II beschrijft de 
eigenschappen en het ontwerp van de experimentele opstellingen die voor onze 
onderzoek zijn gebruikt. 

Hoofdstuk III is besteed aan de theoretische beschrijving van de fysische eigen- 
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schappen van de EMMen. Een essentieel aspect is de magnetische anisotropie, 
d.v.z. de voorkeur van de cluster spin om zich te orienteren in een van de twee 
tegenovergestelde richtingen langs een bepaalde voorkeursas. Om de spin om 
tc klappcn moot er een bepaalde energie worden toegevoerd, bijvoorbeeld via 
thermische excitatie. Als het systeem tot een temperatuur veel lager dan de 
anisotropie energie wordt afgekoeld, dan zijn de cluster spins "bevroren", wat 
suggereert dat het mogelijk zou zijn om EMMen te gebruiken als magnetische 
gchcugcns, op voorwaarde dat er geen andere mechanisme is dat de richting van 
de spins kan inverteren. 

Dat extra mechanisme bestaat feitelijk wel, als een compleet kwantummecha- 
nisch model van de cluster spin wordt toegepast. Omdat twee tegenovergestelde 
richtingen van de spin dezelfde energie hebben, is er con mogelijkheid voor de 
spin om te tunnelen door de barriere, d.v.z. zonder dat extra energie wordt 
toegevoegd. Dit fenomeen bestaat niet in de klassieke fysica, en het is vrij won- 
der lijk om het te kunnen waarnemen in een zo groot systeem als een EMM. De 
Prokof'ev-Stamp theorie van het "spin-bad" is een antwoord op de vraag hoe 
kwantumtiinncling in een EMM kan plaatsvindcn ondanks zijn koppcling met 
de omgeving, die normaal gesprokcn dc mogelijkheid voor tunnelen geheel uit- 
sluit. De totale elektronspin van een EMM is magnetisch gekoppeld aan de kern- 
spins in zijn omgeving, maar deze koppeling is sterk beinvloed door de dynamica 
van de kernspins. Bovendien, zijn dc dipolaire koppeling tussen clusterspins en 
de invloed van roostertrillingen (fononen) ook van belang voor een nauwkeurige 
beschrijving van het kwamtumgedrag van EMMen. 

Een van de meest aantrekkelijke eigenschappen van de EMMen voor de studie 
van kwantummechanische verschijnselelen op moleculaire schaal, is dat het mo- 
gelijk is om de tunnelwaarschijnlijkheid te veranderen door een extern magnetisch 
veld aan te leggen, loodrecht op de anisotropie as. Deze "instelbaarheid" is niet 
alleen handig voor het fundamentele onderzoek, maar geeft ook de hoop om 
coherente kwantummechanische oscillaties van de cluster spin tc kunnen produ- 
ceren, waardoor de EMMen als basis element voor een kwantumcomputer zouden 
kunnen worden gebruikt. 

We hebben de bovcngcnocmdcn aspecten van de interactie tussen EMMen 
en hun omgeving bcstudeerd door middel van experimenten aan de kernspin dy- 
namica van ^^Mn kernen in de Mni2-ac molecuul, beschreven in hoofdstuk IV. De 
Mni2-ac moleculen bevatten clusters van 12 mangaan atomen, magnetisch gekop- 
peld tot een cluster met totale spin S = 10. Dankzij de hoge anisotropiebarriere, 
is het bij temperaturen lager dan T ~ 2 K mogelijk om kernspinresonantie (NMR) 
experiementen uit te voeren zonder dat een extern magnetisch veld aangelegd 
hoeft te worden, d.w.z met behulp van het bijna-statische hyperfijne veld dat door 
de electronspins geproducccrd is. De metingen van de kern spin-rooster relaxatie 
duiden dat de magnetische fluctuaties veroorzaakt worden door kwantumtunnelen 
van de cluster spin, dat zich manifesteert in een temperatuur-onafhakelijk gedrag 
van de relaxatietijden voor T < 0.8 K. Als extra bewijs voor de anwezigheid van 
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kwantumtunnel verschijnselen, vinden we dat er een maximum van de kernspin- 
relaxatie bij veld nul is, wat klopt met de resonantievoorwaarde voor de cluster 
spin. Verder, bewijzen we dat het mogelijk is om "flip-flop" overgangen tussen 
kernspins in verschillende clusters te hcbbcn. Dit betekent dat de kernspinpolari- 
satie via spin-diff'usie over het hele sample vcrdceld kan worden, wat van belang is 
voor het interne evenwicht van het kernspinsysteem. In feite, is de gemeten spin- 
rooster relaxatietijd verrassend snel, in aanmerking nemend hoe klein de kans 
op tunnelen wel is. We hebben daarom voorgesteld dat het kwantumtunnelen in 
een kleine fractie zogenaamde "snelle moleculen" plaatsvindt: dat zijn moleculen 
waar, vanwege een lokale verstoring, de anisotropic barriere veel lager is dan nor- 
maal, zodat de cluster spin een veel hogere tunnelwaarschijnlijkheid heeft. De 
eigen kernspins van de snelle moleculen worden door het kwantTimtunnelen gere- 
laxeerd, waarna de kernspin diflcusie zorgt voor de overdracht van energie naar 
de kernen in de langzame moleculen. Deze interpretatie van het mechanisme 
van kernspinrelaxatie wordt bevestigd door de vergelijking tussen het gedrag van 
de kernspins in verschillende niet-equivalente Mn ionen in het molecuul. Ander- 
zijds, zoals voorspeld door Prokof'ev en Stamp, stiumuleert de dynamica van de 
kernspins het kwantumtunnelen van de cluster spin. We hebben NMR metin- 
gen uitgevoerd op een sample waarin de water stof atomen gedeeltelijk vervangen 
waren door deuterium atomen, waardooe de kernen een kleiner magnetisch mo- 
ment hcbbcn. Indcrdaad vondcn wij daarop ccn vccl langzamcre kcrnspin-roostcr 
en spin-spin relaxatie, waarmce is aangetoond dat de tunnelwaarschijnlijkheid 
kleiner wordt als de elektron-kern koppelingen worden verzwakt. 

Additionele metingen aan de afhankelijkheid van de kern spin-rooster relaxatie 
van een veld aangelegd loodrecht op de voorkcursrichting bevestigen vroegere 
resultaten van soortelijke warmte experimenten in onze groep. 

Ten slotte hebben we een aspect van de kernspindynamica onderzocht, dat 
nooit eerder in beschouwing was genomen, namelijk het verband tussen de kern- 
spintemperatuur en de kwantumtunnel fluctuaties van de cluster spins. Het is in- 
dcrdaad helemaal niet triviaal om te voorspellen wat de evenwichtspolarisatie van 
het kernspinsysteem zal zijn, als de dynamica gedreven wordt door temperatuur- 
onafhankelijke kwantum verschijnselen. Wij vonden tot onze verrassing uit onze 
experimenten dat de kcrnspintemperatuur in evenwicht blijft met de roostertem- 
peratuur. Dit resultaat vraagt om een uitbreiding van de bestaande theorieen 
van het gekoppeld systeem "kwantum spin -|- spin-bad" met een koppeling naar 
het rooster, en is van groot belang voor een nauwkeurigc bcrckening van de de- 
coherentie snclhcid. We sluiten hoodfstuk IV dan ook af met een voorstel voor 
een model dat al deze observaties in beschouwing neemt. 

In hoofdstuk V beschrijven we een aantal experimenten aan een vrij bijzon- 
dere EMM, de moleculaire cluster verbinding Mn^, die clusters bestaande uit 
6 mangaan ionen in een zeer symmetrischc structuur bevat, zozeer dat dc to- 
tale anisotropic van de cluster bijna verwaarloosbaar is. Hier verwachten we dus 
geen kwantummechanisch tunnel verschijnsel, maar wel dat de cluster spins in 
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evenwicht met het rooster blijven tot zeer lage temperaturen. Dit wordt inder- 
daad bevestigd door onze metingen van de ac-susceptibiliteit en de soortelijke 
warmte. Dankzij de snelle relaxatie krijgen de cluster spins de mogelijkheid 
om de thermodynamischc cvcnwichtstoestand te vinden, die bij voldeoende lage 
temperatuur (beneden 0.16 K) correspondeert met een ferromagnetische lange- 
afstands-ordening. Het bijzondere van deze lange afstands ordening is dat zij 
uitsluitend wordt veroorzaakt door dipolaire interacties tussen de cluster spins, 
zoals bevestigd wordt door numcrickc bcrckcningcn m.b.v. Montc-Carlo simu- 
laties. De metingen van de veld-afhankelijkheid van de soortelijke warmte laat 
zien hoe de kernspins buiten thermsiche evenwicht kunnen worden gebracht als 
de cluster elektronspins geprolariseerd zijn door een sterke aangelegde veld. We 
hebbcn dit vcrder kunnen bevestigen en onderzoeken door NMR metingen aan 
de kernspin-rooster relaxatie van ditzelfde systeem. 

Samenvattend kunnen we stellen dat het proefschrift een grondig onderzoek 
beschrijft naar de kwantumdynamica van moleculaire spins, in het bijzonder naar 
de interactie tussen dc kwantumspin en zijn omgeving. Dit is van groot fun- 
damenteel belang voor dc kwantummechanische beschrijving van mesoscopische 
systemen, en zou zelfs essentieel kunnen zijn voor toekomstige toepassingen van 
de EMMen op het gebied van magnetische geheugens of kwantumcomputers. 
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